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ABSTRACT 

A  procedure  to  approximate  data  given  at  arbitrary  intervals  in 
two-,  three-  and  four  dimensions  using  polynominal  expressions  is 
described.   The  programming  of  the  problem  is  explained  in  each  case 
and  a  user  manual  is  given  for  software  implemented  on  the  TPL 
Hewlett-Packard  computer  system.   The  method,  which  is  general,  was 
derived  and  has  been  applied  to  represent  the  calibration  of  flow 
probes  which  have  multiple  sensors. 


1.   INTRODUCTION 

This  report  describes  algorithms  to  approximate  data  patterns  such 
as  those  which  occur,  for  instance,  when  pneumatic-velocity  probes  are 
calibrated.   The  functional  value  (Y)  can  depend  on  either  one  (XI), 
two  (XI  and  X2)  or  even  three  (XI,  X2  and  X3)  parameters.   The  approxi- 
mation Y  =  f(Xl)  will  be  referred  to  as  the  two  dimensional,  the  approxi- 
mation Y  =  f (XI,  X2)  as  the  three  dimensional  and  the  approximation 
Y  =  f(Xl,  X2,  X3)  as  the  four  dimensional  approximation.   These  three 
options  meet  the  requirements  for  probes  used  in  the  Turbopropulsion 
Laboratory. 

Using  the  least-squares  criterion  to  obtain  the  coefficients  in 
assumed  polynominals  leads  to  a  system  of  linear  equations,  which  are,  in 
principle  easily  solved.   Numerical  problems  may  arise  however,  since 
the  software,  as  described  herein  and  as  implemented  in  the  NPS 
Turbopropulsion  Laboratory  Hewlett  Packard  21MX  computer,  uses  32-bit 
real  constants. 

The  author  wishes  to  express  his  thanks  to  Professor  Ray  Shreeve, 
whose  constant  and  critical  interest  was  a  crucial  help  to  solve  this 
mathematical  problem  and  link  it  to  engineering  application. 


2.   TWO  DIMENSIONAL  APPROXIMATION 


2.1.  Problem: 


4 


Xi 


Xi 


A  data  set  of  NPNTS1  data  points  is  given,  where  Y  depends  on 
parameter  XI.   The  data  pattern  is  to  be  approximated  by  a  function 
Y  =  f(Xl),  so  that  the  error  between  data  points  and  analytically 
determined  points  is  lowest. 


2.2.  Approach: 

Polynorainals  are  commonly  used  in  order  to  approximate  data 
patterns.   Y  =  f(Xl)  is  a  function  which  approximates  the  data  value, 
V,  at  each  value  of  the  variable  XI,  we  look  for  our  expression  of 
the  form 

2  (L-i) 

Y  =  C   +  C  *Xi  +  C  *Xi   +  ...  +  C  *Xi 
i     2        3  L 


or 

Y  =  EC  *Xi<W) 

i-i  i  X1      .  (2.1.) 


in  which  the  coefficients  C.  are  to  be  determined  by  the  method  of 

1 

least  squares.   As  Ref .  1  shows  the  least  squares  criterion  leads  to 
a  linear  equation  system  that  can  easily  be  solved.   We  define  the  error 


NPNTSi  2 

R  =      E  Cf<Xi  )-Y  ] 
r  =  i      r    r 

(2.2.) 

where  the  index  r  denotes  the  individual  data  point.   Using  equation 
(2.1. )  R  becomes 


NPNTSi  2  (L-i)        2 

R  =      £  CC  +  C  *Xi   +  C  *Xi   +  ...  +  C  *Xi       -  Y  1 
r  =  i   i    2    r  2    r  L    r         r 


2.3.   Solution: 


R  depends  on  the  selection  of  the  coefficients  C  , . . . ,C  . 
In  order  to  minimize  the  error,  R  is  differentiated  with  respect 

to  C  ,  ...,C   and  the  partial  derivatives  are  set  to  zero.   Thus 

L  1j 


3  R        ! 

=    °  i    =    i  ! 

i 


or 

3R  NPNTSi  2  (L-i) 

=  E  2*CC    +C    *Xi    +C    *Xi    +     ,.,     +    C    Xi  -Y    3* 

2C  r  =  i  12r3r  lr  r 

i 

2                                 (L-i) 
^  [C    +C    *Xi    +C    *Xi    +     .  .  .    +    C    Xi               -Y    ] 
i       2         r       3         r                         Lr                 r 
* =    o 

^   C 

i 

Assuming  that  the  summations  extend  over  all  the  data  points,  Z 

NPNTSi 

should  be  understood  to  mean     22        •      Performing  the  differen- 

r  =  i 

tiations  and  rearranging  the  equations,  we  get 


NPNTSi*C  +£Xi    *C       +  ZXi2*C       +  ...    +£Xi(L~i)*C       =  £Y 
1                 r2                 r3  Lr 

£xi  *c  +27x^*0    +  £xi  *c    +  ...  +  £xiL        *c    =Ey  *xi 

ri  r2  r3  r  lrr 

Cxi  *c  +  £xi  *c    +27xi  *c    +  ...  +  Z7xi(L+i>*c    =£y  *xi2 

ri  r2  r3  r  Lrr 


ZTxi'^c    .27xiL»c    +  2Tx1<L*1)*c  +  .  .  +  £xi<2L-2>»c  =£Jy  ixi^1 


r         r 


In  matrix  notation 


NPNTSi 


CXi 

rxi 


CXi 
CXi 
CXi 


L-l  L 

CXi  CXi 

r  r 


(^ 


CXi 

CXi 

CXi 

I 

! 
CXi 


L  +  i 


Exi 
CXi 

CXi 
CXi 


L-i  y 

r 
L 
P 
L  +  i 


01  _o 


C 
^        l' 


Ey 


EY    *Xi 

P  p 


EY    *Xi 

p  p 


L-i 
EY    *Xi 
\      r        p 


*      C:         =    X*         (2.3.) 


c:: 
x* 


Systen  Matrix 
Coefficients  vector 
Right  hand  side  vector 


2.4.   Equation  System 
2.4.1.   System  Matrix  A 

Since  all  Matrix  elements  on  diagonal  lines  running  up  from 
left  to  right  (2.  order,  as  defined  in  Appendix  A)  are  identical,  the 
elements  are  renamed  for  simplification  as  follows: 


*L 

-v 

L-l}2 

^  L-2jS 

ai« 

%z 

"  aL-t,J 

fw 

'-«LA 

a 


21-1 


a 


2;L-i 


a 


',L 


**,*-!   -      *2jL 


°f.    L-l  '     & 


'-    Q. 


so  that  the  set  of  elements,  a  ,  can  be  written  as 


H**nz> 


ak  -  ZV/ 

r-\       r 


k-i 


*./,...,££-/ 


(2.4.) 


This  operation  is  programmed  in  REAL  FUNCTION  S2  (NPNTS1 , IP0WR1.IY) . 
The  data  XI  and  Y  are  known  to  this  function  through  a  common  block 
named  DTA2.   So  only  NPNTSI,  IP0WR1  (=  k-1)  and  IY  (=0)  have  to  be 
passed  to  the  function,  and  the  value  of  a,  is  returned  through  the 


function  name  S2. 


SUBROUTINE  MAT2  presets  the  system  matrix  in  the  following  way: 


i)   Preset  edge  section  elements  (using  REAL  FUNCTION  S2) 


CL.  <=$  a,   =0  a,  <=$ 


a 


L-n 


& 


a 


2L-I 


7op  odor   secho-r> 

/tV^>/    /?cu-</  ec/ly?  scrAo-** 


ii)      Copy   defined   elements   diagonally 


On      *   *>Z       />a<2 


"it 
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&,,  °>1 
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^22       y?On 


tr 


Oil 
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sl. 
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10 
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73 
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22 
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mm  •   •   •  *   •    #  >\    *    *    * 
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tf. 


£4 
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LL 


2.4.2.   Right  hand  side  vector  B 


The  elements  of  the  right  hand  side  vector  B  can  be  written  as 

4*"    EfX-x.^        k-  I  -A 


To  calculate  b,  ,  again  REAL  FUNCTION  S2  is  used.   Data  XI  and  Y  are 
k 

available  through  COMMON  block  DTA2.   NPNTS1,  IP0WR1  (=k-l)  and  IY  (=1) 
have  to  be  passed  to  the  function,  and  the  value  of  b,  is  returned 
through  the  function  name  S2. 
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2.5.   Software: 

The  software  to  compute  the  coefficients  for  a  two  dimensional 
approximation  is  described  in  Appendix  B  and  is  implemented  in  the 
TPL  HP-21MX  computer  system 

To  work  correctly  with  these  program  modules,  the  user  must 

conform  to  the  following  conventions: 

i)   Provide  the  data  in  two  arrays  (Type:  REAL)  of  256 
elements  through  a  COMMON  block,  named  DTA2 . 

COMMON  /  DTA2  /  XI, Y 
REAL  X1(256),Y(256) 

ii)   Dimension  an  array  (Type:  REAL)  of  7  (seven)  elements 
to  contain  the  coefficients. 

REAL  C0EF(7) 

iii)   Define  the  parameters  NPNTSl,  L  and  IPRINT  (Type:  all 
INTEGER) 

NPNTSl    ...   //of  data  points 
L  $  NPNTSl  £  256 

L        ...   (desired  order  of  polynominal)  +  1 

1  £  L  S  7 

IPRINT    . . .   controls  quantity  of  print  out 

2  ...   Print  system  matrix  and  right 

hand  side  vector  before  and 

after  Gauss  Jordan  Elimination 
1    ...   Print  equation  system  after 

Gauss  Jordan  Elimination 
<0   ...   No  print  out 
>0   ...   Print  equation  (1,1)  with  the 

actual  parameters 

iv)   When  loading  a  program,  that  uses  the  subroutine  MAT2, 

the  binary  library  file  has  to  be  searched  for  externals, 

Suppose,  the  source  file  &USER::26  contains  a  user 
program  named  USER.   This  program  calls  MAT2.   After 
the  compilation,  the  relocatable  binary  file  for  this 
user  program  is  %USER::26.   To  load  the  program,  the 
following  procedure  is  recommended: 


12 


Type  (from  FMGR) 

:RU,LOADR 

and  the  loader  program  will  respond 

/ LOADR :   RE,%USER: :26 

where  the  underlined  information  already  is  the  user's 
input.   This  causes  the  loader  to  load  all  program 
modules  of  the  user  program;   a  load  map  is  listed  on 
the  terminal.   Upon  completion  the  loader  prompts 

/LOADR:   MS,%TPLBL: :26 

where  the  underlined  information  already  is  the  user's 
input.   Now  the  loader  conducts  a  search  for  all  unsatis- 
fied externals  of  the  user  program  in  the  binary  lib- 
rary file  %TPLBL::26.   Since  some  library  programs 
have  externals  themselves  the  search  has  to  be  repeated 
(MS  . . .  multiple  search)  as  many  time  as  is  necessary 
for  all  externals  to  be  satisfied.   The  loader  prompts 

/LOADR:   END 

where  the  underlined  information  is  the  user's  input. 
The  loader  now  loads  all  system  programs,  outputs  a 
load  map  and  generates  the  program.   Upon  completion, 
the  loader  outputs  a  ready  message 

/LOADR:USER   READY  AT  1:28  PM  WED.,  10  SEPT,  1980 
/LOADR:  END 

Now  the  program  can  be  run 
:RU,USER 
Again,  the  underlined  information  is  the  user's  input. 
If  all  these  requirements  are  met,  the  correct  call  for  the  sub- 
routine is: 

CALL  MAT2  (NPNTS1 ,L,C0EF, IPRINT) 
Upon  completion,  the  array  COEF  contains  the  coefficients  (C0EF(1)=C. 
, . ..,C0EF(L)=C  ).   Externals  used  by  MAT2  are:   AB2,  DTA2,  S2  under 
no  circumstances  may  the  user  use  any  of  these  names  for  modules  of 
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his  own  user  program.   In  some  cases  the  program  may  not  be  able  to 

perform  a  Gauss  Jordan  Elimination  to  the  system  matrix  and  the  right 

hand  side  vector  and  thus  cause  the  program  to  stop.   If  this  happens, 

an  error  message  is  displayed. 

It  is  highly  encouraged,  to  use  the  system  function  FNP  to 

calculate  the  value  of  a  polynominal  at  one  specified  XI  rather  than  a 

loop  such  as  the  following: 

S=.0 

DO  08  1=1, L,l 
08  S=S+X**(I-1)*C0EF(I) 
YCALC=S 

To  call  FNP  is  easier  and  faster,  because  FNP  is  a  programmed  Horner 

Scheme.   The  recommended  call  is  then; 

L1=L-1 

YCALC=FNP ( COEF , X , Ll ) 
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2.6.   Sample  User  Program 

2.6.1.      Listing 

PAGE  OOOi       FTN.       12:40    PH      THU . ,     18      SEP.,     1980 

0001  FTN4,L 

0002  PROGRAM  DEM02  (3,99) 

0003  C  

0004  C 

0005  C  THIS  IS  A  DEMONSTATION  PROGRAM  AND  IT  SHOWS  THE  CORRECT  USE 

0006  C  OF  THE  TPL  BINARY  LIBRARY  (ON  TYPE  6  FILE  XTPLBL)  USING 

0007  C  .   EXAMPLES.   FOR  FURTHER  QUESTIONS  READ  THE  TPL-LIBRARY-BINDER 

0008  C  .   OR  CONSULT  THE  SOFTWARE  MANAGER  . 

0009  C 

0010  C  

0011  *,  DEMONSTRATE  THE  USE  OF  THE  TPL  BINARY  LIBRARY. 
0012 
0013  COMMON  /  AFLD  /  A 

ooi4  immmmmmm 

0015 

0016  REAL  A(2S6) 

0017 

0018 

0019 

0020  INTEGER  NOLF ,NOCR <2> , ICLR <3> 

0021 

0022  DATA  NOLF   /006537B/ 

0023  DATA  NOCR   /00 0033B, 040433B/ 

124  DATA  ICLR   /0 1S524B ' 0 15515B , 0 06S37B/ 

„02S  DATA  PI     /3.141S93/ 
0026 

0  027  C     I"  O  R  M  <**  T  S5   I>  IE!  M  O  21       S3  T  r^  R  T 

0028  101  FORMAT  (///"   HELLO!  THIS  IS  PROGRAM  DEMO  AND  WE'LL  SEE  H 

0029  *OW  TO  USE  THE  MARVELLOUS  TPL  BINARY"A2/"   LIBRARY!   YOUR  INTERES" 

0030  *WILL  BE  GREATLY  AWARDED  BY  EASIER  PROGRAMMING."///) 

0031  102  FORMAT  <"   GENERATING  DATA  P0INTS"/5X"I "SX"X( I ) "5X"Y ( I) " ) 

0032  103  FORMAT  <2X , 14 >2F9 . 3, A2) 

0033  104  FORMAT  <"   DATA  POINTS  GENERATED!"/) 

0034  105  FORMAT  ("   INITIALIZING  THE  PLOTTER") 

0035  106  FORMAT  (■   PLOTTER  INITIALIZED!"/) 

0036  107  FORMAT  ("   DEFINE  PLOTTER  AND  USER,  AREAS  AND  DRAW  AXES") 

0037  108  FORMAT  (9X""6X"XMIN  "6XUXMAX  "6X"YMIN 

0038  *"6X"YMAX"/9X,4(""10X"  ")/ 

0  039  *"   HP9872",4("  "F10 . 3" "  )  ,  "  "/ 

0  0  40        *"   USER 4  (  " — "{-  i  o — j  »  "  )     "  "  /  ) 

0041  109  FORMAT  ( "  'PLOTTER  AREAS 'DEFINED  AND  AXES  DRAWN!"/) 

0042  110  FORMAT  < "   DRAWING  DATA  POINTS  INTO  COORDINATES  SYSTEM") 

0043  ill  FORMAT  ("   DATA  POINTS  DRAWN!"/) 

0044  112  FORMAT  ("   CALCULATING  A  CURVE  FIT  THROUGH  THE  DATA  POINTS 
0  045  *") 

0046  113  FORMAT  ("   ENTER  DEGREE  OF  POLYNOMIAL  (NORDER)  TO  FIT  THRO 

0047  *UGH  THE  POINTS   "2A2) 

0048  114  FORMAT  <"   ENTER  IPRINT   "2A2) 
"149  115  FORMAT  (/"   CURVE  FIT  DONE!"/) 

J50  116  FORMAT  ("   CURVE  FIT  DONE!"/) 

0051  117  FORMAT  <"   PLOTTING  THE  CURVE  FIT") 

0052  118  FORMAT  ("   CURVE  FIT  DRAWN!"/) 

0053  149  FORMAT  (H"((3A2))) 

0  054  C      F  O  R  M  ^  T  3   I>  in:  M  (D  ZZ       <:>  T  O  P 

0055 

0056 

0057  C  

0058  C 

0059  C  GET  THE  LU  OF  THE  TERMINAL. 

0060  C 

0061  C  

0062  LI  =  LOGLU(I) 

0063  WRITE  (LI,  101)  NOLF 
0064 

0065 
0066 

0067  C  

0068  C 

0069  C  .   GENERATE  DATA  POINTS  SCATTERED  AROUND  A  POLYNOMIAL. 

0070  C 

0071  C  

0072  CALL    INITG    (13) 

0073  WRITE  (LI,  102) 

0074  XSTART  =  -6. 
0  075  XSTOP   =  +2. 
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18   SEP 


1980 


0076 
0077 
0078 
007? 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
"198 
J99 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
Oil! 

lift 

0114 
0115 
0116 
0117 
0118 
0119 
0120 
0121 
0122 
0123 
124 
ui25 
0126 
0127 
0128 
0129 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 


01 


C 

c 
c 
c 
c 


I 

c 
c 

c 


c 
c 
c 
c 
c 


DX 

NORDER 

COEF(i) 

C0EFC2) 

C0EF(3) 

COEF(4) 

C0EF(5) 

A 

DO  01 

Xi(I) 

Y(I) 

WRITE 

ISTOPi 

DO  02 

DUMMY 

Y(I+i> 

Y(I) 

WRITE 

WRITE 


=  <XSTOP-XSTART)/(NPNTSi-i) 
=  4 

=  -1  . 

=  -0 .500000 

=  -0  .  000340 

=  +0.300000 

=  +0 . 055000 
=  <2.0*PI>*0 . 03125 
1  =  1  .NPNTSivi 
=  XSTART+(I-i)*DX 
=  FNP(CQEF,Xi(I)  , NORDER) 


+  0 .25*SIN(X1(I)/A) 


(LI,  103)  I,Xi(I) ,Y(I),NOLF 

=  NPNTSl-i 
1=1, ISTOPI, 2 

=  Y(I+i) 

=  Y(I) 

=  DUMMY 
(LI,  149) 


(LI,  104) 


(ICLR, 1=1,3,1) 


INITIALIZE  THE  PLOTTER 


WRITE  (LI,  105) 
LP  =  13 

CALL  INITG  (13) 
WRITE  (LI,  149) 
WRITE  (LI,  106) 


ICLR 


DEFINE  PLOTTER  AND  USER  AREAS;   DRAW  AXES 


WRITE 

XPMIN 

XPMAX 

YPMIN 

YPMAX 

XUMIN 

XUMAX 

YUMIN 

YUMAX 

WRITE 

ALPHAX  = 

XA 

XB 

XL 

ALPHAY  = 

YA 

YB 

YL 

CALL 

CALL 


(LI,  107) 


L4 

=  12 

=  12, 
=  22 


(LI 


-6. 
+2. 
-5. 
10. 
,  108) 

0.00 
(XPMAX 


XPMIN , XPMAX , YPMIN , YPMAX , XUMIN , XUMAX , YUMIN , YUM ' 


CALL  AXIS 
WRITE  (LI, 
WRITE  (LI, 


-XPMIN       )/(XUMAX-XUMIN) 
(XPMIN*XUMAX-XPMAX*XUMIN)/(XUMAX-XUMIN) 
(XPMAX       -XPMIN       ) 
90.00 
■  (YPMAX       -YPMIN       )/(YUMAX-YUMIN) 
=  (YPMIN*YUMAX-YPMAX* YUMIN )/(YUMAX-YUM IN) 
=  (YPMAX       -YPMIN       ) 
SETSM  (113,1. ) 
AXIS  (XPMIN, YPMIN, XL, ALPHAX, 2HX1, +2, XUMIN, XUMAX, 4HF4 


i>4 


(XPMIN, YPMIN, YL, ALPHAY  2H  Y, -2; YUMIN, YUMAX , 2HI2,2 , 3 
149)  (ICLR, 1=1, 6,1) 
109) 


PLOT  DATA  POINTS 


WRITE  (LI,  110) 
CALL  SETSM  (113,3. ) 
DO  03  I=l,NPNTSl,i 
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XPLOT   =  Xi(I)*XA+XB 
YPLOT   =  Y(I)*YA+YB 
CALL  PLOT  (XPLOT, YPLOT, 2) 
03  CALL  SYMBL  (1) 

WRITE  (LI,  149)  ICLR 
WRITE  (LI,  ill) 


C      

C 

C      .   CALCULATE  CURVE  FIT  THROUGH  DATA  POINTS. 

C 

C 


WRITE  (LI,  112) 
04  WRITE  (LI,  113)  NOCR 
READ   (LI,    *)  NORDER 
WRITE  (LI,  149)  ICLR 
IF  (  NORDER   .LT.   0    .OR.    NORDER   . GT .   6  )  GO  TO  04 


WRITE  (LI,  114)  NOCR 

WRITE  (LI',  149)  (ICLR, 1  =  1, 2,1) 


C      

C 

C         PLOT  POLYNOMIAL  FIT. 

C 

C 


WRITE  (LI,  117) 
CALL  SETSM  (113,2. ) 
DO  09  I=l,NPNTSi,i 
XPLOT   =  Xi(I)*XA+XB 

YPLOT   ■  FNP(COEF,Xi(I)  .NORDER )*YA+YB 
IF  (I.EQ.i)  CALL  PLOT  (XPLOT , YPLOT, 1 ) 
09  IF  (I.GT.i)  CALL  PLOT  (XPLOT , YPLOT ,3) 
CALL  SETSM  (1 IS, DUMMY) 
WRITE  (LI,  149)  ICLR 
WRITE  (LI,  118) 


C      

C 

C      .   TERMINATE  GRAPHICS 

C 

C      


CALL  STOPG 
STOP  7777 
END 


FTN4  COMPILER:  HP92060-16092  REV.  1926  (790430) 

**   NO  WARNINGS  **   NO  ERRORS  **    PROGRAM  =  01409       COMMON  =  00000 
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2.6.2.      Load  map 


0EMG2 
FNP 


10042  12642 
12643  12746 


DEMONSTRATE  THE  USE  OF  THE  TPL  BINARY  LIBRARY 

Calculate  n-th  order  polyn 


to 

load    th 

ese   programs,    type    (from  LOADR) :    MS , %TPLBL 

mam. 

12747 
14747 
17000 
17132 

14746         2D-Approxj.«ation       /    DTA2      /                      (DATA2         DTA2) 
16777         2D-Appr oxinat ion    arrange    system    Matrix    and    vector 
17131         2D-Appr oxinat ion                                       cowpute    supinations 
17311         2D-Approxination       /    AB2         / 

AFLD 

AXIS 

INITG 

PLOT 

SETSh 

STOPG 

SYMBL 

FCTR 

FACTR 

LOGLU 
READF 
OPEN 
CLOSE 
CLRIO 
QVRD. 
$SMVE 
LURQ 
.DADS 
.DMP 
.DDI 
SESSN 
R/W* 
P. PAS 
.DNG 
PAUSE 
SALRN 
FMTIO 
ERR0 
TAN 
ABS 
.SNCS 
.DDE 
.DIN 
.RTOI 
.FPWR 
.SBT 
.  .FCM 
PAU.E 
ER0.E 
.ChRS 
$OPEN 
RUSUB 
RWND* 
FRMTR 
FMT.E 
IB4A2 
REIO 
RMPAR 
PNAME 
LUTRU 
IPUT 
SSETP 
.CFER 
.LBT 


17312 
20312 
22632 
23645 
24376 
2S337 
2S402 
26137 
26143 

26163 
26241 
27203 
27501 
27711 
27720 
27721 
30010 
30373 
30503 
30651 
31152 
31170 
31327 
31356 
31366 
31467 
31605 
33104 
33174 
33301 
33310 
33452 
33464 
33472 
33566 
33630 
33671 
33706 
33707 
33710 
33774 
34151 
34517 
34642 
40300 
40301 
41602 
41727 
41772 
42040 
42147 
42170 
42215 
42273 


20311 
22631 
23644 
24375 
25336 
25401 
26136 
26142 
26162 

26240 
27202 
27500 
27710 
27717 
27720 
30007 
30372 
30502 
30650 
31151 
31167 
31326 
31355 
31365 
31466 
31604 
33103 
33173 
33300 
33307 
33451 
33463 
33471 
33565 
33627 
33670 
33705 
33706 
33707 
33773 
34150 
34516 
34641 
40277 
40300 
41601 
41726 
41771 
42037 
42146 
42167 
42214 
42272 
42323 


GSP 
GSP 
GSP 
GSP 
GSP 
GSP 
GSP 
GSP 
GSP 


AFLD 
read 


/  A(256) 
control 


array 

no  ye 


/  FCTR   / 


Print 
Y 


Draw  and  label  axes 
Aj   initialize  plotter 
pen  to  a  defined  point 
change  plotter  nodes 
terminate  grafics 
synbol  at  location  of  pen 


wary  size  of  the  plot 


92067-16268  REV. 1903 

92067-16125  REV. 1940 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

750701   24998-16001 

92067-16125  REV. 1903 

92067-16268  REV. 1903 

92067-16268  REV. 1903 

780818   24998-16001 

780818   24998-16001 

781021   24998-16001 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

780818   24998-16001 

771122   24998-16001 

92067-16268  REV. 1903 

24998-16002  REV . 192 

771122   24998-16001 

780424   24998-16001 

750701   24998-1600 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

24998-16001 

REV. 1903 

REV. 1903 


790228 
790719 
790215 
781229 

780526 
790202 
790223 


780413 
781214 
740801 


770715 
6  790417 


780424 

780818 

780818 

780921 

781106 

770518 
750701 
750701 
750701 

780424 
92067-16125 
92067-16125 


92067-16125  REV. 1903 
24998-16002  REV. 192 
24998-16002  REV. 190 
59310-1X013  REV  1940 
92067-16268  REV. 1903 
781106   24998-16001 

771121   24998-16001 

92067-16268  REV. 1903 

92067-16125  REV. 1903 

781106   24998-16001 

750701   24998-16001 

770518   24998-16001 


790103 

781003 

780801 
6  790503 
1  781107 

790802  1153 

790316 


790223 
740801 


15  PAGES  RELOCATED     15 
LINKS:BP     PROGRAM:LB 
/LOADR :DEM02   READY  AT 

/LOADR :*END 


PAGES  REQ'D     NO  PAGES  EMA 
LOAD:TE    COMMON :NC 
12:44  PM   THU.,  18   SEPT,  1980 


NO  PAGES  MSE 
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2.6.3.      Results 


Affect   of    IPRINT   on  quantity   of   print   out 
i)      printed   output 
System   natrix    A    and    vector    B    before    Gauss    Jordan    Elination 


^ 


i 
2 
3 

4 

5 


510E+02 
102E+03 
487E+03 
211E+04 
104E+05 
i 


102E+03 
487E+03 
211E+04 
104E+0S 
S25E+0S 


487E+03 
211E+04 
104E+05 
525E+05 
274E+06 
3 


211E+04 
104E+05 
S25E+0S 
274E+06 
146E+07 
4 


5 
104E+05 
52SE+Q5 
274E+06 
146E+07 
787E+07 

5 


.578E+02 
.998E+02 
. 113E+03 
.944E+03 
.  113E+05 


Equation  system  after  Gauss  Jordan  Elimination 


i 
2 
3 

4 

5 


i 
100E+01 
OOOE+00 
000E+00 
OOOE+00 
OOOE+00 

i 


2 
000E+00 
iOOE+Oi 
OOOE+00 
OOOE+00 
000E+00 


3 
OOOE+00 
OOOE+00 
iOOE+Oi 
OOOE+00 
000E+00 

3 


4 
000E+00 
OOOE+00 
OOOE+00 
iOOE+Oi 
OOOE+00 

4 


5 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 
iOOE+Oi 

5 


.970E+00 
.416E+00 
.162E-01 
.27SE+00 
.514E-01 


>    IPRINT  -2 


Coefficients  COEF(I) 


=  1  - 

=  2  - 

=  3  - 

=  4 

=  5 


970E+00 
416E+00 
162E-Q1 
275E+00 
514E-01 


J 


Equation  system  after  Gauss  Jordan  Elimination 


i 
2 
3 
4 


iOOE+Oi 
OOOE+00 
OOOE+00 
OOOE+00 


OOOE+00 
iOOE+Oi 
OOOE+00 
OOOE+00 


3 
OOOE+00 
OOOE+00 
iOOE+Oi 
OOOE+00 


4 
OOOE+00 
OOOE+00 
OOOE+00 
iOOE+Oi 


5 
OOOE+00 
OOOE+00 
000E+00 
OOOE+00 


.OOOE+00   .OOOE+00   .OOOE+00   .OOOE+00   .iOOE+Oi 
12  3  4 


^ 


970E+00 
416E+00 
162E-01 
275E+00 
514E-01 


>         IPRINT-  1 


Coefficients  COEF(I) 


=  1  - 

=  2  - 

=  3  - 

=  4 

=  5 


970E+00 
416E+00 
162E-01 
275E+00 
514E-01 


J 


Coefficients  COEF(I) 


=  1    - 

=  2    - 

=  3    - 

=  4 

=  S 


970E+00 
416E+00 
162E-01 
27SE+00 
514E-01 


IPRINT  <  I 
>Z 
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ii)   graphic  output 


L6.  0     -4.  0     -2.  0  ■      0.0 


A  j. 


2.  a 


Grafic  output  from  PROGRAM  DEM02 
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3.   THREE  DIMENSIONAL  APPROXIMATION 


3.1.   Problem: 


A  data  set  of  NPNTS1*NPNTS2  data  points  is  given,  where  Y 
depends  on  the  parameters  XI  and  X2.   The  data  pattern  is  to  be 
approximated  by  a  function  Y=  f (XI,  X2) ,  so  that  the  error  between 
data  points  and  analytically  determined  points  is  lowest. 
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3.2.   Approach: 

As  for  the  two  dimensional  approximation,  a  polynominal,  now 
with  two  independent  variables  XI  and  X2 ,  is  used. 


y-  ffxi,  x2 ) 

*  CCt,  ->  C3ZX2+  CftXZ2*  ...    <c3M  x2C"~ '  1-x/  2  * 


^ 


Cm->)^  fL-i 


t  CCU  ,CUXZ<  C^XZ  *  ...   <  CLH  XZ        1-X  I 


The  data  arrangement  for  the  three  dimensional  approximation 
already  yields  to  the  application  of  this  method:   at  each  constant 
Xl  (i.e.:  r  =  constant,  too)  there  is  a  data  set  for  various  X2. 
This  is  the  case  for  example,  when  a  pneumatic-velocity  probe  is 
calibrated.   At  each  Mach  number,  the  probe  is  balanced  in  yaw  and  the 
pitch  angle  is  varied  over  the  range  of  application.   The  probe  out- 
put (each  of  two  separate  non-dimensional  pressure  differences)  is  a 
function  of  the  Mach  number  and  the  pitch  angle. 
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The  least  squares  criterion  becomes 

r*  /    r-  / 

where  the  indices  r  and  s  denote  the  individual  data  points.   Usinj 
equation  (2.1.)  R  becomes 


r-'     *■'  "-/•<  w 


H-l  I-' 


<c„ * cu xzt/  ...iC^XZ^) **/„ 


-  Yn  Jzj 


The  term  in  the     |  -  bracket  shall  be  called  B, 
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3.3.   Solution: 


To  minimize  the  error,  R  is  partially  differentiated  to  the 
coefficients  C  .  and  then  the  partial  derivatives  are  set  to  zero 


SET  ~°  ''•<■..>;  </'*'>->" 


Assuming,  that  the  summations  extend  over  all  data  points,  II 

nphtsi      _  t*m% 
should  be  understood  to  mean    '   . j     J^        .   Performing  the  dif feren- 

/"*.  /      i-t 

tiations  and  rearranging  the  equations  in  matrix  notation  (which, 

while  lengthy,  is  logically  straight  forward),  we  get,  in  matrix 
notation, 
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/< 


(2-l)'/1 


(a-0*M 


I-  ( 


L*h 


Nphtsi  f  HfHTSZ    £  H  XZ . 


ETxz 

ztxz' 


ee«; 

ts.xz':, 


MM  'Zi)fMil 

Zxi„        ZT.xi.xz       TTxi-xz1 


T\Txz„       SE/^"     2TEW 


LTxl 


TExi-xz^ 


TTxi-n       ZTxi-ja' 


TT\xi-xl~'     T.Txi-xz" 

rs      n  rs     TZ 


EEjft/AzJ 


EE 


«„•*&., 


ITx^xi      TTx,ln.X2? 


ee; 


<  .*/  •« 


LT^XhXz',. 

TTxi-xz" 


E&f/ "-A2 "*    E  E*/'~'*2  M 


EL^V""2  E2^>„     SZ//>4 


TLx'n-xzl 

TT/i-xi.' 

rs       n 

TT\*,l-xt 


W-X&-&K 

TTxL-xi 

TLxiL-xz„ 


L 
\_ 


(L-l)*  MH 

ci>;"  sex£«m  :.; 
cs/G««   szaC^z'.  ... 

SSx£w"     TTxfrxC    ... 

L*M 

TZxC-xi:  ■ 
ZTx£.xz": 

TTx£         SZ^-«H    ...              2E<-«;; 

,- 



EZx£  ZZ<"*zB  ... 
22<V,  EEy/^< ... 

SEX".*C    EE<V".  - 

„„        21-1         IH 

ZLX/a  -A2„ 

szx/r-< 
ce<vc 

Q, 


EEv„ 

EE  \*ia 

EE  );<" 

EEY„ 

•XI 

EE  r..*^ 

;*'* 

ELY-xC 

*»„ 

cey0  .  <' 

CEYKiV2B  -< 


ESY^-xC 


*    <C 


<v..?v> 


3.4.   Equation  System: 
3.4.1.   System  Matrix  A 


This  system  matrix  needs  some  examination.   A  crucial  step 
is  the  introduction  of  submatrices  1.  class  into  the  system  matrix. 


/A 


A, 


/A 


12 


A 


IL 


/At         *« 


4L 


A 


u 


A 


LZ 


A 


O- 


The  asterisk  *  denotes  a  submatrix  1.  class  where 


A* 


NP/V7J/  *  NPHTSZ  L-iLiXZ 


ILxi 


rr 


2 


LZx£ 


f*-/ 


LL*z 


£n~Z 


'Z 


'n 


ESXLX2 


a'     a 


*-/ 


ZT>xlx2 


rs    "n 


IS  »b«b 


L2>i,wr    IlJya£ 


»-/ 


EE^, 


x«-z 


EE>7„ 

EEr/,*^ 

a  = 

ZEx/X 

ht.K~*K 

iu/X" 

ex  it*: 

M-l 


E  E  <A2, 

EE  tftf/ 


E.EW 


etc 
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Now  it  can  seem  that  all  submatrices  on  the  same  diagonal  bands 

*  * 

2.  order  are  identical.   The  submatrices  A.,  subsequently  are  renamed  A, 


*\  -  K  -  a*    •  Ks  -  K 


1J 


A*  -  A*,   .  A 


A 


LAI 

/•  LA 


&z 


I 


L-l  I 


--   A 


'J1 


-A 
-  A* 


I -2  J 
J 


■-  A 
Z  H 

~-A]    L-f 


A 
n3,L 


* 

where    A, 


can  be  written  as 


-A 


*A 


A 


£  Ex/  J1  xz° 


En    fc-/    J 


rs 


*-»    Af 


k-t 


M-J 


Like  the  system  matrix  for  the  two  dimensional  approximation  this 
submatrix  1.  class  is  not  only  symmetrical  to  its  main  diagonal  line 
1.  order,  but  also  the  elements  on  each  diagonal  line  2.  order  are 
identical.   The  elements  a,  g . .  of  the  submatrix  A.  therefore  are 

K.  J1J  K. 
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renamed  to  a  #  ,  where  k  specifies  the  place  of  the  submatrix  in  the 

K.  m   i-J 

system  matrix  and  i  and  j  define  the  place  of  a     in  the  submatrix  A,  . 

k. ;  i  j  tc 


k;#  =  ak;„ 


[*/3    "k77t  ^k)ZZ 


akjM   "    %»,      '-    aK;H-lZ    '   Q*]M-Z3 


'k-JHl 


%*ifi*L 


ak 


:M2,        -  a*;M-t     3 


'-  <*K  ,  *1* 


ex 


«;ZM-/ 


a 


*;Mfl 


where  a,  _  can  be  written  as 
k;L 


x»iri<   wrwrSZ      Li  i    i 

ak.t   ■  DZ>„'  X2„? 


(3.4.) 


Equation  (3.4.)  therefore  is  the  fundamental  step  towards  the  pro- 
gramming of  the  system  matrix  and  is  programmed  in  REAL  FUNCTION  S3 
Since  the  data  XI,  X2  and  Y  are  already  provided  through  a  COMMON 
block  named  DATA3,  only  NPNTS1,  NPNTS2,  IP0WR1  (=  k-1)  and  IP0WR2 
(-  L-l)  have  to  be  passed  to  the  function  that  returns  the  value  of 
a,    through  the  function  name  S3. 
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SUBROUTINE  MAT31  presets  a  submatrix  A   in  the  following  way 


i)   Preset  edge  section  elements  (using  REAL  FUNCTION  S3) 


4 


•  •  • 

4 


^/in-i 


<**,,-*  «*,M 


7b/o  <*xor     &orf/Oii 
h<Lj  *«**(  <rci<je  set+iori 


ii)   Copy  defined  elements  diagonally 


a««  ^<?*a.     -a< 


a, 


H'.im 


XV  y/^ 


«*, 


aKjH       %iz      4k; /j 


ax^/        Qir#z    ~<*K}U 
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IN 


a, 


A;/Vf 


a*.,     a,y/2     ^z, 

'V'         a^l        4fcj2j 


77Z*"/: 


*;fcv> 
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H;in 


Gkiii         Qkfl        a«}33        •  -  ■  Gtyt 


Qh;HI        Qk,Kl        ^)h3        -.  Cfk; 


nn 


This  subroutine  that  needs  NPNTS1,  NPNTS2,  M  and  IPOWR1  (=  k-1)  as 
input  parameters,  returns  A   (SUBM1  (7  ,  7))  to  SUBROUTINE  MAT3,  which 
presets  the  entire  system  matrix  in  the  following  way 

i)   Preset  edge  section  submatrices  (using  SUBROUTINE  MAT31) 


A*  =M* 


=*>... 


A* 

5 


A 


2L-) 


7^o      pot**    Qpcho-v* 


LH      ~C'    Ait--I 


7 


a/  &>fbf    sarhot^r 


31 


ii)   Copy  defined  submatrices  diagonally 
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A 


A 
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The  parameters  to  call  this  subroutine  are  NPNTS1  and  NPNTS2  (#  of 
data  points),  L  and  M  (<the  order  of  the  polynominal  fit  +  1>  in  XI 
and  X2  -  direction)  and  IPRINT.   The  coefficients  C.  are  returned. 
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3.4.2.   Right  hand  side  vector  B 


The  right  hand  side  vector 
subvectors  1.  class. 


will  be  divided  up  into  L 


B  - 


*: 


Again,  the  asterisk  (*)  denotes  a  subvector  1.  class.   Now  the 

subvectors  1.  class  B   can  be  written  as 

k 


B. 


ZT,Y„xi 


k-l 
K-l 


ZZra-xCxzMa 


k- !,...,  L 


Finally,  the  elements  of  the  subvectors  1.  class  B,   can  be  written  as 

bk/l  -  zz,Yn.xCxzi: 


k-l,...,  L;    L-i,...,M 


To  compute  b    ,  REAL  FUNCTION  S3  is  used.   Data  XI,  X2  and  Y  are 
available  through  a  COMMON  block  named  DATA3.   NPNTS1,  NPNTS2,  LP0WR1 
(=  k-l),  IP0WR2  (=  L-l)  and  IY  (=1)  have  to  be  passed  to  the  function 

and  the  value  of  b,  T  is  returned  through  the  function  name  S3. 

k;^ 
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Using  FORTRAN  programming  language,  the  allocation  of  b    can  be 
performed  by  two  stacked  DO  loops. 


1=0 

DO  08  L1-1,L,1 
IP0WR1=L1-1 
DO  08  M1=1,M,1 
IP0WR2=M1-1 
1=1+1 
08  B(I)  =  S3(NPNTS1,NPNTS2,IP0WR1,IP0WR2,1) 


3^ 


3.5.   Software: 

The  software  to  compute  the  coefficients  for  a  three  dimensional 

approximation  is  described  in  Appendix  A  and  is  implemented  in  the 

TPL  HP-21MX  computer  system.   To  work  correctly  with  these  program 

modules,  the  user  has  to  conform  to  the  following  conventions: 

i)   Provide  the  data  in  three  arrays  (Type:  REAL)  of 
16  *  16  elements  through  a  COMMON  block  named 
DATA3 . 

COMMON  /  DATA3  /  XI,  X2,  Y 

REAL  XI  (16,16),  X2(16,16),  Y(16,16) 

ii)   Dimension  an  array  (Type:  REAL)  of  7*7 
elements  to  contain  the  coefficients. 

REAL  COEF  (7,7) 

iii)   Define  the  parameters  NPNTS1,  NPNTS2,  L,  M  and 
IPRINT. 


NPNTS1 


NPNTS2 


//  of  XI  variations 
L  £  NPNTS1  £  16 

//  of  X2  variations 
M  $  NPNTS2  £  16 


M 


IPRINT 


(desired  order  of  approximation 
polynominal  w.r.t.  XI)  +  1 

(desired  order  of  approximation 
polynominal  w.r.t.  X2)  +  1 

controls  quantity  of  print  out 

2    ...   Print  system  matrix  and 
right  hand  side  vector 
before  and  after  Gauss 
Jordan  Elimination 

1    ...   Print  equation  system  after 
Gauss  Jordan  Elimination 

<0   ...   No  print  out 

>0   ...    Display  equation  (3.1.) 

with  the  actual  parameters 


35 


iv)   If  a  user  program  uses  subroutine  MAT3,  the  software 
modules  have  to  be  loaded  using  the  procedure  out- 
lined in  section  2.4.  iv. 

If  all  these  requirements  are  met,  the  correct  call  for  the  subroutine 

is: 

CALL  MAT  (NPNTS1,NPNTS2,L,M,C0EF,IPRINT) 

Upon  completed  execution  of  this  approximation  routine  the  array 

COEF  contains  the  coefficients.  Externals  used  by  MAT3  are:   AB3, 

DATA3,  MAT31,  IEL3,  S3.   Under  no  circumstances  may  the  user  use 

any  of  these  names  for  modules  of  his  own  program. 
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3-D.    Sample  User  Program 

3.6.1.   PTNM  Compiler  listing  of  sample  user  program 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
001? 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 


0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
004? 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076 


FTN4,L 

PROGRAM  LAB  (3,99) 
*,Plot  calibration  points  and  calculate  coefficients. 

COMMON  /  AFLD  /  PLOTR 
COMMON  /  DATA3  /  Xi,X2,Y 
COMMON  /  AB3    /  A,B 

REAL  PLOTR (256) 

REAL  Xl(16.16).X2(16,i6))Y(16,i6) 

REAL  A<49,49) ,B(49) 

REAL  YT(16,16)  ,R(16)  ,C0EFF(7,7) 

INTEGER  IDCB(144) ,IFILE(3) ,N0CR(2) ,ICLR(3) ,YTITLE(4) ,ZTITLE<4) 

DATA  PI  /3. 141593/ 

DATA  NOLF  /006537B/ 

DATA  NOCR  /000033B , 040433B/ 

DATA  ICLR  /015S24B,01SS15B,006S37B/ 


i;~  a  r  m 

121  FORMAT  ( 

122  FORMAT  ( 

123  FORMAT  ( 
*.   Veloc 

124  FORMAT  ( 

*  =  <  ,  ) 
* 

*nd  Delta 

125  FORMAT  ( 
*X<  ,  >  } 
* 

*  C  X=X( 

126  FORMAT  ( 
*s   "2A2) 

160  FORMAT  ( 

131  FORMAT  ( 

132  FORMAT  ( 

133  FORMAT  ( 

134  FORMAT  ( 

127  FORMAT  ( 

128  FORMAT  ( 

129  FORMAT  ( 

130  FORMAT  ( 

150  FORMAT  ( 

151  FORMAT  ( 

152  FORMAT  ( 

153  FORMAT  < 

154  FORMAT  ( 
*"   nax , 

155  FORMAT  ( 
*"   wax , 

156  FORMAT  ( 
*x.  #  par 

157  FORMAT  ( 
*x  ,  #  oar 


A  V  3   l„  A  l*       i: 


Enter 
3A2) 

"  Select  i 
ity  X    U2A 
"  Select  1 
>"/ 

2 
C   =  (  . 
"  Select  i 
V 

2 

,  )  3   "2A2) 
"  Enter  *  of 


>  T  A  R  T 

none 


raw  data  file 

Pitch  angle 
Pitch  anqle 


)  J 


,iH", 

*943) 


V/S 


"  Enter  1  to 
/"  Calibrate 
/"  Calibrate 
/"  Calibrate 
/"  Calibrate 
24X"  "9X" 
24X"  "9X" 
24X"  "9X" 
24X"  "9X" 

"   M  1  n  i  M  U  M 

"  MaxiMUM 
"  of  user 
"  Enter  1 
"  Enter 


Veloc  ity 
.  ",1H"," 
Mach  Numbers  and 
data 


print  input 

=  (  ,  )v'/) 

=  (  ,  )  "  / ) 

X=X(  ,  )"/) 

X=X(  ,  )"/) 

"9X"  ") 

"<?X"X") 
"9X"XU) 

values", 3(1X,F9 
values-  ,3(1X ,F9 


"2A2) 

or   2  . 
v/s  Ga««a  and  Beta    C 
",1H","      w/s  Gatnna  a 
Ganwa  and  Beta    C  X= 

v/s  Ganna  and  Delta 

*  of  pitch  angle 

"2A2> 


3)  ) 
3)  ) 


Min 


ar t i t ions 

er   «in 

par  t  it  ions 


Ear 
nt 


158  FORMAT 

101  FORMAT 

102  FORMAT 
107  FORMAT 
105  FORMAT 

%  "  3 

103  FORMAT 

104  FORMAT 
149  FORMAK 

611  FORMAT 

612  FORMAT 

613  FORMAT 

614  FORMAT 

615  FORMAT 

616  FORMAT 

617  FORMAT 

618  FORMAT 

602  FORMAT 

603  FORMAT 


"  Enter 
tit  ions 
"  Enter 
tit  ions 
1X,3(F8. 
"  Enter 
"  Enter 
"  Enter 
Select 


coordinates"3X"Xi"8X"X2"9X"Y 
to  redefine  these  data  "2A2) 
u    nin,  ", 

(>)Min,  ", 

Min , 


Max  , 

& 

Max  , 

& 

Max 

X") 

Max 

4  X  "  ) 

3',1X,F8.3,5X)  ,314) 
order  M  of  XI 
oredr  N  of  X2  - 
IPRINT  "2A2) 
1  ...  to  calculate 


Min 


Min 


Min 


Max 
Max 


Min  , 
nin  , 


«a 


-  appr ox iMat i on 
appr ox iMat  i  on 


the 


n2A2) 
"2A2) 


2   "2A2) 
and  YO") 

•f-%  r-   o  ri 


absolute 
,4H"   ",9X"the  relative  error    "2A2) 
"  Change  coefficients    . . . 
"  Change  axes  .  .  , 

"  Enter  ALPHAX,  ALPHAY,  XO 
"< (3A2) ) ) 

"  1 " ////  "  T  »•■»  «.=>  «i  tie   <:l  ex  -v  *:x 
/16X,I2".  Mach  nuMber"/) 
6X"GaMMa"6X"Beta,,7X"Ph  i.  "/) 
6X"GaM*a"5X"Delta"7X"Phi"/) 
6X"GaMMa"6X"Beta"5X"X  vel "/) 
6X,,GaMMa"SX"Delta"5X"X  vel"/) 
IX, 4( 1X.F9.3) ) 
/9XHX1,,8X"X2"9X"Y") 
(3X,i6(6X,I2)> ) 
lX>I2,16(iX,F7.3)/3(3X,16(iX,F7.3>/)) 


error"/8X 


■r  a.  :i. 


3  A 
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0077  604  FORMAT  (.///"    Absolute  error  at  each  point    C  error  =  (calculated 

0078  *value  -  Measured  value)  )"/) 

0079  644  FORMAT  (///"  Relative  error  at  each  poxnt  in  %    C  error  ■  (calcu 

0080  *ated  value  -  Measured  val ue ) /Measured  value  3"/) 

0081  606  FORMAT  ( IX , 12  ,  18 ( IX , F6 . 1 ) /3( 3X , 18( IX ,F7 . 3) /) ) 

0082  1101  FORMAT  ( 8H   Beta   ) 

0083  1102  FORMAT  ( 8H   Delta  ) 

0084  1103  FORMAT  ( 8H   Phi     ) 

0085  1104  FORMAT  < 8H   X  vel  ) 

0086  LI  =  LOGLU(I) 

0087  LO  =  6 
0088 

0  0  89 
0090 

0091  C  

0092  C 

0093  C  .   read  raw  data  file. 

0094  C 

0095  C  

0096  41  WRITE  (LI,  121)  NOCR 

0097  READ   (LI,  122)  IFILE 

0098  WRITE  (LI,  149)  ICLR 

0099  42  WRITE  (LI,  123)  NOCR 

0100  READ   (LI,    *)  IPX 

0101  WRITE  (LI,  149)  ICLR 

0102  IF  (IPX.LT.l  .OR.  IPX.GT.2)  GO  TO  42 

0103  43  CONTINUE 

0104  IF  (IPX.EQ.l)  WRITE  (LI,  124)  NOCR 

0105  IF  (IPX.EQ.2)  WRITE  (LI,  125)  NOCR 

0106  READ   (LI,    *)  IBD 

0107  WRITE  (LI,  149)  < ICLR ,1=1,2,1 > 

0108  IF  (IBD.LT.i  .OR.  IBD.GT.2)  GO  TO  43 

0109  IF  (IPX.EQ.l)  IRECY=13 

0110  IF  (IPX.EQ.2)  IRECY=9 

0111  IF  (IPX.EQ.l  .AND.  IBD.EQ.l  )  WRITE  (LI,  131) 

0112  IF  (IPX.EQ.l  .AND.  IBD.EQ.2  )  WRITE  (LI,  132) 

0113  IF  (IPX.EQ.2  .AND.  IBD.EQ.l  )  WRITE  (LI,  133) 

0114  IF  (IPX.EQ.2  .AND.  IBD.EQ.2  )  WRITE  (LI,  134) 

0115  CALL  OPEN  ( IDCB , IERR , IFILE , IOPTN , 0 ,28 , 144 ) 

0116  IF  (IERR.LT.0)  GO  TO  41 

0117  CALL  READF  ( IDCB , IERR , XI .512 . LEN , 1 ) 

0118  IF  (  IERR   .LT.   0  )  STOP  0002 

0119  CALL  READF  ( IDCB , IERR ,X2 ,512 ,LEN ,5) 

0120  IF  (  IERR   .LT.   0  )  STOP  0003 

0121  CALL  READF  ( IDCB , IERR ,YT ,512 , LEN , IRECY > 

0122  IF  (  IERR   .LT.   0  )  STOP  0004 

0123  CALL  CLOSE  ( IDCB, IERR ,0) 

0124  44  WRITE  (LI,  126)  NOCR 

0125  READ   (LI,    *)  NMACH,NPITCH 

0126  WRITE  (LI,  149)  ICLR 

0127  IF  (NMACH.GT.i6)  GO  TO  44 

0128  IF  (NPITCH.GT. 16)  GO  TO  44 

0129  DO  801  I=l,NMACH.i 

0130  DO  801  J=l,NPITCH,i 

0131  801  Y(I,J)=YT(I,J) 
0132 

0133 
0134 

0135  C  

0136  C 

0137  C  .   Correct  input  data  and  optional  output  of  data,  coefficient 

0138  C  .   are  based  upon. 

0139  C 

0140  C  

0141  IF  (IBD.EQ.2)  GO  TO  48 

0142  DO  47  I=1,NMACH,1 

0143  DO  47  J=1,NPITCB,1 

0144  47  X2(I,J)=X2(I,J)/Xi(I,J) 

0145  48  WRITE  (LI,  160)  NOCR 

0146  READ   (LI,    *>  IDUM 

0147  WRITE  (LI,  149)  ICLR 

0148  IF  (IDUM.NE.i)  GO  TO  55 

0149  WRITE  (LO,  611)  IFILE 

0150  DO  54  I=1,NMACH,1 

0151  WRITE  (LO,  612)  I 

0152  IF  (IPX.EQ.l  .AND.  IBD.EQ.l)  WRITE  (LO,  613) 
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0153  IF  (IPX.EQ.i  .AND.  IBD.EQ.2)  WRITE  ( LO ,  614) 

0154  IF  (IPX.EQ.2  .AND.  IBD.EQ.I)  WRITE  ( LO ,  615) 

0155  IF  (IPX.EQ.2  .AND.  IBD.EQ.2)  WRITE  ( LO ,  616) 

0156  DO  54  J=l,NPITCH,i 

0157  54  WRITE  (LO,  617)  XI ( I , J ) ,X2( I , J ) , Y ( I , J ) 

0158  WRITE  (LO,  618) 
0159 
0160 
0161 

0162  C 

0163  C 

0164  C      .  .Calculate  Mininun  and  noxinun  values  for  all  coordinates  and 

0165  C         optional  redefinition  of  nininun  and  naxinun  values  by  user. 

0166  C      . 

0167  C      

0168  55  XUMAX=Xi(l ,1) 

0169  XUMIN=Xi( 1, 1) 

0170  YUMAX=X2(l,i) 

0171  YUMIN=X2(i,i) 

0172  ZUMAX=Y(i,i) 

0173  ZUMIN=Y(l,i) 

0174  DO  45  I=l,NMACH.l 

0175  DO  45  J=i.NPITCH,i 

0176  IF  (Xi(I,J) .GT.XUMAX)  XUMAX=X1 ( I , J ) 

0177  IF  (Xi(I,J> .LT.XUMIN)  XUMIN=Xi ( I , J ) 

0178  IF  (X2(I,J) .GT.YUMAX)  YUMAX=X2( I , J ) 

0179  IF  (X2(I,J) .LT.YUMIN)  YUMIN=X2( I , J ) 

0180  IF  (Yd, J)  .GT.ZUHAX)  ZUMAX=Y(I,J) 

0181  IF  (Yd, J)  .LT. ZUMIN)  ZUMIN=Y(I,J) 

0182  45  CONTINUE 

0183  IF  (IPX.EQ.i  .AND.  IBD.EQ.i)  WRITE  (LI,  127) 

0184  IF  (IPX.EQ.I  .AND.  IBD.EQ.2)  WRITE  (LI,  128) 

0185  IF  (IPX.EQ.2  .AND.  IBD.EQ.I)  WRITE  (LI,  129) 

0186  IF  (IPX.EQ.2  .AND.  IBD.EQ.2)  WRITE  (LI,  130) 

0187  WRITE  (LI,  150)  XUMIN , YUMIN , ZUMIN 

0188  WRITE  (LI,  151)  XUMAX , YUMAX ,ZUMAX 
0189 
0190 
0191 
0192 

0193  WRITE  (LI,  153)  NOCR 

0194  READ   (LI,    *)  I 

0195  WRITE  (LI,  149)  1CLR 

0196  IF  (I.NE. 1)  GO  TO  49 

0197  IF  (IPX.EQ.i  .AND.  IBD.EQ.I)  WRITE  (LI,  154) 

0198  IF  (IPX.EQ.I  .AND.  IBD.EQ.2)  WRITE  (LI,  155) 

0199  IF  (IPX.EQ.2  .AND.  IBD.EQ.i)  WRITE  (LI,  156) 

0200  IF  (IPX.EQ.2  .AND.  IBD.EQ.2)  WRITE  (LI,  157) 

0201  READ   (LI,    *)  XUMIN  , XUMAX , YUMIN , YUMAX , ZUMIN ,ZUMAX , INCX , INCY , INCZ 

0202  WRITE  (LI,  158)  XUMIN, XUMAX, YUMIN, YUMAX, ZUMIN, ZUMAX , INCX , INCY , INCZ 
0203 

0204 
0205 

0206  C 

0207  C 

0208  C      .   Initialize  plotter;   define  plotter  area;   calculate  scaling 

0209  C        coefficients  XA,  XB,  YA,  YB,  ZA  and  ZB. 

0210  C 

0211  C 

0212  49  CALL  INITG  (13) 

0213  XPMIN   =  0. 

0214  XPMAX   =  10. 

0215  YPMIN   =  0. 

0216  YPMAX   =  10. 

0217  ZPMIN   =0. 

0218  ZPMAX   =  10. 

0219  ALPHAX  ■  330. 

0220  ALPHAY  =   30. 

0221  ALPHAZ  =90. 

0222  XA      =  (XPMAX       -XPMI.N       ) /(XUMAX-XUMIN) 

0223  XB      =  (XPMIN*XUMAX-XPMAX*XUMIN)/(XUMAX-XUMIN) 

0224  XL      =  (XPMAX       -XPMIN       ) 

0225  YA      =  (YPMAX       -YPMIN       ) /(YUMAX- YUMIN) 

0226  YB      =  (YPMIN*YUMAX-YPMAX*YUMIN)/(YUMAX-YUMIN) 

0227  YL      =  (YPMAX       -YPMIN       ) 

0228  ZA      =  (ZPMAX       -ZPMIN       )/ (ZUMAX-ZUMIN ) 
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WRITE 

(LI, 

150) 

WRITE 

(LI 

151) 

WRITE 

(LI 

152) 

INCX  = 

=  4 

INCY  = 

«  4 

INCZ  =  4 

WRITE 

(LI, 

153) 

READ 

(LI, 

*) 

WRITE 

(LI, 

149) 

IF  (I, 

ne.  i: 

•  GO 
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0229 
0230 
0231 
0232 
0233 
0234 
023S 
0236 
023? 
0238 
0239 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
0247 
0248 
0249 
0250 
0251 
0252 
0253 
0254 
0255 
0256 
0257 
0258 
0259 
0260 
0261 
0262 
0263 
0264 
0265 
0266 
0267 
0268 
0269 
0270 
0271 
0272 
0273 
0274 
0275 
0276 
0277 
0278 
0279 
0280 
0281 
0282 
0283 
0284 
0285 
0286 
0287 
0288 
0289 
0290 
0291 
0292 
0293 
0294 
0295 
0296 
0297 
0298 
0299 
J300 
0301 
0302 
0303 
0304 


=  (ZPMIN*ZUMAX-ZPMAX#ZUMIN)/(ZUMAX- 
=  (ZPMAX       -ZPMIN       ) 
=  8, 


85 


86 


87 

88 
301 


302 
303 


01 


02 


91 


92 


(ALPHAX*PI)/180 
(ALPHAY*PI)/18Q 
(ALPHAZ*PI)/180 
XO 
YO 
XA 
XB 
YA 
YB 
ZA 
ZB 
IBD 


ZB      =  (ZPMIN*ZUMAX-ZPMAX#ZUMIN)/(ZUMAX-ZUMIN) 

ZL 

XO 

YO      =  12 

PL0TR(61)   = 

PL0TR(62)   = 

PL0TR(63)   = 

PL0TR(64)   = 

PL0TR(65)   = 

PLUTR<66)   = 

PL0TR(67)   = 

PL0TR(68)   = 

PL0TR(69)   = 

PLOTR(70)   = 

PL0TR(71)   = 

GO  TO  (86,87) 

CALL  CODE 

WRITE  (YTITLE,1101) 

GO  TO  88 

CALL  CODE 

WRITE  (YTITLE.1102) 

GO  TO  (301,302)  IPX 

CALL  CODE 

WRITE  (ZTITLE,1103) 

GO  TO  303 

CALL  CODE 

WRITE  (ZTITLE,tl04) 

CALL  SETSM  (113,1.  ) 

CALL  AXIS  (X0,Y0,XL,ALPHAX,6H  Gantna,  8 . XUMIN , XUMAX , 4HF6 . 2 ,6 > INCX) 

CALL  AXIS  (X0,Y0,YL,ALPHAY,YTITLE,-8,Y0MIN,Y0MAX,4RF6.2,6,INCY) 

CALL  AXIS  (X0,Y0,ZL,ALPHAZ,ZTITLE,-8,ZUMIN,ZUMAX,4HF6.2,6,INCZ) 


Plot  Measured  calibration  surface. 


CALL  SETSM  (113,2. ) 

DO  01  I=1,NMACH,1 

DO  01  J=i,NPITCH,i 

CALL  THRTW  ( XPLOT , YPLOT >X1 ( I , J ) . X2( I , J ) , Y( I  ,  J ) ) 

IF  (  J   .EQ.   .1  )  CALL  PLOT  ( XPLOT ,  YPLOT, 2 ) 

IF  (  J   .GT.   1  )  CALL  PLOT  ( XPLOT , YPLOT , 3) 

CONTINUE 

DO  02  J=1,NPITCH,1 

DO  02  I=1,NMACH. 1 

CALL  THRTW  (XPLOT , YPLOT . XI ( I . J ) ,X2( I , J )  ,  Y ( I ,  J) ) 

IF  (  I   .EQ.   1  )  CALL  PLOT  t XPLOT , YPLOT ,2) 

IF  (  I   .GT.   1  )  CALL  PLOT  ( XPLOT , YPLOT ,3 ) 

CONTINUE 


Calculate  calibration  surface  coefficients 


CALL  SETSM  (113,0. ) 

WRITE  (LI,  101)  NOCR 

READ   (LI,    *)  MORDER 

WRITE  (LI,  149)  ICLR 

IF  (MORDER. GT. 6)  GO  TO  91 

WRITE  (LI,  102)  NOCR 

READ   (LI,    *)  NORDER 

WRITE  (LI,  149)  ICLR 

WRITE  (LI,  107)  NOCR 

READ   (LI,    *)  IPRINT 

WRITE  (LI,  149)  ICLR 

IF  (NORDER .GT. 6)  GO  TO  92 

M=MORDER+l 

N=NORDER+l 

CALL  MAT3  ( NMACH ,NPITCH ,M ,N ,COEFF , IPR INT) 
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0305 
0306 

0307  C  

0308  C 

0309  C  .   Overwrite  data  array  with  calculated  data,  based  on  the 

0310  C  .   just  obtained  coefficients. 

0311  C  . 

0312  C  

0313  411  WRITE  (LI,  105)  NOCR 

0314  READ   (LI,    *)  IAR 

0315  WRITE  (LI,  149)  ( ICLR  ,  1  =  1 ,2 ) 

0316  IF  (IAR.LT.i  .OR.  IAR.GT.2)  GO  TO  411 

0317  IF  (IAR.EQ.l)  WRITE  (LO,  604) 

0318  IF  (IAR.EQ.2)  WRITE  (LO,  644) 

0319  WRITE  (LO,  602)  ( J , J=l , NPITCH , 1 ) 

0320  DO  20  I=l,NMACH,i 

0321  DO  19  J=l,NPITCH,i 

0322  SU«=0. 

0323  DO  96  11=1 ,M,1 

0324  IEXPi=Il-i 

0325  IF  (IEXP1.EQ.0)  GO  TO  401 

0326  XiEXP=Xi(I,J)**IEXPl 

0327  GO  TO  402 

0328  401  X1EXP=1. 

0329  402  CONTINUE 

0330  DO  96  Ji=l,N,i 

0331  IEXP2=Ji-l 

0332  IF  (IEXP2.EQ.0)  GO  TO  403 

0333  X2EXP=X2(I,J)**IEXP2 

0334  GO  TO  404 

0335  403  X2EXP=i. 

0336  404  CONTINUE  » 

0337  96  SUH=SUM+C0EFF(I1,J1)*X1EXP*X2EXP 

0338  GO  TO  (412,413)  IAR 

0339  412  R<J)»SUM-YCI.J) 

0340  GO  TO  19 

0341  413  R(J)=((SUM-Y(I,J>)/Y(I,J))*100. 

0342  19  Y(I,J)=SUtt 

0343  20  WRITE  (6,603)  I , ( R ( J ) , J=l , NPITCH , 1 ) 

0344  WRITE  (6,602)  ( J,  J=l , NPITCH, 1 ) 
0345 

0346 
0347 

0348  C  

0349  C 

0350  C  .   Plot  calculated  calibration  surface. 

0351  C 

0352  C 

0353  CALL  SETSM  (113,3.) 

0354  DO  81  I=i,NMACH,l 

0355  DO  81  J=i,NPITCH,i 

0356  CALL  THRTW  ( XPLGT , YPLOT ,Xi ( I , J )  , X2( I , J ) , Y ( I ,  J ) ) 

0357  IF  (J.EQ.l)  CALL  PLOT  (XPLOT ,  YPC.OT,  2) 

0358  IF  (J.GT.l)  CALL  PLOT  (XPLOT , YPLOT ,3 ) 

0359  81  CONTINUE 

0360  DO  82  J=i, NPITCH. 1 

0361  DO  82  I  =  i,NHACH.  1 

0362  CALL  THRTW  (XPLOT , YPLOT ,X1 ( I , J ) . X2( I , J ) , Y ( I  ,  J ) ) 

0363  IF  (I.EQ.i)  CALL  PLOT  ( XPLOT, YPLOT , 2> 

0364  IF  (I.GT.l)  CALL  PLOT  ( XPLOT , YPLOT ,3 ) 

0365  82  CONTINUE 
0366 

0367 
0368 

0369  C 

0370  C 

0371  C  .   Redefine  Y. 

0372  C  .                                                                   .  : 

0373  C  ^ 

0374  DO  802  I=l,NMACH,i 

0375  DO  802  J=i, NPITCH 

0376  802  Y(I>J)=YT(I, J) 
0377 

0378 
0  379 

0380  C      
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0381  C 

0382  C      .   Next  step  ? 

0383  C 

0384  C      

0385  URITE  (LI, 103)  NOCR 

0386  READ   (LI,   *>  IDUM 

038V         URITE  (LI,  149)  ( ICLR  .  1  =  1 , 2 , 1  ) 

0388  IF  (IDUM.EQ.i)  GO  TO  91 

0389  IF  (IDUM.EQ.2)  GO  TO  89 

0390  CALL  STOPG 

0391  STOP  077 

0392  89  URITE  (LI, 104) 

0393  READ   (LI,   * )  ALPHAX , ALPHAY ,XO , YO 

0394  URITE  (LI,  149)  ICLR 

0395  GO  TO  85 

0396  END 


FTN4  COMPILER:  HP92060-16092  REU.  1926  (790430) 

**   NO  UARNINGS  **   NO  ERRORS  **    PROGRAM  -  04465       COMMON  =  00000 
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3.6.2.   Load  map 


LAB     10043  20622   Plot  calibration  points  and  cnlrulnt*  c oof f ic len ts 
to  load  these  program  modules,  enter  (from  LOADR):   MS , STPLBL 


DATA3 

MAT3 

MAT31 

S3 

AB3 


20623  23622 
23623  27455 
27456  27651 
27652  30070 
30071  41534 


3D-Appr oximat ion'  /  DATA3  / 

3D-Appro  ximation  arrange  sysren  matrix  and  i/*ctor 

3D-Approx  mat  ion        arrange  subnatrix  1 


3D-Appr o  xinat ion 
3D-Aqdpox  a wat ion 


/    AFt3 


/ 


c  lass 
compute    supinations 


AFLD 

AXIS 

INITG 

PLOT 

SETSM 

STOPG 

THRTU 

FCTR 

FACTR 

LOGLU 
READF 
OPEN 
CLOSE 
CLRIO 
CURD. 
♦SMVE 
LURQ 
.DADS 
.DMP 
.  DDI 
SESSN 
R/U* 
P.  PAS 
.  DNG 
PAUSE 
$ALRN 
FMTIO 
ERR0 
TAN 
ABS 
.SNCS 
.DDE 
.DIN 
.RTOI 
.FPWR 
.SBT 
.  .FCM 
PAU.E 
ER0.E 
.CURS 
$OPEN 
RWtUB 
RUND* 
FRMTR 
FMT.E 
IB4A2 
REIO 
RMPAR 
PNAME 
LUTRU 
IPUT 
*SETP 
.CFER 
.LBT 


41535 
42535 
45055 
46070 
46621 
47562 
47625 
50112 
50116 

50136 
50214 
51156 
51454 
51664 
51673 
S1674 
51763 
52346 
52456 
52624 
53125 
53143 
53302 
53331 
53341 
53442 
53560 
55057 
55147 
55254 
55263 
55425 
55437 
55445 
55541 
55603 
55644 
55661 
55662 
55663 
55747 
56124 
56472 
56615 
62253 
62254 
63555 
63702 
63745 
64013 
64122 
64143 
64170 
64246 


42534 
45054 
46067 
46620 
47561 
47624 
SOiii 
50115 
50135 

50213 
51155 
51453 
51663 
51672 
S1673 
51762 
S2345 
52455 
52623 
53124 
53142 
53301 
53330 
53340 
53441 
53557 
55056 
55146 
55253 
55262 
55424 
55436 
55444 
55540 
55602 
55643 
55660 
55661 
55662 
55746 
56123 
56471 
56614 
62252 
62253 
63554 
63701 
63744 
64012 
64121 
64142 
64167 
64245 
64276 


/  AFLD   /  A<256) 

read  control  array 


GSP 

GSP 

GSP 

GSP 

GSP 

GSP 

GSP  converts  3D  user  coordinates  to  plotter  units 

GSP   /  FCTR   /  FX,FY 

GSP  vary  size  of  the  plot 


Draw  and  label  axes, 
A;   initialize  plotter 
nove  pen  to  a  defined  point 
change  plotter  nodes 
terminate  grafics 


92067-16268  REV. 1903 

92067-16125  REV. 1940 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

750701   24998-16001 

92067-16125  REV. 1903 

92067-16268  REV. 1903 

92067-16268  REV. 1903 

780818   24998-16001 

780818   24998-16001 

781021   24998-16001 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

92067-16125  REV. 1903 

780818   24998-16001 

771122   24998-16001 

92067-16268  REV. 1903 

24998-16002  REV . 192 

771122   24998-16001 

780424   24998-16001 

750701   24998-1600 

780424   24998-16001 

780818   24998-16001 

780818   24998-16001 

780921   24998-16001 

781106   24998-16001 

770518   24998-16001 

750701   24998-16001 

750701   24998-16001 

750701   24998-16001 

780424  24998-16001 
92067-16125  REV. 1903 
92067-16125  REV. 1903 
92067-16125  REV. 1903 
24998-16002  REV. 192 
24998-16002  REV. 190 
59310-1X013  REV  1940 
92067-16268  REV. 1903 
781106   24998-16001 

771121   24998-16001 

92067-16268  REV. 1903 

92067-16125  REV. 1903 

781106   24998-16001 

750701   24998-16001 

770518   24998-16001 


790228 
790719 
790215 
781229 

780526 
790202 
790223 


780413 
781214 
740801 


770715 
6  790417 


790103 

781003 

780801 
6  790503 
1  781107 

790802  1153 

790316 


790223 
740801 


24  PAGES  RELOCATED 


24  PAGES  REQ'D 


LINKS:BP     PROGRArhLB    LOAD:TE 
/LOADR: LAB     READY  AT   3:10  PM 

/LOADR : SEND 


NO  PAGES  FMA 


COMMONjNC 
FR1  .  ,  26   SEPT, 


NO  PAGES  hSEG 


1980 
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3.6.3.   Results 

i)   printed  output 
T  l-i  kk*.  «.r»  i»   c:l  ex  t  ex       -F  r»  o  r«i   -P  .i.  :i.  &       C  A  I.-  X  A  A 

1 .  Mach  nunber 
Ganna       Beta        Phi 


.810 

- 

.005 

-35, 

000 

.519 

- 

.005 

-25 

,  000 

.280 

- 

.005 

-15, 

000 

.182 

— 

.005 

-10 

,000 

.103 

- 

.006 

-5, 

000 

.059 

— 

.006 

-2, 

,500 

.018 

— 

.006 

0 

000 

.019 

— 

.  006 

2 

,500 

.054 

— 

.006 

5. 

000 

.133 

— 

,  006 

10, 

000 

.197 

— 

.006 

15, 

000 

.360 

- 

.  006 

2S, 

,  000 

.605 

- 

.006 

35, 

000 

2.  Mach  nunber 
Ganna      Beta       Phi 


-.781 

-.  016 

-3S. 000 

-.506 

-.016 

-25.000 

-.280 

-.  016 

-15.000 

-.178 

-.017 

-10.000 

-.092 

-.018 

-5.000 

-.048 

-.018 

-2.500 

-.004 

-.018 

0.  000 

.032 

-.018 

2.500 

.061 

-.018 

5.  000 

.130 

-.019 

10.000 

.190 

-.019 

IS. 000 

.346 

-.020 

25.000 

.583 

-.018 

35. 000 

3.  Mach 

i  nunber 

Ganna 

Beta 

Phi 

-.787 

-.024 

-35.0  00 

-.510 

-.024 

-2S. 0  00 

-.276 

-.024 

-15.000 

-.  172 

-.  026 

-10. 000 

-.085 

-.027 

-5.000 

-.050 

-.  027 

-2.500 

-.003 

-.027 

0.000 

.  037 

-.027 

2.500 

.066 

-.027 

5.000 

.  131 

-.027 

10. 000 

.191 

-.029 

15.000 

.34S 

-.029 

2S.000 

.592 

-.027 

35.000 

4.  Mach  nunber 
Ganna       Beta        Phi 


790 

- 

.035 

-35 

,  000 

519 

- 

.035 

-25, 

000 

276 

- 

.  035 

-IS, 

,  000 

174 

— 

.037 

-10, 

000 

055 

- 

.  039 

-S, 

000 

043 

— 

,039 

-2, 

500 

003 

- 

.  038 

0  , 

,  000 

037 

- 

.038 

2, 

500 

064 

— 

.  039 

S, 

000 

133 

— 

.039 

10, 

000 

196 

- 

.  041 

15, 

,  000 

349 

- 

.042 

25, 

000 

596 

- 

.  039 

35, 

,  000 

5.  Mach  nunber 
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GOMMQ 


Beta 


Phi 


79S 

-.045 

-35 

000 

521 

-.  045 

-25 

000 

282 

-.045 

-15 

000 

175 

-.  048 

-10 

000 

082 

-.050 

-5 

000 

040 

-.  050 

_2 

500 

010 

-.050 

0 

000 

034 

-.  049 

2 

500 

067 

-.049 

5 

000 

134 

-.  050 

10 

000 

198 

- .  053 

15 

000 

350 

-.  055 

25 

000 

593 

-.050 

35 

000 

6.  Mach  nunber 


Ganna 


Beta 


Phi 


803 

- 

.  06S 

-35 

000 

528 

- 

.065 

-25 

000 

284 

- 

.  065 

-IS 

000 

173 

- 

068 

-10 

000 

085 

- 

.  071 

-5 

000 

032 

- 

072 

-2 

500 

004 

- 

.  071 

0 

000 

024 

- 

.071 

-> 

500 

058 

- 

.  071 

S 

000 

131 

- 

.071 

10 

000 

199 

- 

.  075 

15 

000 

349 

- 

.078 

25 

000 

603 

- 

.  073 

35 

000 

XI 


X2 
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ii)   graphic  output 
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FOUR  DIMENSIONAL  APPROXIMATION 


4.1.   Problem: 


A  data  set  of  NPNTS1  *  NPNTS2  *  NPNTS3  data  points  is  given, 
where  Y  depends  on  the  parameters  XI,  X2  and  X3.   The  data  pattern  is 
to  be  approximated  by  a  function  Y  =  f(Xl,X2,X3),  so  that  the  error 
between  data  points  and  analytically  determined  points  is  lower. 
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4.2.   Approach 

Again,  as  in  the  two  and  three  dimensional  case,  a  polynominal 
method  is  used.   For 

we  assume 

V-   fa,  *  c*,X3*  •■■   J  C„0X3     J* 


or 


XX 
y? 


*f«w  <  <W**  ■-    *  ^  ** '  "J '  «  "  V  * 


*/^,.  *C.aAT*...    "C-^W     7< 


J  fl=L,  -  C«W  ....  *C^.  *?  "1-XZ  "]*l  * 


N-l 


"  c cui  +  clu*3 ■>  -  *  Cu.»XS""7-*2 


<[cLn  +CLHlx+  ...  ^CLf4f/x/"j-xZ^/x/L 


(4.1.) 


50 


This  data  arrangement  which  discriminates  between  data  points  taken  at 
constant  XI,  X2  and  X3,  respectively  is  needed  in  the  application  of 
this  method  to  the  description  of  pneumatic-velocity  probe  characteristics 
The  various  surfaces  in  the  above  sketch  for  example  can  be  seen  as  data 
taken  at  different  Mach  number  (XI).   At  each  constant  Mach  number 
(XI  =  constant)  yaw  angle  (X2)  and  pitch  angle  (X3)  are  varied,  and  the 
probe  is  measured  to  be  Y.   Even  though  probes  are  usually  balanced  to 
the  average  flow  yaw  angle  in  some  cases  this  is  not  possible  and  the 
four  dimensional  approximation  is  then  needed.   A  similar  method  for 
describing  probe  calibration  data  is  used  at  the  institute  for  Jet 
Engines  and  Turbomachines  at  the  Technical  University  Aachen  (Ref.  2). 
The  least  squares  criterion  is  then 

where  the  indices  rst  denote  data  points.   Using  equation  (4.1.),  R 
becomes 

XfMBi      MPffltt      MPNfiJ  m-| 

r.  i        **i       t* ' 

rt-i 


4  (c    4c    yj    -*      *  C      Yl    t  )'*2  i  + 


#    « 


tU        ^M        rsi       •••  UN 
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The   term   in   the    [      ]    -   bracket    shall   be   called   B 

J**TS/     f+HKh    NPNX3 

r»/      e»i      t-/ 
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4.3.   Solution: 


To  minimize  the  error,  R  is  partially  differentiated  with  respect 
to  the  coefficients  C...  and  the  partial  derivatives  are  set  to  zero, 
thus: 


5£ 

dCgk 


=  0 


/'-  /...  L   •  /--  / 


H    ;  J,*     /      ...        /V 


^  P  ttfinm      NtHFSZ      N*SIKS  — v    p 

O^^k        r*l         -"/    t*l  (-X-iiK 


Assuming,  that  the  summations  extend  over  all  data  points,  ZZI  should 

UPHTTI        *****    M»WTJ 

be  understood  to  be    /_i       z_j   /.    •   Now  even  a  very  large  sheet 

r*/   ft      t*t 

of  paper  can't  take  the  system  matrix,  unless  we  skip  to  submatrix 
notation. 


A 


Vi 


x 


K.       ft, 


a 


K    K, 


A) 


li 


K 


< 

B,* 

c* 

= 

S* 

3     t>    £> 

< 

(4.3.) 


The  asterisk  indicates  a  submatrix  1.  class.   On  the  following  pages 
the  submatrices  and  subvectors  1.  class  are  listed. 
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4.4.   Structure  of  Equation  System: 
4.4.1.   Svstem  Matrix  A 


Ir.  can  be  seen,  that  all  submatrices  1.  class  on  the  same 
diagonal  bands  2.  order  of  the  system  matrix  are  identical  and  therefore 


will  be  renamed  from  A 


i,j 


* 


At  -  A 


I4 
9 


A 


13 


2*  * 


A 


A 


t*i 


-A* 


■/A  * 


-  A 

-  Ah 


K 


-A* 


The  submatrices  1.  class  A,  (k*l, . . . ,2L-1)  themselves  are  divided 
up  into  submatrices  2.  class. 


K- 


-A 


** 


?A 


A 


A 


A** 


/A** 
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The  double  asterisk  indicates  a  submatrix  2.  class.   All  submatrices 

2.  class  on  the  same  diagonal  bands  2.  order  of  the  submatrices  1.  class 

A,   (k  =  1,...,2L)  are  identical  and  therefore  will  be  renamed  from 
**         ** 

'\;i,j    °  \;l 


it 


** 


K-*k 


/A 


<■* 


*+ 


-  A 

^,2,1  ^V/,-2 


/A  **    -  JA  *  *      -A*' 


--  A 


*4 


A 


44 


K 


A 


/A*;«.Z 


. -*« 


\kf/4  .J 


A**  A** 


4  A 


A**  A  A  4 


** 


&K-2M-I 


-A 


*4 


where 


** 


can  be  written  as 


A 


2 


EEDC«    EEDC<«, 


EEEx£;OC 


EM4X 


EI^OOC 


k-ljmtZL-l;t'i....,2P-l 


63 


Note,  that  this  submatrix  2.  class  is  not  only  symmetrical,  but  also  the 

elements  on  each  diagonal  line  2.  order  are  identical.   The  elements 

\,  ,  .  of  the  submatrix  2.  class  A,    therefore  are  renamed  to  a,  ,   ; 
;i;i,j  Tc;l  k;L;m 


** 


•ix 


L  specifies  the  place  of  the  submatrix  2.  class  A^#   in  the  submatri: 

KJ 1 

*  * 

1.    class  A,  ,   k  specifies   the  place  of   the   submatrix  1.    class  A,     in   the 

system  matrix  A  and   i,j    specifies   the  place  of   the  element  a,         .    .    in 


the  submatrix  2.    class 


Qk,l;l     '       °-lr,l,'./ 


\l;Z         %tf2,i  '    aM;'J 


%t}2     '   ak.L;Zi    '  al.;i,z.2-~  %,'J 


-  a 


-a, 


k,l,H  k;l;N,l  ^k;l;N-l,Z  k;l}tH,2 


=  %t,2.*-i  ''      ht?i,» 


%h** 


akjljH,Z    '%l;W 


\i;2/<~<   "  a*><.j2-.« 


%L;U*l 


=  a 


k;L,  «.! 


s%i;+,>"  'aW 


%L}2tH 


where 


%i, 


can  be  written  as 


-a 


M/  Mi* 


kj  t;m*      rtf  1 1\   L  &     "**   "^ 


(4.4.) 


k.l,...jL-i  j  C-'y  ....JM-i  j    m-Ij^./N-I 
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This  operation  is  programmed  in  REAL  FUNCTION  S4  (NPNTS1,NPNTS2,NPNTS3, 

IPOWRl,IPOWR2,IPOWR3,IY) .   The  data  XI,  X2,  X3  and  Y  are  known  to 

this  function  through  a  common  block  named  DATA4.   So  only  NPNTS1, 

NPNTS2,  NPNTS3,  IP0WR1  (=»  k-1),  IP0WR2  (=  L-l),  IP0WR3  O  m-1)  and 

IY  (=0)  have  to  be  passed  to  the  function  and  the  value  of  a,  ,    is 

tc;l;m 

returned  through  the  function  name  S4. 


** 


SUBROUTINE  MAT42  presets  the  submatrices  2.  class  A     in  the 

*  >  1 

following  way 

i)   Preset  edge  section  elements  (using  REAL  FUNCTION  S4) 


%/v/v 


F><Z.  ,  ,0. 


ak;t,  N+t 

•    •    • 


ak;L,l  ~*   ak;l;/y 
ftp  «(b/   S*cfrs*~ 

ak-.l,r***    **    aX;l;2N-f 


ii)      Copy  defined  elements  diagonally 


a«M>*  aw*s?  aM>{  -  w     -  ak,L-  '•» 


J^    k>P?  P^  ^> 


•  •  • 
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G*;/;/,/         *;l:i,2      ^*,(,',r 


aW;  l,N 


■k;l;  V,/V 
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This  subroutine  that  needs  NPNTS1,NPNTS2,NPNTS3 , IPOWRl(=k-l)  and 
IP0WR2(=L-1)  as  input  parameters,  returns  the  submatrix  2.  class 

A,    (SUBM2(4,4)J  to  SUBROUTINE  MAT41,  which  foresets  the  submatrices 

* 

1.  class  A,  in  the  following  way 

K 

i)   Preset  edge  section  submatrices  (using  SUBROUTINE  MAT42) 


C**  v*  *V/*- 


4* 

AA* 

•  •  * 

A 


A* 

k;2H-l 


ii)      Copy  defined  submatrices   diagonally 
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k,i.ff 

4* 

** 


A** 


A 


*;/.' 


A 


A 


A 


Km 
** 


** 


\l,S 
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A*4, 


«-K 


%2,t      Ink;22 


/ 


A 


A 


kl.M 


A 


■f-K 


A 


** 


K 
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44 


^lr;U 
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hJ.r 


'   k  ,n,n 
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S  S  iS  iS  , 


Sk 
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4-t 
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This  subroutine,  that  needs  NPNTS1,NPNTS2,NPNTS3,M,N,IP0WR1(=L-1) 

* 
as  input  parameters,  returns  the  subroutine  1st  class  A^  (SUBM1(16,16))  to 

SUBROUTINE  MAT4,  which  presets  the  entire  system  matrix  A  in  the 

following  way. 
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i)   Preset  edge  section  submatrices  (using  SUBROUTINE  MAT41) 


A  *=$  h[  =£  A*  =*  ...  =£  AL 

Ah 


^2 


O 

£ 
**&-' 


4  -  ^ 


Z". 


2L-t 


I^Ofitj    /*o~</  &(ar    s+d-) 


ii)      Copy  defined  submatrices   diagonally 
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<       A* 
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4.4.2.   Right  Hand  Side  Vector  B; 

* 
The  right  hand  side  subvectors  1.  class  B  generally  can  be 


written  as 


a 


E  2  EXi*&*3* 


E  £  JV^'A 

E  S  SY^xint-x2^ 
E  2  SY«(XiBt'X2^-X?Kt, 


w 


.n-i 


EEEYniXi;;'X2Kl-X!ot 


E  S  2Xi-«I 


SEEXfXi^XT^ 


k-»  ^^  v;""1 


EZEfcxCoa£x& 


k*'. 
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The  subvectors  1.  class  B,  (k 

k 

up  into  subvectors  2.  class: 


=  1,...,L)  themselves  will  be  divided 


*:- 


e. 


*  * 


k.L 


** 


where  subvectors  2.  class  B,  .  can  be  written  as 

k;L 


'kd 


Finally,  the  elements  of  the  subvectors  2.  class  B    can  be  written  as 


A..    -  SSSX-.xtf-xz^x?^ 


*;</< 


**/. 


wl  "  wfc 


,Z;  *•',...,";  m 


--  I 


N 


To  calculate  b,  ,   ,  REAL  FUNCTION  S4  is  used,  again.   Data  XI,  X2,  X3 
lc;l;m 

and  Y  are  available  through  COMMON  block  DATA4.   NPNTS1,  NPNTS2,  NPNTS3, 
IP0WR1  (-  k-1),  IP0WR2  (=»  L-l)  IP0WR3  (»  m-1)  and  IY  (=1)  have  to  be 
passed  to  the  function  and  the  value  of  b  .    is  returned  through  the 
function  name  S4.   Using  FORTRAN  programming  language,  the  allocation 

of  b,  ,    can  be  oerformed  by  three  stacked  DO-loops. 

<;  1  ;m 
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1=0 

DO   07   L1-1,L,1 

IPOWR1-L1-1 

DO   08  M1-1,M,1 

ipowr2=mi-i 

DO  08  N1=1,N,1 
IPOWR3-N1-1 

1=1+1 
08   B(I)-S4 (NPNTS1 ,NPNTS2 ,NPNTS3 , IPOWR1 , IPOWR2 , IPOWR3 ,1) 
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4.5.   Software: 

The  software  to  compute  the  coefficients  for  a  two  dimensional 

approximation  is  described  in  APPENDIX  A  and  is  implemented  in  the  TPL 

HP»21MX  computer  system.   To  work  correctly  with  these  program  modules, 

the  user  has  to  conform  to  the  following  conventions. 

i)  Provide  the  data  in  four  arrays 
(Type:  REAL)  of  5*5*5*  elements 
through  a  COMMON  block,  named  DATA4. 

COMMON  /  DATA4  /  X1,X2,X3,Y 

REAL  X1(5,5,5),X2(5,5,5),X3(5,5,5),Y(5,5,5) 

ii)  Dimension  an  array  (Type:  REAL)  of 

5*5*5  elements  to  contain  the  coefficients, 

REAL  COEF  (5,5,5) 

iii)   Define  the  parameters  NPNTS1,  NPNTS2, 

NPNTS3,L,M,N  and  IPRINT 

NPNTS1...   #  of  XI  variations 
L  S  NPNTS1  H   5 

NPNTS2...  If   of  X2  variations 
M  S  NPNTS2  i   5 


NPNTS3. 


M 


#  of  X3  variations 
N  S  NPNTS3  £  5 

(desired  order  of  approximation 
polynominal  w.r.t.  XI)  +  1 

ULi4 

(desired  order  of  approximation 
polynominal  w.r.t.  X2)  +  1 
1  i,   M  i   4 

(desired  order  of  approximation 
polynominal  w.r.t.  X3)  +  1 
i  S  N  U  4 
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IPRINT  .  .  .   Controls  quantity  of  print  out 

2  ...   Print  system  matrix  and  right 

hand  side  vector  before  and 

after  Gauss  Jordan  Elimination. 
1  ...   Print  equation  system  after 

Gauss  Jordan  Elimination 
<0. . .   No  print  out 
>0...   Print  equation  (4.1.)  with 

the  actual  parameters. 

iv)   If  a  user  program  uses  subroutine  MAT4,  the  software 
modules  have  to  be  loaded  using  the  procedure  in 
section  1.4. iv. 

If  all  these  requirements  are  met,  the  correct  call  for  the  subroutine  is: 
CALL  MAT4  (NPNTS1,NPNTS2,NPNTS3,L,M,N,C0EF, IPRINT) 

Upon  completed  execution  of  this  approximation  routine  array  COEF 
contains  the  coefficients.   Externals  used  by  MAT4:   AB4,  DATA4,  MAT41, 
MAT42,  S4.   Under  no  circumstances  may  the  user  use  any  of  these  names  for 
a  modules  of  his  own  user  program. 
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^•6.  Sample  User  Program 

M.6.1        FTNU   Compiler   listing  of  sample   user  program 

0001  FTN4,L 

0002  PROGRAM  FITR  (3,99) 

0003  C 

0004  C   APPROXIMATION  PROGRAM  TO  FIT  A  SURFACE  THROUGH  DATA  POINTS  Y, 

0005  C   UHICH  DEPEND  ON  THREE  INDEPENDENT  PARAMETERS  XI,  X2  AND  X3 . 

0006  C 

0  0  07  C   AUTHOR:  HANS  ZEBNER 

0008  C   DATE   :  AUGUST  26,  1980 

0009  C 

Q010  *,  Hans  Zebner :      4D-Approxinat ion ;       Y»f (XI ,X2,X3) . 

0012  COMMON  /  DATA4  /  Xi,X2,X3,Y 

0013  COMMON  /  AB4    /  A,B 
0014 

0015  REAL  X1(S,S,5) ,X2(S,S,S) ,X3(S,S,S) ,Y(S,S,S) 

0016  REAL  A<64,64),B(64) 
0017 

0018  REAL  C0EF(S,5,5) 

0019  INTEGER  IDC8(  144) , IFILE(3) ,NOCR ( 2> , ICLR (3) 
0020 

0021  DATA  NOCR  /0000333 , 040433B/ ,  ICLR  /015524B , 01S51SB, 006S37B/ 
0022 

0023  101  FORMAT  ("  ENTER  NPNTSi  NPNTS2  NPNTS3  "2A2 

0024  *) 

0025  102  FORMAT  <•  ENTER  L  M  N  "2A2> 

0026  103  FORMAT  <"  ENTER  DATA  FILE  NAME  "2A2> 

0027  104  FORMAT  (3A2) 

0028  131  FORMAT  <•  ENTER  IPRINT  "2A2) 

0029  132  FORMAT  <•  OUTPUT  RAU  DATA  YES  OR  NO  "2A2) 

0030  10S  FORMAT  </•  XI") 

0031  106  FORMAT  (/1X.I2".  DATA  SET"/) 

0032  107  FORMAT  (3X,i2I10) 

0033  108  FORMAT  ( IX, 12, 14F10 . 3) 

0034  109  FORMAT  (/"  X26) 

0035  110  FORMAT  (/"  X3") 

0036  ill  FORMAT  </"  Y") 

0037  149  FORMAT  C"((3A2))) 

0038  1101  FORMAT  ("HANNES") 

0039  C 

0040  C   READ  RAW  DATA  FROM  FILE 

0041  C 

0042  01  URITE  (1,101)  NOCR 

0043  READ   (1,   *)  NPNTSI ,NPNTS2 ,NPNTS3 

0044  URITE  (1,149)  ICLR 

0045  NMAX=5 

0046  IF  (NPNTSI .GT.NMAX)  GO  TO  01 

0047  IF  (NPNTS2. GT.NMAX)  GO  TO  01 

0048  IF  (NPNTS3. GT.NMAX)  GO  TO  01 

0049  02  URITE  (1,102)  NOCR 

0050  READ       (1,       *)    L.M.N 

0051  URITE    (1,149)    ICLR 

0  0S2  IF  (L*M*N.GT.64)  GO  TO  02 

0053  03  URITE  (1,103)  NOCR 

0054  READ   (1,104)  IFILE 

0055  URITE  (1,149)  ICLR 

0056  31  URITE  (1,131)  NOCR 

0057  READ   (1,   »)  IPRINT 

0058  URITE  (1,149)  ICLR 

0059  URITE  (1,132)  NOCR 

0060  READ   (1,104)  IDUM 

0061  URITE  (l.i49">  ICLR 

0062  CALL  OPEN  ( IDCB, IERR , IFILE) 

0063  IF  (IERR.LT.0)  GO  TO  03 

0064  IL*2*NPNTS1*NPNTS2*NPNTS3 

0065  CALL  READF  < IDCB, IERR . XI , IL , LEN , 1 ) 

0066  IF  (IERR.LT.3)  STOP  0001 

0067  CALL  READF  ( IDCB , IERR . X2, IL , LEN , 54) 

0068  IF  (IERR.LT.0)  STOP  0002 

0069  CALL  READF  ( IDCB .IERR ,X3„ IL, LEN , 107) 

0070  IF  (IERR.LT.0)  STOP  0003 

0071  CALL  READF  ( IDCB .IERR .Y , IL , LEN , 160 ) 

0072  IF  (IERR.LT.0)  STOP  0004 

0073  CALL  CLOSE  (IDCB.IERR) 

0074  IF  (IERR.LT.0)  STOP  000S 

0075  IF  (IDUM.NE.2HYE)  GO  TO  17 

0076  URITE  (6,105) 
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PAGE  0002   FITR    3:47  PM   FRI.,  26   SEP.,  1980 

0077  DO  04  Kl  =  l  .NPNTS1,1 

0078  WRITE  (6,l66)  Kl 

0079  URITE  (6,107)  < K3 ,K3=i ,NPNTS3 , i ) 

0080  DO  04  K2=l .NPNTS2.1 

0081  04  URITE  (6,1<$8)  K2,  (  XI  (K  1  ,K2,K3>  ,K3»i  ,NPNTS3  , 1) 

0082  WRITE  (6,109) 

0083  DO  OS  Kl=l .NPNTSi,i 

0084  URITE  (6,106)  Kl 

0085  URITE  (6,107)  (K3 ,K3=i ,NPNTS3 , 1 ) 

0086  DO  OS  K2=l.NPNTS2,i 

0087  OS  URITE  (6,108)  K2, (X2(K1 ,K2,K3) ,K3»i ,NPNTS3 , i ) 

0088  URITE  (6,110) 

0089  DO  06  Kl  =  l .NPNTSl,i 

0090  WRITE  (6,106)  Kl 

0091  URITE  (6,107)  (K3 ,K3=i ,NPNTS3 , i ) 

0092  DO  06  K2»i .NPNTS2.1 

0093  06  URITE  (6,108)  K2, (X3(Ki ,K2,K3> ,K3»1 ,NPNTS3, i ) 

0094  URITE  (6,111) 

0095  DO  07  Kl=l ,NPNTSl,i 

0096  URITE  (6,106)  Kl 

0097  URITE  (6,107)  (K3,K3»1 ,NPNTS3, i > 

0098  DO  07  K2=l,NPNTS2,i 

0099  07  URITE  (6,108)  K2, ( Y(K1  ,K2 ,K3)  .K3=l .NPNTS3 , 1 ) 

0100  17  CALL  MAT4  (NPNTSl  ,NPNTS2,NPNTS3,  LjM^COEF,  IPRINT) 

0101  CALL  CODE 

0102  URITE  (IFILE.liOi) 

0103  CALL  OPEN  ( IDCB , IERR , IFILE, IOPTN, 0 ,26, 144) 

0104  IF  (IERR.LT.0)  STOP  0006 

0105  CALL  URITF  ( IDCB , IERR  ,COEF ,2S0 ) 

0106  IF  (IERR.LT.0)  STOP  0007 

0107  CALL  CLOSE  (IDCB, IERR) 

0108  IF  (IERR.LT.0)  STOP  0011 

0109  STOP  0077 

0110  END 

t  FTN4  COMPILER:  HP92060-16092  REV.  1926  (790430) 

«*   NO  UARNINGS  **   NO  ERRORS  **    PROGRAM  ■  012S9       COMMON  -  00000 
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U.6.2.   Load  map  of  3ample  user  program  FITR 


FITR 

10042 

12414 

Hans  Zebner:      40-Aooroxiininon;   Y=f  (  XI  .  X2  ,  X3)  . 

to  load  these  progr 

DATA4   12415  14364 
MAT4    14365  20251 
MAT41   20252  20666 
MAT42   20667  21066 
S4      21067  21370 
AB4     21371  41570 

am  nodules,  enter  (from  LOADR):   HS.JTPLBL 

4D-Aqpr oxination   /  DATA4  / 

4D-Approxj.naTi.on  arranqe  systen  natnx  ind  vector. 
4D-Approxj.nati.on        arrange  subnatrix  1.  class. 
4D-Approx mation        arranqe  subnatrix  2.  class. 
4D-Approx mation                conpute  sunnations. 
4D-Aqpr oxinat ion   /  AB4    / 

LOGLU 

41571 

READF 

41647 

OPEN 

42611 

CLOSE 

43107 

CLRIO 

43317 

CURD. 

43326 

JSMUE 

43327 

LURQ 

43416 

.DADS 

44001 

.DMP 

44111 

.DDI 

44257 

SESSN 

44560 

R/U* 

44576 

P.  PAS 

44735 

.DNG 

44764 

PAUSE 

44774 

♦ALRN 

45075 

FMTIQ 

45213 

ERRO 

46512 

ABS 

46602 

.DDE 

46611 
46623 

.DIN 

.RTOI 

46631 

.FPUR 

46725 

.SBT 

46767 

.  .FCM 

47030 

PAU.E 

47045 

ERO.E 

47046 

♦OPEN 

47047 

RW1UB 

47224 

RUND* 

47572 

FRMTR 

47715 

FMT.E 

S3353 

REIO 

53354 

RMPAR 

53501 

PNAME 

53544 

LUTRU 

53612 

*SETP 

53721 

.CFER 

53746 

.LBT 

54024 

41646 
42610 
43106 
43316 
43325 
43326 
43415 
44000 
44110 
44256 
44557 
44S75 
44734 
44763 
44773 
45074 
45212 
46511 
46601 
46610 

46630 
46724 

46766 
47027 
47044 
47045 
47046 
47223 
47571 
47714 
53352 
53353 
5350  0 
S3543 
S3611 
53720 
53745 
54023 
540S4 


92067-16268  REV. 1903 

92067-16125  REV. 1940 

92067-16125  REV. 1903 

92067-1612S  REV. 1903 

750701   24998-16001 

92067-1612S  REV. 1903 

92067-16268  REV. 1903 

92067-16268  REV. 1903 

780818   24998-16001 

730818   24998-16001 

781021   24998-16001 

92067-16125  REV. 1903 

"2067-16125  REV. 1903 

92067-16125  REV. 1903 

780818   24998-16001 

771122   24998-16001 

92067-16268  REV, 1903 

24998-16002  REV. 1926 


771122 
750701 

780818 

780818 

780921 

781106 

770518 
750701 
750701 
750701 
92067-16125 
92067-1612S 


24998-16001 

24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
24998-16001 
REV. 1903 
REV. 1903 


790228 
790719 
790215 
781229 

780526 
790202 
790223 


730413 
781214 
740801 


77071S 
790417 


790103 
781003 
92067-16125  REV. 1903  780801 

24998-16002  REV . 1926  790503 

24998-16002  REV. 1901  781107 
92067-16263  REV. 1903  790316 
781106   24998-16001 

771121   24998-16001 
92067-16268  REV. 1903  790223 

781106   24998-16001 
750701   24998-16001 

770518   24998-16001 


20  PAGES  RELOCATED 
LINKS :BP     PROGRAM: LB 
/LOADR: FITR    READY  AT 

/LOADR : SEND 


20  PAGES  REQ'D 


LOAD:TE 
3:53  PM 


NO  PAGES  EMA 


COMMON :NC 
FRI.,  26   SEPT,  1980 


NO  PAGES  MSEC 
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4.6.3.   Results 


XI 


printed  output 


1 . 

DATA  SET 

I 

2 

3 

4 

s 

1 

.100 

.100 

.100 

.100 

.100 

2 

.100 

100 

100 

too 

too 

3 

.100 

.100 

.100 

.too 

.  too 

4 

.100 

100 

100 

too 

too 

S 

.100 

.100 

.100 

.100 

.100 

2. 

DATA  SET 

1 

2 

3 

4 

5 

1 

.200 

.200 

.200 

.200 

.200 

I 

.20  0 

.200 

200 
.200 

200 
.200 

138 

20  0 

.200 

4 

.200 

200 

200 

200 

200 

5 

.200 

.200 

.200 

.200 

.200 

3. 

DATA  SET 

I 

2 

3 

4 

5 

1 

.300 

.300 

.300 

300 

.300 

2 

.300 

300 

300 

300 

300 

3 

.300 

.300 

.300 

300 

.300 

4 

.300 

300 

300 

300 

300 

5 

.300 

.300 

.300 

300 

.300 

4. 

DATA  SET 

i 

.400 

.400 

.40  0 

4 
400 

400 

2 

.400 

400 

400 

400 

400 

3 

.400 

.400 

.400 

400 

.400 

4 
5 

Ml 

400 
.400 

Ml 

333 

138 

5. 

DATA  SET 

1 

2 

3 

4 

5 

i 

.SCO 

.S00 

.500 

500 

.500 

2 

.500 

500 

500 

500 

500 

3 

.500 

.500 

.S00 

500 

.500 

4 

.S00 

500 

500 

500 

500 

S 

"3 

.500 

.500 

.500 

500 

.500 

C 
L. 

DATA  SET 

I 

2 

3 

4 

5 

i 

-20.000 

-20 

.000 

-20 

.000 

-20 

.000 

-20 

.000 

2 

-10.000 

-to 

000 

-10 

000 

-10 

000 

-to 

000 

3 

0.000 

0 

.000 

0 

.  000 

0 

000 

0 

.000 

4 

10.000 

10 

000 

10 

000 

10 

000 

10 

000 

5 

20.000 

20 

.000 

20 

.000 

20 

.000 

20 

.000 

2. 

DATA  SET 

1 

2 

3 

4 

S 

1 

-20.000 

-20 

.000 

-20 

.000 

-20 

.000 

-20 

.000 

2 

-10.000 

-10 

.000 

-10 

.000 

-to 

000 

-10 

000 

3 

0.  000 

0 

.  000 

0 

.  000 

0 

.  000 

0 

.  000 

4 

10.000 

10 

.000 

10 

.000 

10 

000 

to 

000 

5 

20.000 

20 

.000 

20 

.000 

20 

.  000 

20 

.000 

3. 

DATA  SET 

1 

2 

3 

4 

5 

1 

-20.000 

-20 

.000 

-20 

.OJJO 

-20 

.  000 

-20 

.000 

2 

-10.000 

-10 

.000 

-10 

.  000 

-to 

000 

-to 

000 

3 

0.000 

0 

.000 

0 

.000 

0 

.  000 

0 

.000 

4 

10.000 

10 

.000 

10 

.000 

to 

.000 

10 

000 

5 

20.000 

20 

.000 

20 

.000 

20 

.  000 

20 

.000 

4. 

DATA  SET 
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i 

-20 

000 

-20 

000 

-20 

.000 

-20 

000 

-20 

000 

2 

-10 

000 

-10 

300 

-10 

000 

-to 

000 

-10 

000 

3 

a 

000 

0 

000 

0 

.  000 

0 

000 

0 

000 

4 

to 

000 

10 

000 

to 

000 

to 

000 

10 

000 

s 

20 

000 

20 

.000 

20 

000 

20 

000 

20 

000 

s. 

DATA 

SET 

i 

-20 

1 
000 

-20 

2 
000 

-20 

3 
000 

-20 

4 
000 

-20 

000 

2 

-10 

000 

-10 

000 

-10 

000 

-to 

000 

-10. 

000 

3 

0 

000 

0 

000 

0 

000 

0 

000 

0 

000 

4 

10 

000 

10 

000 

10 

000 

10 

000 

10 

000 

5 

20 

000 

20 

000 

20 

000 

20 

000 

20 

000 

X3 

• 

L. 

DATA 

SET 

I 

o 
c 

3 

4 

5 

i 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

2 

-20. 

000 

-10. 

000 

0 

000 

to 

000 

20. 

000 

3 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

4 

-20. 

000 

-10. 

000 

Q 

ooo 

10. 

000 

20. 

000 

5 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

2. 

DATA 

SET 

I 

3 

4 

s 

1 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

2 

-20. 

000 

-10 

000 

0 

000 

10 

000 

20. 

000 

3 

-20 

000 

-10 

000 

0 

0  00 

10 

000 

20 

000 

4 

-20. 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

5 

-20 

000 

-10 

000 

0 

000 

10 

000 

000 

3. 

DATA 

SET 
1 

2 

3 

4 

5 

1 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

2 

-20. 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

3 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

4 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

5 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

4. 

DATA 

SET 

1 

-20 

1 
000 

-10 

1 
000 

0 

3 
000 

10 

4 

ooo 

20 

000 

2 

-20 

000 

-10 

000 

0 

000 

to 

000 

20 

000 

3 

-20 

000 

-10 

000 

0 

.300 

10 

000 

20 

000 

4 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

5 

-20 

000 

-10 

000 

0 

.000 

10 

000 

20 

.000 

5. 

DATA 

SET 
1 

2 

3 

4 

5 

i 

-20 

000 

-10 

.000 

0 

000 

10 

000 

20 

.000 

2 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

3 

-20 

.000 

-10 

.000 

0 

.000 

10 

000 

20 

.000 

4 

-20 

000 

-10 

000 

0 

000 

10 

000 

20 

000 

s 

-20 

.000 

-10 

.000 

0 

.1700 

10 

.000 

20 

.000 

1.    DATA   SET 


I 

2 

3 

4 

5 

i 

.400 

.700 

.800 

.700 

.400 

2 

.450 

.750 

.850 

.750 

.450 

3 

.500 

.900 

.900 

.800 

■  S00 

4 

.450 

.750 

.850 

.750 

.450 

5 

.400 

.700 

.800 

.700 

.400 

2. 

DATA  SET 

i 

2 

3 

4 

5 

1 

.463 

.763 

.863 

.763 

.463 

2 

.S13 

.813 

.913 

.813 

.513 

3 

.563 

.363 

.963 

.363 

.563 

4 

.513 

.813 

.913 

.813 

.513 
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5 

.463 

.763 

.863 

.763 

.463 

3. 

DATA  SET 

1 
2 
3 
4 
5 

1 
.650 
.700 
.750 
.700 
.650 

2 

.950 

1.000 

1.  050 

1.000 

.950 

3 
1.050 
1.100 
1.150 
1.100 
1.0SO 

4 

.950 

1.000 

1.050 

1.000 

.950 

S 
.650 
.700 
.750 
,700 
.600 

4. 

DATA  SET 

1 
2 
3 

4 
5 

1 

.963 

1.013 

1.063 

1.013 

.963 

2 
1.263 
1.313 
1.363 
1.313 
1.263 

3 
1.363 
1.413 
1.463 
1.413 
1.363 

4 
1.263 
1.313 
1.363 
1.313 
1.263 

1 

1 
L. 

S 
.963 
,013 
.063 
,013 
.963 

5. 

DATA  SET 

i 

3 
4 

s 

1 
1.400 
1.450 
1.500 
1.450 
1.400 

2 
1.700 
1.750 
1  .800 
1.750 
1.700 

3 
1.800 
1.850 
1.900 
1.850 
1.800 

4 
1.700 
1.750 
1.800 
1.750 
1.700 

1 
1, 

1 

1, 

1 

5 
.400 

450 
.500 

450 
.400 

Co«ffici«nts  COEF(I,J,K> 


K  »  1  2  3 

J  *  1  -.418E+04  . 96SE+03  .691E+02 
J  »  2  -.134E+01  -.39OE-01  .282E-02 
J  -  3   . 121E+02  -.284E+01  -.203E+00 

K  »  1  2  3 


K  ■  1          2  3 

J  »  1   .4S8E+0S  -.106E+0S  -.757E+03 

J  =  2   .143E+02   .417E-f00  -.300E-01 

J  =  3  -.133E+03   .312E+02  .222E+01 

K  =  t          2  3 

I  »  3 

K  =  1          2  3 

J  »  1  -.994E+0S   .230E+0S  .165E+04 

J  »  2  -.238E+02  -.696E+00  .499E-01 

J  =  3   .289E+03  -.678E+02  -.484E+01 

K  =  1          2  3 
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5.   CONCLUSIONS  AND  RECOMMENDATIONS 

Compared  to  the  probe  calibration  approach  described  in  Ref  3. 
the  method  presented  here  provides  a  much  simpler  way  to  both  calculate 
the  probe  calibration  surfaces  from  measured  data  and  to  apply  the 
calibration  to  on-line  data  reduction.   The  iteration  required  in  the 
method  of  Ref  3  is  completely  omitted.   It  is  noted  however,  that  the 
Gauss  Jordan  Elimination  lacks  some  sophistication.   The  system  matrix 
for  example  undergoes  the  elimination  without  prior  conditioning. 
Particularly  in  the  case  of  the  3-D  and  the  4-D  approximations,  numerical 
round  off  errors  have  an  influence  on  the  accurracy  of  the  coefficients. 
Consequently  it  is  recommended,  that  the  Gauss  Jordan  Elimination 
routine  be  revised  to  use  double  precision  constants  and  a  routine  that 
conditions  the  system  matrix  be  added.   Time  constraints  made  these 
steps  impossible  for  the  author. 

Following  these  refinements,  a  calibration  procedure,  similar  to 
McGuire's  (ref  3)  should  be  formalised.   Since  the  4-D-  approximation 
requires  large  arrays  and  therefore  extensive  CP  -  memory,  extended  memory 
access  (EMA)  is  necessary  if  the  HP  21MX  computer  is  used  to  perform 
the  calculations. 
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APPENDIX  A:   Some  Useful  Matrix  Conventions  and  Operations 
Al.   Submatrix  Notation 


A- 


Om     °,z    atI    <*,4 


Qz,     °&    azs    *u 


%     art     ***    **» 

°4,       Qai       Q+l      *4* 


II* 


*** 


Om     <2+t    <?«*    *W  a^^r 


&*r- 


K     A* 


A 


lr» 


a: 


/**»» 


defined  through  elements  a, .  or  submatrices  A. . 

ij  ij 


A2.   Diagonal  Lines  and  Diagonal  Bands  1.  order 


Diagonal  lines  1.  order  and  Main  diagonal  line  1.  order 
Diagonal  bands  1.  order  and  Main  diagonal  band  1.  order 


A  - 


*   Diagonal  lines  1.  order 
w.r.t.  elements 
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>  y> 


A 


*     Diagonal  bands  1.  order 
w.r.t.  subraatrices 


*  ** 


A3.   Diagonal  Lines  and  Diagonal  Bands  2.  order         0($>  ^ 
Diagonal  lines  2.  order  and  Main  diagonal  line  2.  order 
Diagonal  bands  2.  order  and  Main  diagonal  band  2.  order 


A  - 


Main  Diagonal  Line 

1.  order  w.r.t.  elements 


*  Diagonal  lines  2.  order  w.r.t.  elements 
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** 


Main  Diagonal  Band 

2.  order 

w.r.C.  Submatrices 


*  Diagonal  bands  2.  order 
w.r.t.  submatrices 
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APPENDIX  B:   Software  Description:   Flow  Charts 

The  following  are  given: 

Flow  Chart  Page 

MAT2  93 

52  97 
MAT3  98 
MAT31  106 

53  109 
MAT4  111 
MAT41  119 
MAT42  124 

54  127 
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7ouv   eAcW        MAT 2 


(entry 


1 


IT/0    /4«Sta^rvicvnh 

L0=  £ 


YES 


2.1.   a*d   2.2. 
wi-Bi  actual  para- 


0/L3- 


1 


Sort  ^fttztl 
edge  Sccfro* 

iTCtjRl  •  0 

Li--2*L-I 


/ 


93 


OO  04  K*l,U,l 


0 


L 


To 


op 


i 


£*fne    location 
of element 

7*X 


o3 


t 


Function  SZ 
aincf  H"-//e   value 
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LO-C 


YS$ 


£ 


Ul'spfay  e<juaf*ons 
4J  ancC   d-2. 
U//3A   GtdoMxL 
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Sut*nfa+rt.x    A  c/ass 
(SvQMl)   iASt»j 
SUSGOV  TINE  M**l 


II' 0 


Q 


4.-12 


DO   04-  £> TSTART.JSTOP,  I  jL-j. 
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APPENDIX  C:   Software  Description:   Listings 
The  following  are  given: 

Listing  Page 

Two-Dimensional  Approximation  132 

BLOCK  DATA  AB2  132 

BLOCK  DATA  DTA2  133 

PROGRAM  DEM02  134 

SUBROUTINE  MAT 2  137 

REAL  FUNCTION  S2  140 

INTEGER  FUNCTION  IEL22  141 

SUBROUTINE  IEL2  142 

REAL  FUNCTION  F2  145 

Three-Dimensional  Approximation  146 

BLOCK  DATA  AB3  146 

BLOCK  DATA  DATA3  147 

PROGRAM  DEM03  148 

SUBROUTINE  MAT3  151 

SUBROUTINE  MAT31  155 

REAL  FUNCTION  S3  156 

INTEGER  FUNCTION  IEL33  157 

SUBROUTINE  IEL3  158 

REAL  FUNCTION  F3  161 

Four -Dimensional  Aooroximation  162 

BLOCK  DATA  AB4  162 

BLOCK  DATA  DATA4  163 
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2S  •  »£    SREATLT    A»ARQ60    *T    cASIER/JS*      ANO    FASTI?    Rt(iS«AMMt"»««///| 

103  FORMAT  (2**      OATA    >3INTS    SENERAreO.) 

10*  FORMAT  U*X,5HX1.1tN.oX,i*»l*AX.  7«.*nrMIH,  :i,«rirnAi/4A,  7-n>l3TTE>  »<•< 
•U.  F10.3I/<«X,<.'*JSE»,*<1<»F10.3I  ) 

SO  103  F19BA;  («.7h      ".jrr  =  5    mo    JSE'    AR«S    OEFt'lCO    ANO    AXES    5R4**.! 

196  FORMAT  (20H      OATA    POINTS    }»aj*.I 

107  FORMAT  (3*M      51TER    }30ER    JF    '0ir*-3NIAL    (HaR.0R.lJ,     1RRINT    »-<0    ISlATl 

108  FORMAT  (22H      CURVE    PIT    :3l»L£TED.) 

109  format   (iSH    :u*vc  fit  ora**.i 
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OEFiNf   »iarrt^   »*o  jse*  *»e»Si     oraw  axes. 
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10 
IS 
20 
23 
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.0 
♦3 
30 
33 
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63 
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S3 


7*/l73      J»T«2      P""ON*  FTH    ».a»<.»4 

subroutine  n»i2  cnMTu.i.»caEf,ii»»  hit.istat.ierror  i 


a  C  / 1  2  / : 


PRESET    SYSTEN    NATRIX    4    »M3    <IShT    tAHD    SIOE    VECTOR    J    F  )R    4 

2!}-»Qirsapt»t  i^»«jJin»ri3N.     :h£  .iiKid  jau*iioi*  system  then 
H   sjl.eo  »r  jsi^i  »   muss   joroam  cuutiin". 


AUTHOR! 
OATEi 

VitlillC 

L 

COEFIT) 

IPRINY 


ISTAT 


moot 


(1(2)0) 
YU36I 


HANS    :E4NER 
3ECEN9ER    2*. 


1439 


Em  an*  no* 

•iui9£*  j?  5»rt  points 


NUNBER    0*    .. 

3Hi     ry  l*r  n        - 


COEFFICIENTS 

-.13*  'oCrwaHii 

POLTN'JNIAl    ^EFFICIENT? 


int.s'ua 


CQHTRfl 
•2    PR" 

4N3    4PfER 


•I 


»0    PR 


siinination. 
print  eojati 
the   jauss    J3    . 
PRINT    JASIC    53UAT 
ACTUAL     PARAIETIrS 

no  output   »t  all. 


FOR    THE 
L    •<OROER< 


,rr  of   pr;nt  iut 

THE  JAJSS  JGPPAN 


1  I 


TYRE 

INTEieR 

real 

INTEGER 


ON  SYSTEM  AFTE 

no*"*  9Lit:'i/.r: 


IJMS     •'ITH 


ON. 


.STSICS    TO   it    PRINTEO 

print  :*5iNAT|i  :?  iAfA 


"«rii:J 


integer 


«o 

control  snrrsi 

/ISTAT/.2   prini    - 

ANO    APPROXIMATED    "3IST3 

ANO   OEVIATIO'IS     (EPR3RI. 
/ISTAT/«1    KIT    PRINT    E>RflP*, 
ISf *r>0        SET    »DSOLyT£    SP.a.n*. 
tSTAMO         GET    RELATIVE    ir'«l!», 

error  fug  *eturnej   ro  the  calling 

PR3GRA1. 

•♦i  no  git  to*,  oetecteo,     oooo  JOS. 
•-I  M  cR*qr  occjreo  In  Tir   :aos:ia« 

ELIMINATION;       COEFFICrZHTS    ARE    MOT 
ACCUKATE.       REVIE*   PAP.A1ETERS. 

XI-C3SR0INATES    OF    THE   OATA    POINTS 
Y-CdOROINAUl    JF    THE    DATA    POINTS 

THE    OATA    POINTS    XI    1NQ    r     lAVe    TO    Si    PR0V13E0    Ti 
THE    LASELSO    C3NN3N    HOCK    OTAJ    4*0    The    5TSTE1 
NAT»U    »    4N3    T-if    RIGHT    -tAMP    JlOt    rfiCTOR    9    ARE 
AVAILABLE    THROUGH    ftE    lASELEJ    CC.1S0N    3L3C<    432. 
THIS    SAVES    A   L3T    3F    THAT    Pf.ECIOUS    THING   CALleO 
NENORT. 


REAL 
<EAL 


101 
SOI 
602 

60S 
606 
60  7 
60S 
609 

■Hi 

616 
613 
616 


COntON   /    0TA2      /    «l.r 
CORAOH    /    AS2        /    ..i 

REAL    »1.236J#Tt2J6l 

SEAL    :3lF(T).ST4r(256» 

FOaNAlS    tAT2    STA6T 
FO»«AT    (///7H      »FTE*i 


I3.*IH    STEPS    NO   SaUUTION    <rTHIM    THE    »CCU«ACT 
T./) 

JF    APP^OItNATtn^    POLTNONIAL    IN    U    0I«ECTI0n» 
5AT4    »3tNTS    <ITH    CONSTAHT    <1>16<«I3/I 
N    1AT*U    A    ANO    /ECT3K    S    SEFuBE     iAUSi    JOIiOAN    i 


•0F,£10.3f 11H    (        I     JUt 

FoaNAf    i///«ih     3aoE* 

•5X»li/»0H        NU«1E*     3P 
F3AN4T     [///Din       STSTE 

•LININATION/I 
F0HP1AT    (><.!. 13taXt  [2)11 

FCRNAT  i2X»I2»i3UX.l»p.2)/((AX,13(«X.:»  =  6,4))  I  I 
PORKAT  l///4aH  ;5uaT;3h  iTSTEfl  4?TE?.  'JuSs  JOHOA 
FOBNAl    I///22H      :3EFF[;iEtrs    C3EMIJ/> 


AN    EL1NCNATI3N/) 


soa«4r   (2x»3Ht   ••t2»t 

FO«NAT    ((5(2X,I3»1«.; 

:3BN4T     I///4AH        jJOIl 

SOIHAT  W//3CH  44S3L 
FO»H»f  (///3JH  »EU'" 
rORNAf    ((6I2X>Z3>1X.2 

fornat  (///iirt    sT4r: 

♦2/7«.19MNAXi*UH   ERROR 
-ORNAT    (///21H      JTATI 

•2»2H     t/7X»  I'S-'NAXt  1UN 
FORNAT    (///7H       AFTER, 


X»l»s9.2I 

4TJ3    3F    3ATA    POINTS      - 


as  ;w{n/i 


•2aRSi?s3;iTS  Itop  3ul 


pe<».2  I  >) 

4fci  as  oata  porsrs  -  as  approxxnateo/) 

JTE    =»R3R    AT    ;AC->    POTNT/1 

IVE    6RR3R    4T    EACH    POINT    IN    './)  -> 

PF6.il 1 1 

:MC5     .  ;»«»»n  //7X  .;3HAVERAGE    ;R«3R»  10I»1P£9. 

AT     I'.i  3.  L-p  lf»€9.;I 
STICS    ;j«NARri//7Xt;'rtAyERAGE    ER«0*»iOi#2PP'». 
ERROR    4T    [•.Ii#ix»i?Pt'.2*2H    S) 
!3«*1H    STEPS    NO   33LUTIQM    .1THIN   THE    ACCURACT 
T./l 
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subroutine  i»fi  r*/i7»     jpt«2     >*o*r  ?tn  *.d»-.<» 


no 


115 


120 


123 


1»J 


no 


LI'LOSLUI  ISESSU1 
L0»« 


.    •« 
EERR0R>1 

ITrnax.50 

<WJ  tP5lUN...wCE-O07 


PRINT    APPROXIMATION    »t«AN6rs»S.     IF    IPRIHT    IS    ;»HTc«    KM    0. 


N0R0Rl»l-l 

♦«  if  <iprint.lt. 01   ;a  TO  01 

WRITE    1L0.    6011    iJR0RL,NPNTSl 


START    TO    PRESET  EOGc    SECTI3N    ELEMENTS. 


100  31    IPOW»1-0 

|*2«I>-1 


J' 


HljJ 


I?    (ll.6T.LI    S3    TO   02 
109  TOF-   EOCf    SECTION   :LENcNTS. 

1*1 

la*  TO    03 

RISNT  1ANQ  EOCc  SECTION  ELEMENT!. 
02  !•<«»;)-(. 

J«t 

PRESET    E00C   SECTION   ELEMENTS    IN    57STF.N    1ATRIX    JJINS    St. 


03  a<i#ji.S2<nrntsi»ipowri.oi 
a*   iPOwRWPom*! 

copt  oefineo  elenents  jisonallt. 


Li«L-l 

oo  as  1.1,11,1 
oo  ;j  i«i.i./i 


05    ill  HiJ-i  )  •»(  I.J) 

w  PRESET  RIUmT  1AN0    SlOe    VECTOR    S. 

130  90  OS    l»l#l,i 

0*    IUl»S2<NPNTSi»lP0*Rl,ll 

PRINT    SYSTEM   MATRIX    »N0    RI0NT    N4.N0    SI5t    VECTOR    S.    IF    IMUT 
135  C  IS   SET    TO    2. 


IF    IIMINT.Nfc.2l    SO    n   o» 

write  (i- 


wRITE  <LQ»    6w2) 

WRITE  <LU»    603)    (J.J-l.t.D 

1*0  oo  or  i'I.cji 

07    «*ITe  <L0»     50»)     I.  (»l  I,  J).  J«l,t.l).3tll 

•  RITE  iij.    aojl    U.J.i.l.K 
IPRINW 

1*5  !         .    JAUSS  JOROAN   ELIMINATION. 


C 


OR 


N1«L 

6SSL0N.€PSLCiN»Nl 

CALL     IEL2     <Nl,r:=R,ir»NAX,EPSl1NI 


150  IF    (ITER. LT. CI    SO    M    Zl 

C  PRINT    EQUATION    SYSTEM.     IF    [PRINT    IS    SET   TQ    1. 

IF    < JPtIHT.Nf.il    S3    T3    10 
155  WRITE    <L0»    oCil 

write   ua.   son  ( j.j.i.Mi.n 

00    09    I.l,»l,l 
09    WRITE    <IJ»    6041    I.UH.JI.J'I.iI.IIfJCI 
WRITE    (L0#    6031    U»J«l»Ml»ll 

C  0EFIN6   COEFFICIENTS   ANO   PRINT   then*   :p    IPtINT    IS    GREATER    IMAM   0. 


c 


10    CONTINUE 

00    11     L.-l.l.l 

IX  ::efu;i.j(l;i 

If    (IPRINf.i.T.01     JO    TO    13 
WRITE    (L0»    6061 

12    28x11    ("ioi'wTI    L1,C3EFUU 
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\ 


C  STATISTICS*    IF    tSTAT    AUTTHIHS    cl.SE    SKT    U, 

If  fisTAr.LT.oi  :e»rj»-i 

ITS  IP    t  ISTAT.4T.0I     r.ErUJR.2 

iSTAl«lA8S<ISTAll 

C  CONFUTE    48S01UTE    ERR'IR. 

C  C    »t »*F«uiI«»f,0    tfALJc      -      SXVEH    »»tuel     ! 

OU   1*    I'l.MRNTll.l 
14    SfATID'frZtCOEF.XKII/HOHOiD-f  (  II 

2  CONFUTE    RELATIVE    ERRJR.     IF    [STAT    <4S    LESS    THkH    ISRO.    I.s.     It**JR«l. 

us  i  t  •«Ap*taxinAreu  <*tje     •     jive*   *alue  i  /(givcn   <«ljei   j 

if  1 1 E»«aa  .Nt.;j  so  ro  i6 

SRROR    3TATI$TlCSi       CJI^UTE    AVERAGE    ERR^R    4*0    fli»0    LOCATION 
C  OF    KAIIXUH    eiKGR. 


ISO 


190 


16  •«!»«•.» 

195  c«»H4»«A8S(STAT(ill 

if  7a8sistaT(i  )T.(.£.cr*"a»»   ;a  to  17 

eRiHAj.AiS  UTAH  ;  )  I 
200  UU.I 

17  6RRAVe«ERRAVe»»8S<STATtI)) 


239 


;rravc>c  Jj»vt  /  -i»m  rsi 

IF    (I- 


.     USt»!.»i.ii    40      1    IS 

20)  C  "It"    3R0INATSS    3F    3aTa   "3INTS    4i»0   AMRflXIllATEO    VALUES. 

C  IF    /ISIAI/    IS    EQUAL    2. 

<»IT5     (10.  60S) 

fclltfc     (LO,  5-j>l     (l»r(I>»I*l*MF«TSl/l) 

210  •RITE    (Li.  61C.) 

«RITE     lUi  <*»)    (I,F2(C3EF,<lll),  «0RPR1),  I*l/KPNTS1»1> 
ISTAT-l 

C  FRINT    ERRORS*     IF    /tSTAT/    its    1    3R    2. 


219  C 


220 


la  if  (isT4r.ME.il  S3  ro  :o 

If  iltftijR.ea.il    ;a  rj  .1 

'•INT  ABSOLUTE  ERRORS. 


-RITE  (L0«  6111 

*<UTS  U3.  »C<»)  (I*JTAT(II.I.l».i»HTSl*l) 

40  ro  lU 

229  C      »RI*T  RELATIVE  ERRORS. 


C 


19    .aire    <L0.    612) 

■RITE    nj!    al3>    iI.srATili.  [.i,-H»HTti,;i 


230  C  MINT    ERROR    STATISTICS. 


20  IF    i  IERROR.cO.  J)    4RITE    <L3»    61V)     i  **wr  ,  L  -t  1 .  c  *  *  »4i 

IF    I  llRRQR.Ea.l)    «RITE      L3«     613)     :  ■Hi*  C  <  I  H> .  •  J  -  'm 
.0    TO    23 

FATAL    ERROR    OETECTED.       3UTRUT    ERROR    SeSSAGE    4N0    PRESET    C3EFFIC  I  i  UTS . 

21  ITER  — ITER 

IERROR  — £  ERROR 


240  8(1 1  •Ml) 

90   22    IC«2»N4.1 
22    KICJ..C 

•  RITE    (LI»    1C1I    ITER.?»SL3N 
VRITE     (LJ.     6161     I TER*  E»  5L3N 
249  SO    TO    10 

C  RETURN. 

230  »  !«UM 
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?*/W:       ;»'•*       »■.;»» 


ni/l  2/ii 


: 


»|4«.  'w«c;:>  :  2  : «•# *  r : ; » :»«••: 


:  r  i 


3! 


♦) 


>£»*:»■  sununons   *:  >»:ic'  •■•*  :.:-£irs  m 
;rsr£i   "»r!;i  1  n3   -»e   -.i*i>?s  'ill    :■  "-€ 
Slot   <£:•-»   ;  ::»   •-=  .:'-i»  ;.-»•:*••  :•:•;-. 


:>    t-it 


MTtji 


lecusu  2»» 


riHULi 


tt 


:j 


SUMS     .'•>. 


"■:'n' 


.jts5   :»   •■!{ 


M2S6I 


r    /iLreS    He    ;-<C«X65    ;i     '-e 
SUWUIIMi 
-»    !i;-ir    ■•its    ::3J    fiZr^t    SLi^.) 
*«SUir    3*    SJMUTIM. 

U-CflOtOPUTfS    3*    P4i    :*T*    »OIMTS 

'-c:-»3:«*rej  3*  >c  s^rt  »3;hts 


TT#s 


tnliii' 

;tre;c« 


:ir«:8« 


iNTjieg 


itU. 


:^»QtTii»ri 


5«r*    J-3JHTS    (1     4H0    T     J>.»-    '3    Je    >»3<I3£}    "i«'J 
.ideCfQ   :3i«a>«    SIX<    3r»2. 


.;•■:*  /  3T*2     /  ii»r 
iiii.    U  Uitl»r(2S»l 
S».0 

:»   :  :»c<i'..es.:)   ;2  T3  n 
:.••.'.•». ; :  ail 


,,     T?    P    32 

:f   t  jr. is  .ci   ir«; 


3J   s>$*Ji»si 

:2«s 

Hi 
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45 


3* 


7A/175        J»T«2      l>10ni» 
NTEiE»     PUNCH')*    I6L2J     (Nl 


rTN    ^•a»-»«»i 


GAUSS  JOROA(  eiIitNATI3N  WITH  *ATRIX 

<««UU:  iLLZ  IS  SET  ' 


*nd  vector  i.     r-tc 
ra  i.   ip  the  £L:n!M»r:u»«  is  cjupttTEO. 
ifcti  . t ss  tmah  o  indicates  a*  e»ro*. 


4(jThOR  t 

OATEi 
VAMlASLfc 


IE12 


4(  7,7) 
9(7) 


5    ISiHti 


HI, 

36Cc»SE*    ?*.    1980 

£x(»tA.>tArn,«  tyre 

MUN8ti»     JF     EQUATIONS  INTelcU 

control  para^ete*.  returmco  ro  r»s  ii»Ts4e« 

CALLING    ROUTINE. 

•i   jAjss   jo«9>'i   ri. i.iiuAr T.nj\   was 

•6RPORN60    CORSeCTLr. 

<o  e^ris   encountered  during  gauss 
joroan  elimination. 

srsrEi  i*r»i<   »  j=*l 

Jli-tT    HANO    SIOE    (tCIM    1  »€AL 


INRJRTANTi    tATRH    A    4N9    VICTOR    S   -iaVC    TO   36    »»CWIO?0    MROuG.- 

'ne  cAatLEo  :3iion  acacK  \ti.     sut-<  »  **u  4  ..III 

•  ILL    UNOE9G3    CHANGES    OURING    EXECUTION*    [•«     Tit 

elXnuuon. 


annaN  /  A82     /  «>a 

EAl     1(7.71,4171 


ill 


oo 
00 

CON 
M 

!P 

91  COM 

3(K 

92  CON 

on 


•a»j 


2 

-030 
#N.I 

»N»i 


98    4(1 

A(? 


9b  I 
92  X 
K.E 

sr  — 

01  J 
•  J|a 

(J.t 

riNu. 

i«a(*»»c3NST«adi 

rings 

l) 

»N»1 
[>4l/C3-IS? 


n  go  ri  g» 

#n,1 
K>W)*C0'4ST«A(I»JI 

I)    4l«.,j»-.   : 


mi'v.i 


11 

J, 

;e« 
1? 

CON 

£0 

IEL 

RET 
ENO 


J 
.J  1 
»!)• 
)•»( 
C»  J 
ROR" 
»  < 
A8S. 
«0«« 
T[NU 

1  us 

iif 
r  inu 
ra  0 
i—i 
zi-l 

URN 


41 
I. 
I) 
•I 

o 

•I 

(4 

IE 

fa 

IE 

»J 
E 

a 

EL2 


/CONST 

/N.l 


(J»H)).;r.E*S)    J3    T3    9* 
8808*1 

(JH.GT.£rSI    ;0   T3   33 

8808*1 

«,J,«1I    S3    TO    37 
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IS 


29 


70 


MV171       )PT«i      PNONR  FlN    *«••%*! 

suntouriMk   itLi  CN.iTER.iTRNAj.tPSLrmi 


GAUSSIAN   ELIMINATION    »ITH    UNEAR    NATION    SYSTtN    V    <  THt    •ild-tr 

hano  sioc  vector   i  is  integrateo  [3GCT-4S*  vith  ths  system  14- 
trix  *  in  r -te  mtrix  vi.     t re*  Is  t-h.  -iuhkck  qf   HERArijNi  «e- 

TURNEO    TO    MS    CALLING    i>RJG!AM.      THE    ITERATION    TO    lN*«Jvei    t-iE 
SOLUTION    SO    «ANV    TIMES.    JtTlL     Thf    JC'HIttNG    SRRJR    IS    .(ThIn 


C  .       THE    NAGNITUOS    3f    S»SLJN.       SH3JL0   THE    'UtMaiR    ilF    ITERATIJN4    it 

10        c    .  ctto  itrnax.  r-46  su3rojtime  ; Si*«iM«.r cr  the  calculate  »no 
c    .  sets  iter  rj  4  negative  value. 


AUTHORS!  KONRAO     VQG61SJ.     HANS    ItSNCR. 

OATfci  OECENSt'R    IS.    U80 

VARIAILE  EXPLANATION  TYPe 

J**     •*•  VjSSIS'gf  "iayirnss"****""***  '■iic-?j 

ITER  NuNStR  OF  ITE^ArlJNS  rtl  8C  'STURNEO     INT£GE« 


20  t  .  >0  GAUSSJ4N  ELIMINATION  »CP.FQRHEO 


II 


CORRECTLY. 
«0  HO  CJNVeRGEiCc  3f    Hi  CDLUTION  ?0R 

iteration  jio  T£»Nt«»r:  4fT£» 

ITRMAX    ITERATION    5T£?. 
RNAX  MAXIMUM    HUNGER     3F     ITERATION     ;TCPS  INTEGER 

SLON  OESlREO    RESOLUTION    Of    T HC    SOLUTION  MEAL 


C  .       4(7,71  SYSTEM    MATRIX    A  REAL 

30  C  •      3(7)  RIGHT    HANO    SIOE    VECTOR    3  REAL 

C  I  INPORTANTl     MATRIX    A    ANO    VECTOR    3   HAVE    r0    3E    'KOVIOEJ    THROUGH       '. 

C  .  THe    LAStLEO    COMMON    AL3CK    4U3,       30TH    4    4N0    )    .ILL 

C  .  *(U    JN0ERG3    CHANGES    OURIHG    CXECUTIGN.    I.e.     Tie 

39  C  .  £LINUTIJN. 

*J5SS5S***A8i*/"v**** •••••••• ••• 

REAL    V(7.d>,G<n).Air.  91.3(41 

OATA    M    /3.1*1J<»2»*/ 

C  F'JRIUTS    lELi    START 

*3  431  fC«H4r    (//S34H      MATRIX    A    JNO    VECTOR    3    VOP.    cLintNATION/l 

602  FORNAJ    ((41,13(J«.I>  I  »» 

603  FuRNAf    (2X»iZ,13(:x,lPE9.2>f(  (R(»L3(1X»  l(»f  9.2  J )  It 

43*  FORMAT  (///2«H      'IVjfUteHTEj    SYSTEM    IN,:3,13n.     IT  jRAriONSiCHR  I T  T/ 
•  ) 

JO  305  FORNAT  (A//2X»I2,2*H    <OEFF HI  ENTEN    AUS    GAUSS/ I 

40*  FQRNAT  ((2X»:2*lx»i*E6.2l> 

407  FORMAT  (///39h      <LelNSTER    11  IL;NVS<TORtCHrtITT«INKeL    ••£12. 3> 

SOS  FCRN4T  (27X.12H      ReiI3IUN    ..-U.i.SH    -(AC.H,  13.  Urt    i 'eRAf I JmInI 

409  FORMAT  (///icx.llH      FcHLSR    V ( . 1 2. in, Tl. 3h >    «.ei2.i/l 

35  610  FQRNAT  (///2JH      !rSTSH    1ACH    E  Lin  [N*r!3M/l 

C  FJRflATS    IEL2    STO> 

L0»6 
N«N»L 
tO  *•*-! 

REARRANGE    VECTOR    3    IN   ARRAT    V.    IF    N   IS    LCSS    THAN    7. 


IF    (N.E0.7)    GO    T3    i)2 
45  90  01    I'I.n.i 

01     VII. HI-     ll.it 


FRINT    5TSTEN   NATRIX    A   ANO    RISNf    HANO    StOC    VECTOR    »    8EF0RE    THt 
ELININATION.    IF    IFRtNT    IS    SET    .TRUE.    . 


92    IF    (IFRINT)    WRIT?  (LO.  »0l » 

IF    (IRRINTI    «RIT:  (LJ»  602)     (i.J«l.N.i) 
00    03    I«l» N.l 

03    IF    (IPRINT)    .RITE  1L0.  40))     I  .  I  v  (  I  ,  J  1  .  J  •'.,  n,  I ) 

79                                          IF    ((PRINT)    .RITE  ILJ.  6i32)     I..4«L.N.1J 


LINEAR    EQUATION   STSTSN      VIN. H )««(N )•»( M#H»ll 


•o  33 IX  \:\:V.\ 

9*     4(1, J  «•»(  I.w) 

oo  os   :«i.N,i 

09    S(I)«.0 
39  i  CHECK    FOR    STABILITY. 


u: 


susftjuii.it   ul2  /wi7*     iir.j     <iov  e^i  <..o»--.« 

c 

f  i-.o 

06  sm-.o 

90  JO    is     1*1, N,l 

00   07    J*t#N»l 

*    ilil«iJ»I(uU  I  ) 

30  io  r-L»«»t 

93  JJ«l»l 

jj.i?     J'JJ.H.l 

30   It    I3«l#N»l 
99    5K*SK»¥<I#I31«V<J»I3> 
100  f-S«/  (G(!  )«1|  j  )) 

IF    <A8S<F).<.r.Aa$(FI>)    ?I*F 

io  continue 

FJ«AaS<FI) 

fNacosif  t  > 

10»  F  I«(f  l»l»0. U*\ 

I  START    ITERATUS    TO   FIND    »10    MFROVE    :OLUTION. 

30    27   <.J*l. 999.1 
110  I  T  £  a  - c  « 

00    13    J»1,K»1 

FIVOTISAllON   1-40    SAJSSIAN    STEft    FOR    COlUtlH    J. 

115  IX«J 

i-VIJ.JI 

00    12    1-j.kiI 

IF    UJSi»II,jll.,r.4iS(!ll    SO    T3    11 

GO    TO    12 

120  ii  :»■[ 

£«.<Z*4I 

12  CO,NTJ«tuE 

123  8.1  .Mi... J* 

13  ru.  i  >•«<  I) 
SCNI..0 

1»    -««j»i  > 

130  00    IS    [«1N.N,1 

n**<  I  .w >SV(J.J ) 

oo  1.3  *.n*i,n 


«lXfMI*V(!#«KW(<J*«RI«« 


13    CONTINUE 
133  t»* 

IF    <V<I*il.cQ..<<>    SO    TO    32 

SOLUTION. 

1*0  00    l»    1«1>*»1 

IS   6<I>..0 

S(N|..tN.ni/V(N,Ni 


00    19    II*l»*.i 
I.N-U 


1*3  IF    (»<I,i  >.fcO..OI    S3    TO    12 


JO   17    4  J-  1.4*1 


IF    U.cO.I)    S3    TO    Is 

ISO  17  in  i.si  :i— ;ui»vi  :.*  > 


il 


CONTINUE 

£(I)«S<Il/v(I>I) 


CORRECT    SOLUTION    ?r    »HUNT    OF    iRRQR.    .MIC*    «AS     IXV01V63    In    THIS 

oration  s" 

RATION. 

20    S<lT«6(ii*«(I> 

CONFUTE    THE   UMItnt   ERROR. 


153  C  III?Ali?N    STEP.       TAKE    MIS    SECTOR    »S    START    VECTOR    FOR    Trie    SEAT 

C 

M 

160 


00   22    I*1«N,1 
»< I »rt  l-.O 

163  90  21    J'l.N.l 

21  «T:,ni«v( :,«i»s»j i»«( i* j  t 

22  »<I»NI«AU#N)-V<I»l1> 
R*.0 

00    23    t-l.H.l 
170  23    **R*i<l«NI«v<t*N> 
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SUSRJUTtnE    1EU  '*/i>5       ."•;      P101F  rH    ».s»o.a 

IF    (R.lT.EFSLUII)    it)    T3    Zi 

IF    <ii.it .  itb'hjh    ;j   u    13 

5u   2*    l-l»n.l 

173  :*  m>«i<i> 

2  PRINT    STSTEN   *»T*IX    t    kttO   RISHT    -UNO    SIOC     /tCTlU    3    XFTSR    rue 

c  hvotisation. 

ISO  IF    CIPRINT)    WRITE     (U.     6J<.I    LX 

IF    flP4INT>    »i»IIc     it  J.    402)     <J»J«l»N»ll 

29  if   ii'ttNT)  .ritc   ai.  603)   ;,i»i  r.Ji.  J«i»«#u#sm 

IF    ll^lKfl     aitlT:     iLu»     IJil     IJiJ'l.flilJ 


119  00    26    l«l»N#l 

26 


i« 


oo  2 a  •<«U'*'  + 

STOF    LOOF     10   FINO    53LUT13N. 
27   CONTINUE 

REARR1NCE    SOLUTION    IN    TH(    tRRXf    1, 


It9  28    CONTINUE 

00    29    l«ltN,l 


200 


29    m»4i«i(  I) 

PRINT    SOLUTION,    IF    IPRINT    IS    SET    .TRUE. 


JO 


209  C 

e 


IF  (IPRINT)  WRITE  (L  )>  60>>  4 

if  Iipkint)  Kir!  uj>  o06)  i  t.ii  ti»  t»:,.'i»D 

IF  «I»*INT>  wRITe  U3»  607)  f  I 

IF  ItPRINTI  -a;ic  llli  6041  4.LX 


2X0  00   31     I«i»t»l 

-   (If-- 


ORUCAcN   0(R    'illlJ   4    JNO   3ES    VECTORS    3    HACH   =  LIP  I  N»  t  IZ*. 

IF    IIM[*n    WRITS  ILO.  610) 

IF    (IPRINT)    MRITe  ItJ.  602)     (J.j-l.fi»i) 
00   31     I«i»t»l 

31    IF    (IPRINT)    «RITj  (LO.  603)     t .  <  »  <  I  .  J  I  .  .' -:,  ■< .  .  I 

if    (IPRINT)    .RIT3  UJ»  602)     (J.J«l.N.;i 
RETURN 

219  £  i«*OR    ENCOUNTERED.       STOF. 

J2    ITER  — I 

«»ITi     (L0»    6C4)     :.I.<lI,i) 
G(l)«S(I)/V(If I) 

220  RETURN 


33    ITER  — ITER 

229  JifuRN 


NO   CONVERSENCE    OF    Tit    SOLUTION. 

ITER 

SiTU 
ENO 
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FUMCrtl*   fi  7*/l7>       JFT«2      P101F  Fin    ».J»«.»-t  c(.y4j 

1  liU    PUNCT'.OS    F2    (CJEF,XI»H.J*0»1) 

2  1     cohfute  function  vuj;  r.«<«ii   of  »  cs-'oltnoiuii  using 

1  ',       AUTHCMl  -(INS    IE3M6*  '. 

10  C  !       <i*I48Lc         £XFL4n4riON  TfFE  ! 

2  :    tW»r      »a£?S8nl«'!3i»!2iEfi« Siitt***  : 

C  .  XI  Xl-*»tUfc»    #16*6    f    «Ul    JC    :n.1»uTC0.  *e»C 

C  .  H0B0«1  WDER    If    TH6    'OLTNO^UL    IN    TM£  IlTeSe* 

13  £  .  <l    JiOcCTIJn. 

C  .  FZ  JftOlN*rE    If    TIE    !3-"lLM0i:il.  >EA(. 

«Eli*caeF(7) * * *** ••••• 

ZO  *C<HFi»>iO*0»l*t 

T-c3Ef<NCa6Fl> 

ao  ci    i:-i,NO»o«i,i 
NncoEFl-ii 

31    r«r*Kl»CJ£r(I) 
25  F2«» 

IETMN 
END 
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sioc*  3*r*  H» 


to/in     itt'i     MO** 


fTi    +.a»**i 


5t/»i/'7.       4-.. 


'Axe 


SLOCK    3»U    »S3 
COfinON    /    4S1        /    «,l 

[MO 
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HOCK    OAT*    ;*fij 


/WWJ      Ji>W      »i(Hi« 


Flrt     *.')»-»-»0 


SLOCK    3*r*    :»r»3 

COkhon   /   jir»3   /   <l»'2»r 

*€iL     Utlo.ioi.iilio.  wl.mo.lot 
EM 
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PHUGfun  aeioi  7*/i7s      :*<•>     ^ion»  ftn   <.o>».j 

1  PROGRAM.  3E«03    i  isi"jt.  Jur»ijr,  r»»f  i,  rii>r'6»r»PEe»  f»»t -«i  ~<-'jr  i 

'.     this  is  4  oshxistratio*  program  *mo   :t  sho*s   r««  ca*»£cr  use 

i  C  .      Of    THE     TPL    JUART    H3RARY     jSHw    curPUS. 


10 


25 


10 


»J 


50 


70 


AUTHOR  I  lA.lj    :mner 

9aT£i  JSCeflStR     H»    1980 


connoN  /  afid     i  mr« 

COMMON    /    OAlAJ    /     <1»<J»T 


19  'SAL    P10TRCZS6I 

REAL.    Illlft«loi»X2(16»16l»Ml6>16l 

real   :atf  t*«i.  3ji«iir.7> 

INTEGER    I»AW<;> 

20  LOGICAL    I  PLOT 


OATA    ICNT       /10-t»  / 

OATA    IPLOT    /.TRUE./ 
OATA    PI  /i.l<.159Z6S*/ 

FORMATS    OEN03    START 
101    FORMAT    I///72H      HeLLJ-      THIS    IS    P«06RA.1   PE«03    ANO    «t»LL    i?«»    -tO*    T 
•0  USE    TmE    nARV€llJUS/73H      TPL    UNARY    l:3»aRY.       rOu*    INfE-lEST    .ill 


»u  Uje    thi    nAPVSLL JUi/ 'an      i'\.    jinai«t    L.a»»KT.       rjux    MrsntST    i 

•  8£    SREATIT    A»A*OcJ    IT    l4I;;*/»h       AMO    fASTCR    PPwiRAM*I)u«///> 
F3RNA?    (22H      PLOTTER     [N ITI AL I l\0.  i 

ijo    FORNAT     UlH       3AIA     >ll-iTS    RE43    ?R]N    T4»t     3.) 
0*    FOKNAT     (i<»«,5H«l.irs»'j«.il«ltAX»ftA.  JHX2i:'!»0Af  Ji^ZfAI.  U  .  .H  r  «  I  *  ,  7  i  . 

•  <.nri»i/Si.  7NRL  TTT  E^.KH.r  10.  3»/-»X#*hijc;?,^(1  J,  fic.i)  I 

SOS    FORMAT    (*7M      PLOTTER    1N0    JSE»    tRES    OErtnCO    ANO    AXES    3RA*N. ) 
34  fgR«A£   Men     hIasurEi)  ->ata  sjrfacc  jpa^m.) 
07    FORMAT    169H      ;hT?5    }R0£»    JF    'CB.YN3NIAL     1 71    U-    ANO    « i— JiReCT  ION.     IP 


•«[tr     ANO     ISTAT) 
109    FORMAT    \l<,*      idRFACE    FIT    :3l»l£TE0.l 
!J5    PuRiAT    IZ9H      APPR3Xl.1ATsO    iJRfkCe    Oft/.«M.) 
601    FORMAT    (1H1/1HT) 

101      POR»Af      !11.»»U.»JI 

»52    FORMAT    liiol 
S03    FORMAT    13HS.0) 

FORMATS    JEiQJ    STOP 


t/0    ASSIGNMENTS. 


l5«« 


■LOGLUtlSESSIl 
•6 

LP«13 

«»IT£     (Li.     Ml) 


INITIALIZE    THE    PUTTER. 


65  ^         !l*!!Kfif!'SX£C*!AtK'{CH 


MCITE    (11 


)    CALL 

f  ;o2> 


VARIATION  OF  PQLTNOM.IAL  3EGREE.   L3HP  STMTS  HERE. 


75  mni:i 

HajTRT.l 
MZSTOP-6 

00    04    M0R9R1»N1JTRT,M1ST0P,1 
to  30   04    NQRDRZ'NZSTRT.NZSTOP.i 


as 
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»«dG»»*    3;1Ji  7WL7S        1PW       »«0«P  FT*    H.3-,,1 

c  .     jeao  *£»su«6o  o»r»  r^ii.*  tape  <  »hp  ::nucriN60dSLr  se*i»c*  r*t 

C  •       «|N{flU«    4NI)    1AXI1U.1     VALJE     3F     S4CH    <1*     (>    »NO     T     (  »e  J  J  I  •<  ^  3    TJ 

c  .     set  rxc  'ujt  at^i  ihc  umonis:  If  a  s>*e;i  :r  o»»c-«i. 

M  J  ! 

JErfiNO    LO 
11"!»    •     *9<)S<)9. 
Hfil     •    -9Q944, 
X2pUN    •    •«|)qgo1 
«9  I2IUX    •    -'»<»<»'»'». 

TrllN  •  .«09«J). 
TBAX  •  -9<»999, 
IVAL-IOH 

oi  »eao     i lo*    son   ival.irah 

100  IF    UVAL.eO.ICiri     SI    M    ;l 

ogcoot   i30.jC2»;uii    i*"<r]i,i»*rsj 

00   03    IW.NPNTSLi  a 

oo  0)   i2»i*r)p*rs2*i 

I¥4L«lpM 

109  02    »6A0       ILO»    Sul)    I»AL»I*A4 

IF    I  IV4L.E0.K  in     iJ    n    32 

3EC03  6    (<.5»oC3»I*A.<)     t 1 1 II  .  J  2  I.  ttill,  {  2),  »  {  U.  12  J 

IF    (Xltil»I2).LT.xirlt1l     UMti'UU  1*12) 

if  (x2tu*i2i.Lr.x2*rii   x2*(.i><2(i:*:2) 

110  IF  (T(i;*(|)   .Li.r.ii*   i  tm*  •ttiipU) 

IF     (Xl<II»l2>.  if.XHAX)      <l14X.«l(Il*I2) 

i!  i*fiil«l«»-<:i»x2^**»   «2*4x-«2(ii,i2) 

IF   (TTll.IiT   .ar.rixx   i    ri4X    •rtu.12) 
03    C3KUNU6 

119  xairs  ili*  loi) 


120 


OtFtftt    »L3TT6»    ANO   USE*    AREAS:       0P.4U    A.XES. 


55 •"IW 

129  TO  •     3.0 


liPfUH  •  .0 

a»«»»  *  1C-.0 

(2F1IN  •  .C 

I2PAAX  •  10.0 
130                                       tphih 


TPBAX  •    10 

UU»IN 
I1JH4X 


:8 
■1:8 


I2Ur)IK  •       -.1 

139  X2UMAA  •         .0 

TUNIM  «-3i.a 

rui»»  •    35.0 

4LF4X1  •    30.0 

4LF4X2  «i:u.O 

1*0  XLFxr  •     H  .0 

PL3T*(6l)  •    (4LF4II*»(M190. 

'LCfJCjil  •    <ALFAX2»»[)UB0. 

PtaTRI63>  •    (4LF4T    ••II/UO. 

»iaT*i»*)   •  xa 

1*9  »LQTai«5)     •    TQ 

•LOT»(0»I     •     (XiP*AX  -«t»*M  )/  (XlUr)AX-«lU"INI 

PLBTR(»7)     •    <xI>"UN»«lU*AX-«l»*AX»XiU.'1IN>/<  UlHAX-«lJ1t"l 

►LStaijai  •  (xiHix  -« I ?<•  t <•  j  MxIuiAx-«iii*f  *) 

FLOTsIb'JI  •  <  iJl>*lN««2'j'1Al-»2PlAX»X2lH:n)/(x2unAX..«iiMl"4l 

190                                       f>L0TRI7C)  •  (    fP^AX  -    r»»t*                    /(    T'j*4X-f up  I*    I 

?IQTR(71>  •  !     "»*r<«rttr14X  -    rP*AX»TUr!lM  )/(     ru<i4x-ruM"«     I 

XIL  •  (X1P1AA  -«lP»t,«  » 

X2L  •  (X2P1AX  -«2P»l*  I 

rt  ■  I    rpixx  -  r?*M  ) 

199  1RITE     (If*     10*)     IlP'»IN»Xi»14X»J2P««IN»X2P,14X»fFr)IN.r»»4X.Xll)«!N.XlU 

•  hax»xzu«I*»x2U.iax»  rgii*.  rj<ux 

IF    (IPL3T)    C4LL    AXe*     (Ad.  TO.  X  iL#ALFAXt,CHXl,»2,  XlU'UN.Xl  J*4X.*i4FT. 
•1»7»4) 

IF    (IPLOT)    CALL    4XES    HO.  r3.X2L* ALFAX2*  2rtX2**2»  X2Urf£N,  X2"Jp|AX.-.riF7. 

''    UPLOT)    CALL   AXES    UO»rO»fL    »ALFAT    .CM*    #-2#rUrlIN    ,TUr»»X    »»HF7. 


IF 
•3*7,4) 
•«ITE    (LI*    103) 


169 


ITO 


3RAM    iEASURENEHT    3*TA    SURFACE. 


149 


HJiKkn    3£iu3  7W17S       ,>:•?       ■"•Ofi»  PIN    1.1x11 

oo*5J*I  w#npmHI#I 

IE  ({ploti  call  rnnrw  uPi3T»»PL3r»x<  (i;»i2)»x;<  n»;2i.n  ii.iji  ) 
its  if  (Iplot   .Ana.  iz.a.ii  :»ll  'lot  (*?l>,tpioJ.  ii 

0«    IF    (IPLOT     .AhC.    I2.iT.ll    CALL    PLOT    ( XPL  IT,  TPL  JT»  J ) 

80    0?     [2M.NPKTS2.1 
0  0}    !>[,f<PNrsi,  i 
IF    (1PLQT)    CALL    MRT*    (XFL3T.r»tOT,XM!;tI2).x;UJ»I2).r(U.I2>  > 

no  if  iipl3t  .ano.  u.£a.n  call  plot  <xpl*t,  tpljt.  n 

99  if  uplot   .ano.  li.sr.ii  call  plot  (xPLQT.rpt.or.2i 

■MITE    (LI.    106) 


114 


1W 


CALCULATE    SURFACE    FIT    M43U6H    DATA    POINTS. 


9*    .UITJ     <LI»    1C7) 

PEAO       (    9,         •)    *OR0Rl»*<ORDR2.tP«tNT,ISTAT 

ng*ii(«i««g«oii 


N3K0a2«ftQR0ei 
193  IPRIHT-O 

ISTAT«2 


l«*o«ori»i 

*«H0»0R2»I 

IF    (L.LT.l    .OR.    L.C-T.7)    S3    T3    36 


200  IF    (H.LT.1    .04.    l.iT.71    S3    "3    36 

IF    (L'H.cT.l    ,3R.    LM.Sr.t9l     S3    T3    04 

CALL    "AT3     <r»PNTSi#NP*rS2,L,.1.C3cF.IPPIMT,ISTAT»Sue,U»U*«ORI 

if   (ierr3R.lt. at   ;a  TO  39 

MUTE    (Li.    IC3> 


COS 


210  C  .      ORAM    APRROXIHATEO    SURFACE. 


OJ   07    I1»1.»PNTS1,1 
30  u7    J2»I»nPHT52#l 
21S  TAPP-FJ(COEF.xllIij  I2)t«2(  II.  121  .N0Pn«l.«flP0«2) 

if  iiPLOri  call  ri»Tw  (xPLar.TPtor.x' i;i.i;),«'<ri,i2>»»Ap»i 

*F    (IPLuT     .i«3.    t2.e).l)    CALL    'tJT    (  ?»C*S  mot,  31 
07    IF    (IPLOT    .ANO.    I2...T.1)    CALL    'LOT    (XPLU7,rPL0T,2) 
00   US    I2«1»nPnTS2>1 
220  00    CS    Il*l.MPNTSl.l 

TAPK««HC3£F.111I1.  12).  X2li;»I2)»MOP04;<MnP0R2) 

if  (IpljI)  call  mSN  («>Lar,rFL3r.xI(!*,  :2),xJiii.i2)»rAP»i 

IF    (IPLOT     .anO.  U.E0.1I    CALL    >LOr    ( X  PLOT.  TPLOf.  i  ) 

00    IF    (IPLOT    .anO.  U.iT.i)    CALL    »L3T    (XPL  JT.  rPLOT.2 ) 
22?                                    ■RITE    iLi,    :o! 

09    IF    (IPLOT)    CALL  SUPAS6 


230 


.       TERHtNATE    GRAFICS. 
233  If'uPLOT )*CALL*STOF«* 


STOP    0077 
2*0  ENO 
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jwi7s     ji»r«2     p.iqnp 


c(  Ui.1  : 


1 
I 

10 

1J 

20 
29 
30 
39 


A3 


99 


60 


49 


70 


T9 


S) 


ij»nuuriH6  mt3  ( npnts i. ** its z.lpn.coef,: print.  isr»r»sjdii.i:»«j'(i 


preset  systeh  i»r<tx 

3U->>0krN0Al4k     4PPR0XII1Ari3N. 
iS    SQXVEO    it    USIMS    A 


4N9    *ICHI    HAMP    SIDE    YfcCTUR    1    tn    4 
TI3N.        THE    t:tlC»«!    d3U4U3N    .TSTeN    Trtfn 

*uss  jjao»N  ec:*iN*riaN. 


^NS    ZS.3XER 

3ECErt8S«    U,    1)80 


UJlHOkl 
OATEl 


'4RI4SLS    :>Pl*f4*rtON 


TYPE 


npntsi 
npnts! 

:jep»« 

k 


nunijer  of  3*r*  point 


9i 


NUN^ER    OF    OAfA    POINT. 
»l)lTNlNjAl    ^EFFICIENT 


IJI 


CJnSTanT  <i 


/■.r  cinstant 

NUNSEH    OP    Cliff  tCIENTS    FOP.    THE    POLY 


IPRINT 


[STAT 


NQNUk  IN  XI  OtRECTION, 
NUN9E*  IP  C3EFF  ICIENTS 
NONIAL  IN  <i  3l»eCTIQM, 
CONTROL    OUANTtTr    3F    PR 

•i  print  equation  syst 


P3P  The  03LT- 


AN3  APTER  T-tE  ;*USS 
"  INlNAriON. 

IN?  SJUATIO'i  STST 

'"rft 


,MT  OUT 


i  aft:» 

.nii.cTiati. 

MS  WXT.-i 


n.*l 


SUSHI! 

tERiat 


an*,  lii 

Jif 16.16* 
Tl  15.1b! 

IHPORTANTI 


«0 

control's  tat  is  tsli5"ra 
dir»[/<;  print   oroinat 

/:STAT/«l    3NLY    PRINT    £f 
ISTAT>3         S<rT    4$S0L'JTE 
ISTAMO         1ST    ^eLATtvE 

;u4i»r»r»    i.   class   in 

■utrix. 

;r*3r  fl4s  returnso   m  the  callin6 

PROSRAN. 

••1    NO    ;RR3R    OETECTEO,       'VJ10    JJB. 

—I  an  ?rr3r  occureo  !'i  thu  GAu»Jr»i« 


iZ    P'tNTEO 

»t>::  of  oata 
iPRipr. 

EPP.T. 
:P.r.tl». 

rnc  ;tstep 


INT^ieR 

IfUiitR 

INTeGcR 
INTEGeR 
I.NT6GE* 


htegsr 


<EAl 

INT-5-R 


ELi*tNAru««r  :3ep 
accurate,     jevj 


FICICflTS    Ai<£    NOT 

papai:t£»s. 


xi-cooroinates  op-  rsf  luri  »oimts 
xz-CJanoiNAres   ■&   rs£  ">a;\\  points 

T-C30R0rUrF$    OF    THE    OATA   »3INTS 
SYSTjl    "4TRIX. 

Rii*f  -*4no  izac  sector. 


<£4k 

«^4L 
ReAL 
•tAk 
<eAk 


COnnaM    /    9AT4J    /    U,tz,T 
CCnnON    /    433         /    4j  1 

>E«k  imi«m.i.»<i 

«£al   :3EP(«9i.:u»«i(i.^i.ST4r<is>i6) 


101 

401 


»oz 

I 

603 

»0* 
609 
466 
607 
609 
609 

til 


S0Hn4T5 
=  :imi 

'OF.  ci'j. 

f  a««4T 

>91>I3/4 
'//♦OH 

»P    3474 

pa»fl4f 

■tI«t'*4T 

pQ*«4r 

soar* r 
P5»«4T 
Fu»n4T 

»o«n*r 

FC«N4T 

Pom»»r 

FOR*  AT 

Fu«n4r 


»4T3    STAB 
i  /  /  /  7H      AF 

J»ii«  t  I 
I///51H      3 

1H       3«0FJ« 
N'J«8E«    :p 

points  jit. 

!///62rt      J 

IOM/) 

l(«((i3C8X 

(///A«H       ; 

(///ZArt      : 

Ilii.iMJ    • 

I2l>3nl  •• 
(  /  /  /  3  9H  3 
(///A4H  J 
(///30rt      4 


rE».J3»»l*    STEPS    NO    SOLUTiaN    WITHIN    t-IE    4CCU84CT 

»Q£«     IP     »P»»0II«4T:T4     'nLTNa«t4k     IN    «l     JfECTia^f 
JP     iPPROt  M4ri3N     '^UTMIItAL     IN     «2     JUSCI  I  3"<»  S  «  .  1  5 
34T4    >3INTS    <ITh    ca'IiTA:iT    H*l64»I3/'.uH      Nu/liE*    3 
X    12NST4<r    ti,  lit ,  11/  ) 
rSf«N    >i;JU    4    .-10    /£CT"!R    3    3EP3RE    24USS    J3R0AN    £ 

!t..p?i>.zt.ii/((*»»i3(:t.iPE<».ziin 

3U»fl5N    SrjTEl    4PTE»    OAUSS    J3H0AN    ELlniNATIOM/) 
JeFFICIENTS    C3EF(I»J)/1 

i.:..:ni.^.3i.ii'3(»*i.:iu>r.».i]/ii 

J3IN4TES  OP  34T4  'OIN''-:   -   AS  ilVEN/) 


,  'Of"IT. 
iSikUTE  s*«3«  4T  E4CH  POINT/ ) 


»S  »PP*0«I^4T60/1 
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SU««JUflN{     nl)  ;»/173        JfT'i        .-HO''  *JN    ..,,♦-<■>  o'...,'-^>.      t-. 

»12  format  (///354    ^tL*rwt  ch»o»  »r  cach  »-H"it  in  s/i 

»1J  FORMAT  <///2in   ittllSriCS  i  JM-M8YJ  //  ?»  ,1  in»<  t«  A»«  r««iJ9  .  I  3«.  1  •"■*. 

•  2/7«.i»HM»iifiyi  e»«n4   »t  :•.  1 3.  ii.mji,  :%  i».  irn.i  > 
814  Fj»N4r   i///jih     statistics   su«-»»«t i//7r,.j.-iAVLi'.»iii   t*«j»»is«.iPP9. 

90  •Z.iM    t/7J,19MN4Xl,1UM    6»9Q»    »r    f«»I3<:Xi»'U«#tj»i»ii?F9.i»iH    CI 

415    FURNAT    (///7M      »FTE*.  I  J,»IH    STEPS    "H   iU'lTiaN    4U*IN    T-tr    ACCURACY 


100 


x*o 


►  GBmai!    1*?3    Stop 


99  LI*l.aCLU(ISESSil 

ItRROR-l 

[TRMAOiO 

EP$L0N».lCi0E-<iO7 


i>«[ht  ^nci[»*r:w  parameter   if  :f>«r<r  is  i»t  ate  *  man  «. 


MCRORl't-l 

NOR0R2-N-1 
109  IF    (IPRINT. LT.OJ    03    ")   31 

«»t  t£     (to.     oCi)     ■<<3.»0»l»,«0R3R2t<<PNTSl,rl»»TS2 


c 
110  01  ipowri«o 


START    TO    catse  T   E36E    SECTION    SUSH A TRI  =  tlS    1.    CLASS. 


,    U»    ««l*Li*l 
If    (K.4T.L1    SO    T3   02 

119  TOP    cDir    SECTION    SUJMAT9ICES    1.    CLASS. 


iifss:;1 


jsTA*r»<K-ii««*i 
120  jsrgp-«»n 

so  Ta  oi 


l 


129  02    ISTART»U-L)»M*l 

IST0P«<  tK-U«U«1 


DISmT    h»«0    EOSE    SECTION    SU8NATRICES    «.    CLASS. 


i|f  ART«<L-1I*N«1 
JST"" 


iI0»«L»1 

130  C,     PRESET  EOSE  SeCTION  SUSNATRICES  1.  CLASS  IN  SYSTEM  MATRIX  USING  1Al3i. 

03  CALL  HAT31  (NPNTH*NPNTS2»N.IPaB«l»SU3M:j 


30'° 


_    (.4    I«lSTARTiIST3».l 

139  II- t 1*1 

00*34  3-iSiA«.r,jsrap.i 

34     »(I, J  l-SU8"lltl»Jl) 

140  09    IPQWR1*|PQ»R1»L 


COPT  3EPINE0  SUSNATRICES  I.  CLASS  3IC0MALL*. 


149  30   )>1    Il«l*Ll»l 

30    37    Jl.i.t.,1 
IJTA«T-<Il-i)»N»l 
!ST0P«I1*N 
jftART»iJl-l)«N»l 

1J0  JST0P.J1»S 


139  C 


30   3o    l«ISTART,IST3P.L 

30   06    i't  J  IA» r.jSfuP. i 

06  »<  !♦*.  J-1  )•»(  In) 

07  CONTINUE 

PRESET    RICMT   HANO    SIDE    VECTOR    I. 


t-0 

30   08    LI«1,L>1 

ip3*ri«li-i 

oa  os  mi»i»n»i  .v 

[pguR2«ni-l 

OS    8<I)-S3INPNTSl,NPNTS2,IP3*Rl#lP0\(a2#ll 

169  Nl.L'N 

PRINT  SYSTEM  MATRIX  AND  RISNT  4AN0  SIDE  VECTOR  S»  IP  IPRINT 
{      IS  SET  TO  2* 

170  If  I  IPRINT. NE. 21  SO  TO  10 
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susroumne  lirj  /wws    w\*i     itioiw  mi 

WRITE  (10,    602) 

WRITE  (tit    ->uj>    (  J.j-  i,  Nl,  II 

30    59  I'uiH,! 

99  write  uo»   sci   t,  ui :.  j».j«i.«i.ii.&<  i) 

179  **lTi  '(.J.    6031    U,.-uil.tl 

ifrint-i 


210 


220 
229 


240 


2*9 


sauss  jo*04N  elimination. 

fsljn«efsl0nw"u 

ul  isl3  (nl,  iter»itr«u»..erslon) 


HO  10    ERSLUN-EFSLONWU 

CALL    I:l3     INI,  IT. 
IF    UTER.LJ.U)    5J    TO    27 


109  C 


c 

c  frint  euurtion  ststen,   if  imiht  is  set  to  1. 


IF    UFRINT.NS.l)    JO   to   12 

VRITE    ilJ.    005> 

WRITE    UJ.    »tn    U.J«>l»Nl,ll 

"ITS    ?L0»    a 


190  11  ««:fc   (L3»  ».•>>  !.<*( :. Ji.jti.u. i)»4i  11 

i«IT6     UJ,    603)    (  J.J'ltitt  II 

c         ocFiiik  coefficients  *mo  »»f*r  rwn*  :f  iprint  is  sreats*  man  3. 

1*5  12    CONTINUE 

30   13    L1«1»N1,1 

13    C3iFUU»«iL;i 

Ifiirrint.lt.o)  so  ro  is 

WHITE:     (LJ.     3C6I 
200  WRITE    ICQ.    ai}7)    (  12.  Ii.l.1.11 

30    !<•    il'l.L.l 

ic$T«T«m-n»^»i 

ICSTJFMlwN 

i»  write  <lo>  »g«)  ri#<c3EF(ic>»i:»ccsrnr#icsTnR>;) 

20!  WRITE  ILJ.  607)  <I2#IZ«I»1.1) 


STATISTICS.  IF  ISTAT  ANVTNINS  ELSE  9UT  ... 


15    I! 


I  ihsfif:L?:3!  ihlUt 

IF    (ISTaT.ST.O)    IERR0R«2 

istat^iassustam 


4  confute   absolute  error. 

219  C  t    •(ARFROIINATjO    VALUE      -      Jt V6N    VALUE)     1 


■  00   1ft    tl«l,MPNTSl»l 
30    Id    I2«l>N?«rS{,l 
16    lTAfTll.l2)«F3(C^eF.Xlt:i.  12), (2<tl>:2)*'IOF3Rl«NOR0eZ)-r<  11,12) 


i  CGNRUTF    JELAIIVE    ERROR,     TF    IJTAT    WAS    L  C  TS    THAU    I;»9.    I.e.     IcRRUa-l. 

C  I    •♦AFPROaINaTSO    VALJE      -       i(vE1    V  ALUfel  /(GIVEN    VALUE)     « 

IF    I  I  ERROR. ME.  1)    50    T3    IS 
So   17   :i«l,NRtT';,i 
30    1/    I2«1.mP*«T$2»1 

17  SUT(li»I2)*STAfTlI,I2l/T(tl,  12) 

J    .  ERROR  STATISTtCSi   CJNPUTE  A'/ERAiE  ERROR  4N0  FlNO  LOCATION 
230  C     OF  •Ailnue  cRRlR. 

C 

18  3RR1VE-.J 

iunAi**isisTArii»l)i 

239  JUKI 

00    19    Il«l,NRNTSi,l 


30    19    i*>l,NPNTS2»l 

If  ,*aicfiti  :i.::n  .le.errnax) 

;nMiti:jn;i,i2) 

INAA-ll 

j«r»  x—  r  2 

19    cRRAVS-SRR AVE «4*S< STAT  ( 11,12) ) 
ERRAV  =  «fc»RAVE/(Ni"<rSl»'<»-»rS2) 
IF    USTAT.f.E.21    50    f/J    22 


;o   TO  19 


RRINT    3R0INATES   3F    3ATA    'OINTS    ANO    *FRROT I  NAT £0    VALUES* 
IF    /[STAT/    IS    S8UAL    2. 

".it"!  m:  *?! 


290  WRITS  <LJ.  eC3>  <I2.I2«"1»NPNTS2»1> 

00  20  lI«l»NFNtSi,l 


20  WRITS  U0,  SO*)  Ii,tr(U,I2)t  I2»l,NPNTS£»l) 
WRITS  (LJ,  603)  (  12,  I2«l»w»NTS2.l) 

299  WRITE  \ia'    603)  [ I2# 12«1»NFNTS2»1I 
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SUBROUTINE    1AT3  7W175       1PT«2      MO.IP  FT(»    i.sxo  --u/ii/i-. 

90   21    U-l.NPNTSl.l 

•l*ftPNTS£* I) 

"<U5     (L3»    5031    (I2»l2«l.*PNTSZ.l) 

C  »«I«r    ERRORS*    IF    /ISTaT/    WAS    I    1R    2. 

22  if  (tsTir.ni.it  sa  n  26 
263  if  (iE«no*.«a.i.)   ij  n  :* 


MINT    ARSOlUTc    if.»a*Z. 
WRITS     (tO.    611) 


270  write   iia.  t>03i   (:2»i2«i»(PHr$2,ii 

23    «*l"£    (lJ.    60*)    ri*(STAT<ll«I2>»l2«L«NPnTSC»l) 
WHITE    <LJ»    5C3I    CMM.iKi'iMI 
£0   TO   26 


2TS 


PRINT  RELATIVE  ERRORS. 

M  "HI!  118;  tti! 


•  RITE  Uj.    slj)    (I2»I2«lt'«PNTS2>l> 
2(0                                         00   29  U»l»N.-NTSlfi 

2*    *!fli  !Ml»    *0,,l    IU(ST*ntl>i2).l2«l»N»HTS2#l> 

•  RITE  UJ.    6C3)    (I2»!2»l.  «J-«r52,.l 


2.5  ! 


PRINT    ERROR    STATISTICS. 

26    IP    (IE8POR.E0.2)    WRITE    in.     sill     E  »»A/E  ,  !  1iJ,  J  fli«  .  -.  •  J"u 
"ERROR.  Ea.l>    WRITE    <U»    ii*l    E»*A*E/1 1AX»  arut.e  R»«AX 
2» 


«a   ii 


2M  C  FATAL    ERROR    OETEZTSS.      3UTPUT    ERROR   KCSA6C    »l»    'RESET    C3eFFlCIe*TS. 

27    ITER  —  ITER 

lERROt'-I ERROR 

2«  JA^.'IiriJNi.i 

23  auci'.o 


300 


WRITE  UU  1C1>  [TER.;»SUN 
WRITE  (I.9*  o!3>  ITER»E.»SL3* 
S3  TO  12 

RETURN. 

29  RETURN 
ENO 
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SUSUCSUTIftE    1AT31  7*/ 17S       1?T.£       "10H»  Fl 

I  sumoumke  "Am   in**isi,**<M%i,*.itn*9;, sosii) 


{  .  'BESET    5U8I1ATBI*    1.    CLASS    S-H-U    !H    thl    jrSTCH    <AT9l»    4    F<3» 

C  I  AUT*0»l  -UNJ    £E4*C« 

c  .  oatei  3tcEnaii  3i»  i?eo 

10  C  *  ¥4«IA8LE         EXPLANATION  TYM 

c  I  sJSHi"      nu*«5*tf*}It  a*  »5iif!*4r*cVisTA*f**<i      !nH£|3 


t  .  *n*isi  iu(ne»  if  3>r»  points  at  ctistant  tl  iutz 

t  .  n  NUNB6»    JF    C3EFF  ICleNTS    £O0    ME    ^QLr-  !,ilc<c« 

13  5  .  N0/1IAL    IN    «>    OlRECTtQM. 

J  .  :fow»i         power  for  me  a  value!  in  the  inteser 

C  •  "ruiiATION. 

c  .  sumkcn)  susiatrix  i.  class   in  me  jtitsh  REAL 


20 


SEAL    SUSHKlt/1) 

3U   03    *«l<fll*l 
23  IF    iK.il. IM    3Q    '1    01 


10 


TV    EOKE    SECTION    SLE1ENTS. 

1"! 

SO   TO  02 

RI6HT    NANO    EDGE    SECTION   ELEMENTS. 


13  01    IHMUl 

C  FRESET    E06E    SECTION    ELENEiTS    IN    SUBNATriX    1.    CLASS    USIN6    S3. 

♦  0  '      32    SuaitUI«JI«S3<  <PNTSl»N»NTS2iIF3Ml»l»3¥R2»0» 

33    IF0hk2<[F0m»2»1 

I  COPY    3EFINE0    SUNCNTS   OlAJONALLf. 

♦  3  Hlaft-1 

30   B*.    t«l»Nl#l 

00      -'        .•i,*.L 


RETURN 
30  ENO 


3*    Su§niU»lf  J-ll-S  J4«lCl.Jt 
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1 

5 

ie 

13 
20 
2J 
30 
3J 
»0 
»S 
JC 
59 


Si 


74/173       JPT.2      PiOIP  FT*    ».s»»i 

REAL    FUNCTION    51    ( titlTZ  I,  *PMT  S2t  r  POVP  i,  IP1WP<>  [  n 


perform  su.i"»ri'jr«s  ro  »«cStr  rut  clccnts  ki.ji    tt  Me 

JYSTEN    NATRU    A    itO    r-te    :Lt«e.-*ti    3(:)     IM    7H-;    Rt„rff     din 
SIOE    VECTOR    J    FOR    ME    lIleAR    SJUAUl.'l    Sr5T  =  i»    [-.»f    ig>T 

at  solved  >-a«  a  »airNjNUL  3o-APPfj».:iATiaN. 


AUTHOR  I 
36TEI 


tA*S     JtiNER 

oece-her  )i»   i^eo 


VARlASLfe         EXPLANATION 

NPNISl 

IPOrfRl 


A.T     CONSTANT    XI 
STANT    <2 

rue 


IP  OUR* 
IT 


S3 

11(16*161 
X2(l6ii6l 
T(16»161 


IIU1SER    3F    3ATA    JOINTS    /.T    C1N5 

iuii?»   3F  2»ta  poihts  jr  criNS 

»Q*eR    SOR    The    U    VALUE!    IM    Th 

:ui*ation. 

»a.eR  fir  the  tz  /ALut:  :m  the 

suination. 

cgirRQL  »a»»*ete* 

•3     r  Values   arS  sot  ikcluoso  in  me 

SUNIATION. 

•;     r  valued  are  NeiuOeO  In  me 
sun-urn*. 

l-»    RISMT    -(ASO    3I0C    /CCT3R    ELE*.I 
RESULT    3F    SJNNATIO*. 

Xi-CJOROISATE*    OF    THE    1ATA    P1INTS 

<i-c:mii*ArES  if  r*t  -urA  points 

T-CJOkOMATES    3F    7HE    DiTA    i"11MfS 


Intuit* 

Lsre^ea 

IiTtGER 
!*TE5ER 


ISMljR 


RS1L 
RcAi 

«c»L 


MPORTAhTl     rwj    J»TA    'TINTS    (1,     «2    A.'!0    T    Htv£    TQ    9e    »<«)V13kO 
THR3Ui-t    TH-    LA4ELE0    CJIidN    3L"jCK    3AU3. 


s« 


NQN    /     5»r»3     /     <1.X2,T 

L    Ull6»i6)»<2t.s»l3l.fU»,l>l 


00   C»    *l«l»NP»TSl,l 

oo  as  <2«i.h»nts2»i 

If  UFO»«i.ta.')»   10  ra  ?i 

|i»xIigI**ii»«tf»awRi 

51-1. 

IF    (  IP0.R2.tJ.il    JO    '1    01 

S2«»2  i*i,i;i««[?:."2 


»j  ra  o« 

33    S2-1. 

36    ST.TII«l,*2) 

.   if  ur.es. 3>  ST»i. 

93    S«S»Sl«S2»SY 
S3. 5 

RETURN 
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1 

s 

10 
19 
20 
29 


7A/17S       0RT«2       010HR 
INTEGER    FUKCIIIN    IEL13    <N> 


FTN    *.«♦*■»« 


ec/. <:/?«. 


SauSS    JORDAN   cLININATTd*    dITH    HATRIX    I    4N0    VECTOR    9.       THE 

"    r-E    SLIrtlNATION    IS    CO"R(.sTEO. 


<»aiist.£  ulj   ;:  set  ra  i# 

leLi    LESS    THAN    J    INOICATES 


:? 


•RRCR. 


AUTHORI 
OATS« 


<*«uai£ 


IEL3 


A(*9>*9) 


*ANS    ZE3NSR 
JtC=19E*     jl.    l»«0 

EXPLANATION 

'W.1BER    OF    ?OUArtO*S 

CONTROL    »»«A«STE»»    RETURNEO    TO    THE 
CALLING    ROUTINE. 

•I    GAUSS     JORDA*    elINlNATI»l    WAS 
MRFORNEO    C3RRSCTLT. 


<0    ERROR    cNCOU*TEREO    OURIrM 
•  OR  DA*    ELIMINATION. 


SYS  TEN    1ATRIX    A 

)i;*r  iano  itOE  yector  3 


1AUSS 


ttre 

?»TE5c5 
I*TeGeR 


"?1L 
»EAL 


IHRORTANTl    1ATRIX    A    AN!    »ECTOR    *   -(AVE    TO    SE    'ROVIOEO    MROUS* 
THE    LASELeO    CONNO*    4L0CK    *93.       4C1TN    i    »uo    4    .ILL 
WILL    UNOERSO    C*A*SES    0URIN5    EXECUTION!    I.e.     Trt6 
ELIIlATiaN. 


30 


39 


♦9 


30 


39 


60 


69 


COMMON    /    AS3      /    *#• 
REAL    Al*9t*<>),  J(S9> 

ERS-l.OCE-OC 

DO    06    I'Lintl 

30    02    <•!,»., 1 

IF    IK.eO.Il    GO    H    02 

C0NST»-4tK.I»/AU*ll 

»T«.j l«Ai«.ji.CO*ST«»(I.j» 

IF      U.cO.I)      A(*,JJ...1 

01    CONTINUE 

41*  I  "«1«  I  ♦CONST'K  II 
92   CONTINUE 

coNsr«A(i>i) 

90   Oi    J«1#N#1 

01  »itiji»»i:.ji/c:.iST 
A(i<i)«i. 

8(1)-»<I>/C0NST 


0* 


93 


00   96    j'l.nii 

IERROR-0 

00    3*    *»l#*»l 

IF    (A8S<A<J»«I).5T.£RSI 

!|R>aa«URROR«i 

CONTInUE 
F     H4S  HI  JM.iT.eRSI 


SO    TO    1* 


If  II6R* 


:o  TO  39 


ROR.iI.**U    CO    TO   97 
06    CONTINUE 

SO    TO    OS 

IEL3--1 

2EL33*!eL3 
Of    RETURN 

6*0 


9? 
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/♦/17'       JPT«2      .MONP  Hi    ■..•)♦<.■•<» 

1  JuiitOUTlNE    UlJ   (  i.  [TER.ITRNAi.SPSLani 

J  .      GAUSSIAN  ELIMINATION    *(T>«    LINEAR    ;3U»TtTt   STSTc*    V    (He    JIG-iT 

j  c  .    HkHit  siot  vector  i  is  intechateo  '«*,?*£*  dis  i»t   >r;re*  ••»- 

c  .     t*u  a   in  me  ir»tx  vi.     iter   is  r-ic  luntjci*  Or   itcrauons  «►- 

C  .      TURNED    ro    TIE    CALLING    PRIBRAM.       T4E    ITERATION    TO    IMP-IVES    ME 

c  .     solution  sa  hn»  Tries,   oiru   me  RcStllmi  error   :;  <imin 

C  .       tut    1A&NITU0E    OF    iPSLQN.       S-tJULO    THE    'tUrtRER    0?     1  T5*A  f I  OnS    :X- 

1C  C  .      CEgQ     [T8NAA,    THE    :  jJSOJf  lit    TcRmIHATC!     TrtC    ;  »LC'JL*I  :J-«    x.tit 

i  .      ScH    ITER    n   4    -HiiU/i    /»tue. 


19 


TO 


T5 


SO 


AUTNQRSi         KONRAO    VOSELER.    -i»*S    :E3NE* 
OATEl  OECeNiER    31.     I >«0 

.ARIASLt         EXPLANATION  rtPI 

■hjnSe***?  "fiilf  ions •••••       t5fj*,£j 


C  .       ITER  NU1SER    OF    ITERATIONS    TD    BE    'ETUPNEO  INTEGER 

20  C  »0    GAJSSIAN    cLINlNAnnM    PEP.FOPNEO 

C  .  CORRECTLY. 

.  <0   NO   CJNVE^SgiCE    IF    I'M?    MLUTta«    FOR 

IHS    DESIRED    SSSgLUTIJi'i    £?SLON. 

.  Iteration  3to  rl^nr-urn  after 

itriax  iteration  src». 

.       ITRNAX  MAXl-UN    NUN8ER    OF    I T fertATj  J»l    STfcPS  INTEGER 

.      EFSLON  OSilReO    RESOLUTION    IF    T  (E    SOLUTION  *cAL 

!      A<44>4«l         SYSTEM    MATRIX    A  REM. 

30  •      S«**»  »Iinr    -UNO    UOE    VECTOR    9  REAL 

I       INFOtTANTl    MATRIX    A    4N0    VECTOR    9   H4VC    TO    9t    PROVtOcO    MROUGN 

THE    LA9ELE0   COMMON    9L0CX    A93.       901-1    i    »10    <    -ILL 
.  JILL    JNOERSJ    .MANGES    OURXMG    cXECUTION,    I.E.     Mt 

35  C  ELlNlATIlN. 

C  . 

C3H1QN     /     A63    /     V 

RcAL     Hot.  ;(,).<;:90I.A<4?.>0>.400) 
oo  LOGICAL    IFAInT 

DATA     IPRINT    /.FALSE./ 
OAT  A    PI    /J. I*139i30*/ 


C  FORMATS    IELi    START 

♦  9  601    FORMAT    (///3*M      1ATRIX    A   JNO    VERTOR    3    VOF.    ELININAT ION/ I 

402    FORMAT    (  (  oi,  i  j(9«»  [2  )  >  I 

603    FORMAT     I 2 1  ,  12,  13  (  1  J.  1  ><■  9.2  1  /  f(  01 ,  1J(  IX,  .PS9.  2  )  I)  t 
40*    FORMAT    I///23N      »lVOf  iStE'TES    STSTeN    IN»;3#i"»N.     I  T£RAT  I3NSSCHR  t  T  T/ 

JO  405    FORNAT  (///2X, 12,241    «OEFFi:iENTEN    AUS    GAUSS/) 

603    FORMAT  ( I 21, 12, IX. l»= J. 2 ) ) 

40?    FGRNAt  I///3*-"      <LEI  ISTER    2E  tL£NVE*TaRrC1NITTwIN<6L    •.ei.2.3) 

60S    FCRNAT  ii?ttiZ-i      'eSIOIU*    •>E1>I3»3H    it.CI,  13.  L2H    ITFRaTIjnEnI 

»0<»    FORMAT  j///2d>U-l      -ErtLeR    V I  .  t  J  •  in,  U,  J-i  |    •,?ii.3/l 

35  410    FQR"AT  4///23M      >TSTSN    na:n   ELININ aTIOM /) 

C  FORNATi    1EL3    STOP 

L0»6 

40  IC»N-l 

C  RE4IIR4N6E    VECTOR    9    IN   ARRAT    V.    IF    N  IS    LESS    THAN    «9. 

IF    (N.E0.44)    SI    TO   32 
6)  00   01    I«I,N,1 

Ot    v«I»(il«v(  !,)U 


PRINT  STSTEN  MATRIX  A  ANO  RISMT  -IANO  SIOC  VECTOR  9  9EF0RE  THE 
ELIMINATION.  I?  IFRI1T  IS  SET  .TRUE.  . 


02 
03 

IF    (IPRINT)    *RITS     (L0» 
IF    (IPRINT)    <R!TE     (LJ. 
00   C3    [«1>™,1 
IF     (IPRINT)    *RITe    US. 
IF    (IPRINT)    -«a;    Hi, 

401) 
4021 

403  1 
402) 

(  J.J'l.N.i) 

t#IV(I»4l#J«:/M»ll 

(  J .  -  •  1  .  1 »  L  J 

LINEAR    EQUATION   STSTEfl 

V(N. 

1)»«(N)«V(fl,N»l) 

04 

OJ 

00  04  !•;,»>. i 

30   04    .-1,1.4 

All.  i  l«V(   l.Jl 

oa  o)  i-i.N.i 

1III-.0 

85  C      CNECH  FOR  STA9ILITY. 
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subroutine  tcLi  rsiin     3f»w     *«<}n*  p rt 

c 

Ft-.O 

00   ca    l*l,ii,l 
3*    SCI)«.0 
49  90   09    lM,".i 

00   07    »«l»N.I 

8J  i«i:!4llr«m,M* 

,9  H.^1'1""1 

00    -0    J«JJ.N,i 

00   u9    13*1»N,1 
39    S*»SK»v(I»I3l»V(J»I3l 
100  F»SK/<4(I  )»GU>> 

[F    (ABSIF).iT.»»StFII  1    '!•» 
10   CONTINUE 

^  I  -  *  3  S  «  -  I  1 


FJ.4C0SIFZ) 

109  f(«(f;«i>«.]/»! 


133 


160 


ST«»T    ITERATION    TO    F  I  *0    4*0    :i»»0»f.    lOLUTION. 


00    27    lX«l#«99.l 

no  irea-ux 


30    13    J-  I-  *»1 

FIVOTISATION    ANO    SAUSSlAN    STS»    FOR    CDLUnri    J. 


IIS  Jx.j 


DO    12    i»J.N,j 

IF    (A8S(V<I«JM.«T.*«S(EU    30    TO 

£0   TO    12 


120  11      i'l 

E«V(t«J> 

12  CONTINUE 
IF    (H.cJ.JI    31    TO    14 

oo  13  :mi».i 

129  iU)«*iu.;i 
<<U.I>«v<J,t) 

13  Kjpixu  I) 
Slftl*,0 

14  nn.j.t 

130  00    IS    I  —  -.n,i 
<•»(  luinij.il 

oo  :j   »*• jt» 

V(I>KK>«vi[,«K)-*<  J^ftKI** 
13    CONTINUE 
l*N 
IF    IV<t>I).ca..Ot    SO    T3    12 

i  SOLUTION. 

140  oo  :»  i«i»n»i 

16    Sill*. 0 

jin  )  •<(  N.m/y  m,n> 

oo  19  ;w»*»i 

T-N-1I 

149  IF    :vi  [,l  KE0..0)    33    TO    32 

J(I].»<I,«> 

00   17    JJ«l,K»l 

3-N-JJ 

IF    (J.E3.I)    33    T3    14 
130  17   Sit l«G« I)-««J>«V<I#JI 

15  CONTINUE 

19  ill  I  -til  il  'VII.  I) 

C3RRECT   30lUTt3N   -»r    »NOUNT   3F    ERROR;    HMTC*    HAS    rN<3lVE9    IN    r«t$ 
139  !  ITERATION    STEP.       TAKE    THIS    VECT3R    AS    START    VECTQS    ?*R    ME     1SXf 

C  ITERATION. 

00   20    l*l*N*l 

20  JlU-31  I» -3C> 

I  COTFUTE    THE   REGAINING   ERROR.  .» 

00  22    I-l»H.l 

y<:>Ni«.s 

169  00   21    J't.N.l 

21  ll{.fl!-VI 1,1). ilj !•»! [»3) 

2!  vil»n)«»i  ;./ii-<i  ;.N) 

R",0 

00     43      l-l.N.l 

170  23    R«R»VtI»NI»V(I#N) 


159 


SUBROUTINE    leLJ  7W175       WT-i      FN01F  FT*    ».j»^»a 

R<SQRT(R> 

IF    IR.LT.EFJLOUt     SO    TO    28 

if  u«.ir.  n»^i«i   id  rj  3i 

50   is    I«i.»»i 
ITS  2*    8(t)«SU> 

»«INI    JTSTEn   N4TRIX    »    »N0    RISHT    -UNO    SIOC    VECT1R    8    iflti    Ih£ 

c 

ISO  IF    IIMIND    WRITE     (Lit    4  JV I    L  * 


190 


200 


205 


Z10 


if  liniND  WRITE 
F  I  IFRINM  WRITE 
0   2)    l-i/N.l 


29    IF    (IPRINT)    .itTE     (l  J.     60JI     U(y(I.JI,J--.jH,i),«ii) 
IF    (IFRINT)    .9 [TE    ILJ.    6021     ( J.J-i.t.iJ 


lM  §8  IS  Wffct 

2ft    IMI»J>«4<I#J) 


5 


193  28    CGNTIMUE 


STOP    1. 00*    TO   FINO    SOLUTION. 
27   CONTINUI 

REARRANGE    SOLUTION    IN    THC    IRR4T    /. 


00   24    fl.N.l 

29  VtI>iOI«i(  II 

PRINT    SOLUTION*    IF    t»RINT    IS    SET    .TRUE.    . 

30  IF    (IPftINT)    WRITE  (LJ>  60S)    i 

IF    (IPRINTI    wRIT§  i\.Z.  60S)    (I»6(t>»:»i*N»l> 

!f  <iprint>   w*it2  ilj»  6on  fI 

IF    (IFRINT)    «*:Td  ILJ.  60S)    l.LX 

ORUCKEN   OfR    H4TRIJC    4     JtO    3ES    VECTORS    8    'UC*   LLI"  ("*  r  I  ON. 

IF     (IFRINTI     WRITE  (Lit  1UI 

IF    (IFKInT)    WRITE  ILJ.  6021    U.W«l»N»U 

So  U   t>i«Nti 

31  IF    (IFRIMT1    »«ITc  )LJ.  6031     ! . i 4 ( I , J ) , J .., *, ; ) 
IF    UfrInT)    WRITE  IUi  6021     <J»J«l»/»»lT 


RETURN 

ERROR    ENCOUNTEREO.       STOP. 
31    ITER«-l 


C 


«linit?i.n:ii' 


t.I.Wd.I) 

220  _  RETURN 

NO  CONVERGENCE  OF-  T*i    SOLUTION.' 


33  ITER— ITER 

225  RETURN 

END 
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FUNCTIdN    FJ  /W17i  li»W        ^10n>>                                                             FTN  l.ltlrt  zU\.Hc*,     I't.rd.'.J 

1  REAL    FUNCTION  F3    (  C  Zi  F.  «  I.  »  2>  «3«  3»  WflJP.Cf  2  I 

c  I  confute  function  <»i.j£  r»Fui.x2>  ?f  a  30-»0i.rNuni4i  jshg  1 

J                         C  .  Mt    tQSN6N  SOeNE    IN    :»£)?»    tl    [NCF.tAin    SFSfcO    >no    4C.jrf4Cr. 

c  I  iuthosi  uns   :;**?«  ! 

C  .  0AfE«  3eCEN*efi     U»    H»3  . 

10                             i  *  <4il48t.£  SXFIANATION  TTFE 

a  :  »«Ft»i  J«tfW^H"!5EFFii!?st» 5?ir**        : 

C  •  II  XI-V4LJE.    4HE«E    r    <£LL    3C    C1NFUTE0.  *E4L 

C  .  12  <2-V4C<JE»    **£«£    T    JILL    iC    CQrtFUTEO.  »c4L 

IS                          C  .  Na«0«l  1«JE«    JF    TH§    JOCfNONIAl    IN    T*e  I»re„Ei< 

C  .  (1    DIRECTION. 

C  .  N0«0«2  ")»3Ea    JF    MS    »3lTN0«Ul.    I'l    TIE  INTEGER 

C  .  ti   0I*SCTiOH. 

C  .  F3  IROlNAfs    3F    THE    J0-»0irs,1.1IAl.  >E4l 


20 


2J 


«E4l    C3EF(*9).3(4<»» 

ncqef  i«nqro»i»i 

nCOEF2«nO«OR2*1 

CONFUTE    COEFFICIENTS    1    3ASE3    ON    COEFFICIENTS    COEF    4N0    tl. 


30   02    I«l»NCOEFi»l 
IC«!»"»CCEF2 
30  r-c  JfcF«  IC  I 

00    Jl    Jl«l»K0R3R2#i 
J-NCaEF2-.il 
IC»<I-lT»tC0£F2»J 
01    T-C0EF(IC)»*2»f 

39  02  9(I)*r 

2  COAFUTc    1    9AS£3   JN   C3EFFU IE1TS    0    AHO    Xl 

T-OtNCOEFll 

»0  00     ;i     Ji-l.«0»1(U»i 

J«NC3£FI-J1 
03    T«0U)*X1»T 

Fl.T 

RCIU4N 
♦'  ENO 
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sioc*  o*r»  49s 


fttin     js»r  •^     ,»«onf 


FI"»    -..*••..:> 


>«/»</<:<».   ituiiH 


SIOCK    0*1*    «3» 
CannON    /    »s*         /    *.d 
'l»l     »<6»»0»t»8(o*l 
6  NO 


162 


SICCA    OAT*   OATkt,  !»M75       JPW       'ici»  FTN    ..(>♦»<■)  SC./t,i/iS.     1*. 

cannoN  /   oat**  /   iui:,ijif 

tiki     HI«i7i»l,lil',;iM.<JH.7.)l.fO,:,JI 
INO 


163 


>«OSK»K    0£.1.1»  ?»/li'5        >»X'i       ?*a**  rt*    «.«»*<*«  e.C/12/.I*.     ii.w.     I 

1  MQGRAM   OtllH    UMFUT.OorPUT.TAFEl.  7AF€6#'AFt7.74i>n»I.'«i>uTl 


this  is  4  oenoNsr*AnaN  frog^a*  amo  rr  shu**  the  correct  use 

OF    ?Hi    r?\.    uiiptr    JHHT    USING   SXAf1P1.ES. 


AUTHORI  -(A1S    IE*1€« 

OaU<  3ECe«S«    11.     1963 


10 


COAION    /    A«LO      /    PLOT* 
C0J1N0N   /    0«T*<.    /    <w<2>l)>r 


19  MAC    P10TR12J6) 

RtAC    IIU.  7»;T.x2<9,  MliDD./ill.rH,!,  5) 

REAL    CaiF<6*l»SU8«l(»9,*9I.JUM2<7»7l 
INTEGER    I»i*(4) 

20  LOGICAA,    I  HOT 


:s 


5»r»    ICnT      /10M»  / 

3»f4  ipiqt  /.true./ 

OAfA     PI  /3.1*1592oJ*/ 

formats  jcngA   ;twi 

181    FORNAT    I///72H      tein-      T-IIS    IS    •RO'SPAI    Sinn*    AND    <t»U    SEE.    *3<    7 
•0    USE     7H«     BARVEU0US/7UH       TPl     3I1ART    LlBPAPT.        TOUR     INrc-«t  ST    tflU 
'UR/25H       AMD    FAST;*    P»GJS AH1ING'  •  • 


Mi    GREATLY    AWA40E3    it    :1SU«/2!H       iiiO    FAST:*    P»GJS  »«1II*6«  /  /  /I 
rO««AT    <22H      CLOTTe*    11 ITlAt I2E0. > 
FQRBAT    (31M      3»TA    »3INT$    *£»3    FS0.1    7A(»C    7.  J 


10  102    FMfUT    <22h      clOTTc*    II jriAU  I2E3. ) 

133    FQRBAT    (31m      3»TA    '3isri    }£»3    FS0.1 
144  f  crha  r    i;,>,n<ii(».M,«-ui*ii.»i.  iH«i"t':,oA»;Hii.i4j.7«,,tr«i«.  /». 

•ArtTf14X/oA»7riPL  177 £  J  .9*  .X.F  13  ,  J I  /  9x  ,»mu;  ;P  >  o(  1  X.  F  It  ,  3  )  I 
109    FQ«rt*r    C7h      i>L37?E4    4.10   JSeR    4RES    JCrt'iEO    AMO    44ES    0*kmM.  I 
39  106    FQRKA7    (31H      Pt:ASJRE3    1ATA    SURFACES    3RA4N.) 

107    FORflAT    (7*M      e  t  T  §  *    040ER    JF    POlTNUnlAC    I't    »!-».  «2-    M9    Xl-OIRSCTZu 
•  N.    IPRtNT    Af<0    IS747I 


108    FORMAT  (2JM      SURFACE    FITS    :3NPt£7E0.) 

159    FqRftAT  (30H      4»P*3xInATE0    SURFACES    3P.AWM.  I 

401    fO»*AT  liHA/lHTI 

701    FORKAT  (A  1,44  "  " 

7J2    FQRrlAT  (3111) 

70J  FganAf  (»Fis. 


*9 


50 


791    FORraT    itl.aini 

fJRHAfS    3  End    STOP 


70 


*0 


I/O    ASSIGNMENTS. 

U'toGtuf  isHsro*"* 


m 


LP*1J 
99  ««17;     (II.    1011 

(13*    Mil 


4*ITe 


60  C  

C  '.     INITIALIZE   THE   PLOTTER. 

♦9  tf'.mSf  »*CAtL*tiifG*(L?r 

MtTS    (LI,    102) 


VARIATION   OF   FOLTHQTItAL    3ESREE.      LOOP   STARTS    HERE. 


T9  niJM«t 

nlJTOP* 

■   *2JTRT« 


*1S73P.« 
12STRT.1 


Willi 

38  to  :8282|:3HTl?::iiT3::t 

00  10  io<pR3->«3srar,<}sr3F,i 
Irn«!naaeii»n»<«jR3»2»iI»«no»OR3»i) 


99  IF   urn. .7. 4*1    jo   73   10 


164 


94 


•«aca»n  oeio*  74/175     3pt«2     pion#  fn  *.<>•**•  ic/i2/2v. 


C  .      16*0    *EASUI)eO    34T4    fail    TaP«7    4*0    St-tULTANfclUSLT    SeA«^    7H£ 

C  .       ilNInui    ANO    m(;vj1    /ALJE     ~.f    ;tC1    Xi,    «2    ANJ    r    i»cJuUdJ    TO 

c         .    sir  the  »lot  <inh  irti  soiwoauie;  5*  a  sheet  of  ?«»c-ii. 

w  c      ; '. 

UVINO    10 
IVAL        •     10M 
X1NIN     •    ♦<»99')<l. 

xihax    •  -<><»<><;<». 
loo  x2«ii«  •  ••»<>9'j,». 

X2HAX     •    •«««. 

thin     •   •«»«<). 

TIHX        •    -9<»9>»9. 

91  «£>o     no.    7ci)   ;v»l,;j»* 
109  IF  nviL.co.lciri   ;5  'a  oi 

DSCOOE    (<.S»702,I»*KI    Hf»NTSi,*Mm.«?NTS3 

oa  S3  IX*I'<*"*rlw 

so  03   I2^i,«i»nts2,i 
lio  ival«iom 

02   *t»o      no.    7ii)   ;v*i.;;t» 
if  i  iv»i.t9.ic.ir)   oa  rn  j; 
iECJJe    ioo.7C3,i*a*i    i:(il,i2,l3)»i2(t:.,:z,!3>»x3ui,i2»ii>,Tt  u,J 

115  IF    1X1(11,12,13!.  LT.Xi.1INl     IL.1I1.«l(Il,i;,I3( 

if  (rai#i2»i3i  .ir.rnjN  i  r»tt  •»<n» :.:.  Hi 
If  txi(Ii,i2,i3i.iT.xii4xi   xL.i»x«xi(;i>  :;,i3i 
if  (x!tII,i2,:3i.>iT.i2ixx>   «-J«4i.i;t!.,  :*, :  oi 

120  If    i»HJ#i2#i3»    .4T.T1AX    )     n»i    •rii;,;:/r3) 

93    CONTINUE 

•  «IT6     ILL    1031 


129 


130 


1M 
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OCFINc    PLOTTE*    A10    JJE'    »«£4S»       OR**    \X£St 


10  *   1>. 

TO  •      3. 

X1FNIN  •         .0 

H"U  •    10. W 

139  X2fHIN  •         .0 

i;»«u  •  io.o 

TfnlN  •         .6 

"»u  •    10.0 

HUfltN  •     IIMI 

1*0  Uu«»l  •    l!l4X 

X2UNIN  •    X21IN 

X2UN4X  •    12H4X 

TONIN  •    TtJN 

TU94X  •    T«AX 

1*9  ALfAXi  •    60.0 

ALFAX2  .J2C.3 

AlfAT  'UU.O 

fLOT«(*ll  •    <4L£4Xl.»T)/}3a. 

PLOTHOZJ  •    (4i.F4xi»(»ti/l35. 

190  CLJTa  i->3  I  •    IALFAT    •»I)/139. 

PL0T»I64I  •    10 

PL3T*<»JI  •    TO 

PUT*(<>6»  •     (J1P1AX  -TIPIT"*  )  MXlU.IAX-XJ'.MTNt 

PL0TAC»7>  .    (Xl»1IN»XlU«*X-Xlf14X»XlU.1I'l»/lxIu.14»-<t'J'"tHl 

199  PLJTalsei  •    (X2P14X  -X2P1I*  l/«»»,jl4x-«2jni»«) 

PL3T»<»9I  •    (X2f  1Il»t2Ul1*X-»2Pi»X»X2UN::ii /(X2U(1»j-«2uilNI 
r»*»x  -    fo>»!N  )/(    r<j*n-nj«(>»    I 

r»«I-|»TU«4X    -    'P»4X»TUr)IH     '/«     TUIAX-HJNll     ) 


PL3i»(7G( 

X1L  •     (XIP14X  -XlP^tN 

I2L  •    IX2P14X  -xJPIl.1 

TL  •    <     TP14X  -    T?m* 

*«ITE  (IJ,  1C4)  <lPHIN»Il»14X,X2PNIN,X2M*X,TfHN.TPPl*X,HU«tNi.XiU 
•■4X,i2u«fi.«2Ul4»,r'jij-»,rji4« 

if  (IPLOI)  CALL  AXES  (X3»T3»XlL,ALFAXl,2HXl#.2>.riUNIN.XiJlAX.*Hf 7. 
•1,1 til 

IF    (IPLOTI    CALL   AXES    (xa.f3»x2L#ALfAX2,2HX2,»2,I2U«lN»X2UNAX#*Hf 7. 


•3,7, <J> 

If  tin 

•3.7.11 


:plot>  call  axes    :xo.to»tl  .4LfAr  , ;ht   ,-i,tj*m  .tunax  .*h<?, 

170  iklii    III,    109) 


165 


175 


>*aZHH    J£tO»  rstili        IPT«2        IMONP  'l«    »..)»»■»*  n/,i/i«.     li. 


ORAM    1EASU«E«eNT    OATA    '>U»f'C£i. 


5S*o*"  IJil  ^M?»rs5 » 1*  •••••••••••••••••••••  •••••••• • ' 

88  %%  u:{:s»Hi:i 

IF    ({PLOT!    CALL    ruRTu    (  (HOT  «  TPLOT.XK  II*  12,  13  )»  X2(  U.  12,  131.  r  (1 1. 

if'iiplot  .ano.  t2.cO.li  :all  put  (xpl:it,tplot,3i 
us  99  if  iiplot  .aho.  u.if.i)  call  pljt  (xptar, rnar.ii 

00   OS    IJ«l,NPNTS3.l 
00   QJ    I2«l,NPNTi2,l 

30  a)   iw-npntsi,  i 

If  (man  call  mrtw  <xPL3T.rPiOT,xi<:i,i:,i3J,x2(ii.i2>i3i.rm. 

if  (Jplot  .ano.  u.eo.ii  call  plot  (xplot.tplot.ji 
9)  if  tmar  .ano.  u.si.ii  call  plot  ixPLiT.n»L3T»2i 
»RtTE    ILL    1061 


no 


190 


199 


:oo 


:zj 


290 


299 


CALCULATE    SJRFaCE    FITS    TMROUS*    JATA    POINTS. 


06  SJfti'ICIJ'ISff' 

REAO       (    4,         •) 


.Oa3«l,NG»0«i.iO«0»3f  IPPI'IT.ISTAT 


NdlORl>»3XORl 
209  NOROR2«»<3»OR2 

NaiMlxioiioii] 

iprinNq 

ISTAT-1 

NN*M«N 

IF    (L.LT.l    .a*.    L.ST.  2)  S3  T3    36 

IF    (H.4.T.1    .OR.    i.il.2)  SO  r5    3* 

219  IF    (N.lT.I    .0*.   1.GT.71  13  rj   gi 

IF  il'in.lt.i   .jr.  L«*.-i.;r  .»»i  so  to  36 

•  C4I»7.     ?4IJ0ij'>NT51"'<*HrS2.'«P'«TS3.L»n»'l,C3CF.IF*INT,IST»T.SU8«t#^N.S 

IF    (SERRQR.LT.OI    SO    TO    99 

220  .HITS    (Ll«    ICO) 


\      : 

C  .      ORAM    APPROXIMATED    SURFACES. 


5a*5Hi:t;s«f??;i 

210  00   07    U»l»NPNTSl»l 

30   37    l2»l.NPNTS2»i 

TArF«F*(catF,xiii:,;:. i3i. «2ii:j£2,i:),x3  hi. ii.i3i.  nor  o«  i.«o»o«2. 

•MRQR3I 
If    Tt  PLOT  I    CALL    THRTrf    t  XPLOT.  TPL3T.X11  1 1,  12,  I  3  t  .  12  I  11/  12.  I  3  I.  t  APP  t 
239  IF    llPLOT    .ANO.    I2.c3.il    CALL    'LOT    tXPLTT,  rPLOT.  3 1 

07  iS  oipVSIi;w8hSfiiT,u  :iLt  "-or  lx'taT'T"-0T»*' 

30  08  I2«1,NPNTS2,1 
30  04  tl'I.NPHTSl.I 
2*0  TAPP.f»<C0EF.xlUl»IZ»I3l,x2C  11,12,13),  X3(  U»  12,  111  .NOR0R1.  *OR0«2» 

_0T)    CALL    THRTW    I  XPLOT,  rPLOT.X'f  II,  12,  I  3)  »X2  (  U,I  2,  13  »,  T  APP  » 

(IPLOT    .ANO.    Il.E3.il    :*LL    PLOT    (XPL2T,fPL3T,31 
OS    IF    (JPLOT    .ANO.    il.ir.il    CALL    PLOT    ( XPLOT,  TPL0T.2) 


F    (tPL 

RlicV 
9«    IF    UPLOT)    CALL   ~«UF4Si 


«»    I?    (IPLOT     .ANO. 
2*9  YRITE    (LI.    1091 

99    IF    UPL0T1    CAL 
10   CONTINUE 


TERMINATE    SRAFICS. 
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mOGZM    06*0*  7*/ 173       JPT-2       ^lO.V 

If    <mOT)    CALL    STO»S 


FTN    ».o»-.Jj  H/ii/;*.     H.-OOd  J»*c  ■> 


2tO 


STO»    0077 
ENO 
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?*/17J       3PT«>       PiOrtP 


FTT»  *.J»»N9 


1 
5 

10 

1J 

zo 

?5 
30 
3) 
♦  0 
49 
50 
JJ 
40 
6S 
70 
75 


•  9 


SUi routine   «»r<.    MP-.r5i.NPNrs2»NPNrs3.i,.i.n.c.j£F.iP«lNr,isrAT.Jj»-« 


PRESET    STSTEN   NATPIX    A    ANO    *IGmT   HAND    S  IOE    VECTO*    i    P-3R    4 
40-P0LTM3NIAL    APPROXIMATION.        TM£    uIMfA*    F9UATIJN    jfiTel    THeK 
IS    SOLVED    it    USING    4    lAjiS    JOROAN    ELIMINATION. 


author  • 

0*TE  I 


-<AN>    IE9NEB 
OECENBE*    U»     1900 


/AilAilE         EiPLA.naMON 


rr»t 


NPNTS1 
NPNTS2 

NPNlSJ 

COEF(a*i 

L 


IPRINT 


ISTAT 


18ER    OF    C3SFF 

lit*  if  c3eff 


NUN9ER    3F    3ATA 

turtle*     OF     3lt« 

nunser    :f    jata 
»0LfN3nut  :jef 

NUN9ER 

sail  4 

HON  t  AC     IN     <>     31 

Huia£»   op  coeff 

iOnIal    IN    (3    31 

CONTROL    JUANTIT 

•Z    PRINT    cOUATI 

4140    AFTE4    TH 

•1  >*Nf   ;3u4r: 

TH«    GAUSS    J3 
XS    PRfNT    *4SIC 

ACTUAL  P4RAN 
<0  NO  JUTPUT  4T 
CONTROL  STATIST 
/UlAlt'l  PRINT 
ANO  4 
AtO    3 


NTS  AT  CONSTANT  II 
NTS  AT  CONSTA.NT  t? 
NTS    AT    CrlnSTANT    «3 

IENTS 

6NTS    Ff1R    THE    POLT- 

TION. 

snts  far  the  polt- 

TtQrt. 

e«TS    POP.    THE    POLT- 

tiom. 

S    'PIN*    3UT 
STSTjl    9CFURE 

auss  jcipoam 


/JJTAr/«l    ONLT     P 
ISTAT>1         GET    A8 

ISTAMO         SET    (ELATIVE    sPPJP, 
SUfUIRNXMI    5U3.NATXIX    I.    CLASl    II    t  (C    STSTEN 


ON    SYSTEM    AFTER 

roan  ELjNiriATlnn. 

iOUATIONS    WITH 
ETERS. 

sic*s"ro  3E  printed 

IROtNATg!    1«    OAT  A 

•P»OXI"ATfO     POINTS 

EVI4TI0MS    (ERROR). 

PRINT    EPRO". 
SSOLur!    EPPOP. 


INTEGER 
[NT?  JER 

RtAL 
INTEGER 

INTEGER 

INTEGER 

INTEGER 


INTEGER 


DM 

SUSN2IN>" 


IERROR 


11(9.7,5) 
12(9.7.}) 

?!«»}«' 

»(»kfi*l 

S14+I 
MPORTAftTl 


matrix. 


0INENSI3MIN5 
'ubNAftJiX     " 


PtAHLE    POP    SUBNl    ■••< 

ass  In  '■hC  ;jaKAT*iA- 


2.    :l 

1.    CLASS. 

ion  =la<  «efjfiN£o  ro  rue  calling 

PROGRAN. 

•♦L    NO    ERROR    It 


•<CAC 


INTEGER 
ReAt 


I.NTEGeR 


TSCTEO.       GOOO    J09. 
CUREO     PI    !■!«     jAUSSIAN 
i      COEFFICIENTS    ARE    NOT 
REVIS*    PAP.AT.ETeilS. 


•-1  4*  ERROR  3C 
eLlNIN4TI3l 
ACCURATE. 

U-C33«OINATeS    3P    THE    CaTA    POINTS 
«*-C30P3tNATES    3F    TIE   DATA    POINTS 

»3-c33s<9iNAr;s  SP-  f><E  3aTa  points 
T-caoRorUTsS  :p  m«  oata  points 
irsTEi  lAiatt. 

PICHT    -IANO   S13E    VECTOR. 


5EAL 
<tAL 

<t4L 
JiAL 

*SAL 
-i»L 


THP     OATA     'OTMTS     «l.     f?.     JI    4M0 

VIOED    T-tPOuiN    f^E    LASCLEe    " 


US 


VISES 

THE    STSTEN    1ATJ  IX    4    ANO 

9    AP.E    AVAIL43Lt    TNPOuG.I    rut 

*3*.     this  saves  a  lot  op  ?-*at 

CALLEO 


_    T   HAVE    TO    9e    P*0- 
OnnON    jl;*<    jil»»    anq 


PI3HT    hanO    SIDE    VECTOR 
L4Sikl9.5'j?!'J'».4jrOCR 


EL59    : 
'RECIO 


OUS     IrtlN'i 


COKNQN    /    0aT4*    /    t\,i:,ii,r 

cannon  /   as*       i  a«4 

PEAL     XJ.(9,7.J),l*[?,7,JI,l3C9.7.H,rt«»7,  «l 
REAL    A<6A<64I>345*I 

■  EAL    C3EP  tf-),  3UB/"liHN./lNl.  SJ8NIIN„N),STAT(9,  7,  51 

C  rORNATS    NAT4    START 

101    FORMAT    I///7M      .rTSP,lJ,*H    STEPS    NO   SOLUTION    \IITKI»    THE    4CCUPACT 

•0P.cl0.3»iiH    |       I    JUIT./I 
401    fO*HAT    1///51H      3R06R    3£_4PP»IJX  I'ATIpN    »gLTNO«TAL    IN    «1    OIBeCTIONt 
•»?»I3/SlH_    0P3ER    OF    APP»0tI.N*rijN    »?tTMONliL    irt    <?    3I3cC  f  ION.  J  t,  13 


•  /40H 


OPOER    OP    APPROXIMATION    P0LTN3NIAL    IN    (3    3IREC T I 3N , >X ,  13 


ici 


NUN3ER    OP    OATA    POINTS    JITH    CCNJTA'IT    XI,  l«Xj<}/4>9H      hij-he*    OP 
3ATA    POINTS    <IT-t    CONSTANT    X2  #iJX.  iJ/^OH       <U.N8ER    UP    OaTA    POINT;    .i 
•Th    CONSTANT    <3,l*X#I3/> 
MX    FORNAI    1///62H      iTSTEN    1ATRIX    4    4N0   VECTOR    3    3EF0RE    GAUSS    J3P0AN    i 
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subroutine  i»r*  7*/i7i     irr*i    »momr  ftn  <..o»*,<i 

•  LIMITATION/! 

i33  =:««»r     (I  *X.13(8«.  12)  I  I 

604  FgR*AT    I2x»12#i3(l<»i?64.4l»ix/(c«.-3(ir,  1P£9.2))I  | 

603  FORMAT    <///***      cOJaTIJN    SfSTEt    xFTET    iiMS*    JOROAN    tLIMHI[]<4/l 

1  606  FORMAT    (///2bH      COEFFICIENTS    Ca£F(I*J»lU/J 

»07  format  (2*.iMi  ••  I?) 

605  FORMAT    (12».  3H<    .,12,  UI9X.I2I.  U/ 3<  7X#  Ul  9X»  I2l/)l 


99 


130 


144 
1*1 


609  FGRNAT  I2J.3HJ    •.[2,i2(Ut:<'.3».I.X/3(»ir,;2llX/C*.3)/)( 

610  format  l///3<>-<    jRofurEs  if  oata  »o1nts     -    *s  ii<ENi 

411    FORMAT  t/9M      T(I,j,.  (2»  Ml  » 

61*  Fjkmj  i///*6rt    jrolsates  3f  jata  »o;nt;     -    i 


61 

1    FORMAT     l/Jn       f(I»...[i.MM 

AS  APPROXIMATED) 

•il*    FORMAT     1/T.ln      jRRJ*     (l ,  j .  .  I  2  .  1*  I  I 

6ij    FORMAT    l///3tM      •iiLATjye    ERROR    «T    EACH    PflTNT    !•)    I) 
1M  616    FORMAT    i///2lH      jTiriSftC*    3  JMMA*r  I  //  7X#:,3HAVERASfc     jRiO».2aX.lPi9. 

•Z/7(> llHMAXIflun   ERROR    AT    I«.I3»3X»2MJ«>  C  l.  3X>2  <K>.C3.1<  »i?E?.2  I 
617    FORMAT    (///21M      STATISTICS    3  JMMART  t // 7X#.3'lAV£i>.  AiE    i<i*  OR  ,  is*,  >»  F«. 
•2»£H    t/7x,HHHAx;,1ui    rRJOR    AT    I  •.  i  3#?X#.MJ»»i3j  3A.2  H*«»  13  •  lit  i  PF,. 

109  613    FORMAT    <///7H      AFTER.  li.AH    STEPS    MO    SOLUTION    UtTMlN    THE    ACCJRACT 

•  0F.E10.3.UH    .       I     VJIT./I 
C  FORMATS    MAT2    STOP 


110  v?T5 


LI-LOGL'JI  ISESSMI 


ERROR-l 
TRMAXOO 
PSL0N».iOOE— 707 


119  PRINT    iPPROXIMATION    PARAMETERS.     IF    IRRtMT    IS    CREATE*    T-tAN    J. 

120 


NORORl«l-l 

-«0«2«M-I 


NORC 

N0R0«3«M-1 

IF-  TlPRINr.LT.O)    30    TO    51 

««I(c     U3.    6011    iO^ORl.  t'3«OR2.»OROR3»'4F«rS:,NPNl,2.  MPNT  ",  ) 


START    TO    PRESET   ESSE    SECTION    SU3N4TRICES    1.    CLASS. 


129  31    I»OKR1«0 

ll«Z»L-l 


1N«M«N 

D 


TOP    eOSE    SECTION    SUdMATRICES    1.    CLASS 


ISTART-1 

IST0R«MN 

US  4START.(i«-l)«MN»l 


JSTOP««»MN 
SO    TO    01 

right  1am0  eose  section  subnatrices  1.  class. 
02  :start.<*-li»mn»i 

ISTOR«<  U-U»l)»i* 
JSTART»<l-U»NM»l 

.)STOP«L»NN 

PRESET    EOOE    SECTION    SUiMATRlCcS    1.    CLAIS    IN    STSTEM    MATRIX    JSINS    *AT<>1. 


03    CALL    MA7A1    (NPNTSl.NPNT52#NP1TS3»N,N»tPJHRl»SU8Nl,NN,  SU3M2) 

n«o 

190  30   0*    I-LSTART.ISTOP.l 

U-Il»l 

000*    J-JSTART.JSTOP.l 
J101»l 
199  9*    A(T,J  ••SU8MKI1, Jl) 

09    IPO*Rl»lf>Oi»Rl»I 

CORY    OEFINEO    SU3MATRICES    1.    CLASS    0I60MALLT. 

1*0  Ll-L-1 

S8 «  $i:$Jthx 

.     IST»RT»(I1-1I»MN»1 
IStOP»ll«MN 
1*9  JSTART'IJl-lMMN*! 

JST0R«41»MM 

00   06    I«ISTART.ISTOI>.L 

30    06    J«JSTARf.JST3».l 

36    h:hi.j->i.im;:,,i 
7    CONTINUE 
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SUBROUTINE    «AT*  7*/ 1 75       1FT-2      f*On*  Fl*    '..o-.'.j 

I  FRESST    RISHT    *A*3    SIOE    VECTOR    S. 

C  1-0 

179  00   OS    (.1-1,1,1 

iraMAi«Li-i 

oo  os  ii»i.«.i 

IP0Ui4«r11-l 

00   08    Xl-l»*»l 

ko  iro»ft}«Ni-i 

08    Jti)isMXFNTSl,»«F,1TS2»«FNTS1.IRa«»Rl#:FO\tf,.2»IFailR3»l> 

Nl«t*NN 

lis  mint  srsTtn  natrix  ino  ris*t  -uno  sioc  vector  a.  if  i»«m 

C  IS   SET    TO    2. 


If    (IFRINT.NE.2)     iO    TO    10 
HRITE     110,     i02( 
190  •kill      10,    <>03>    <J,J-1,>U«1) 

30    04    t-I»Nl»l 

09  <«,[£   ilo,  so*)  i.nii.Ji.j.uu.u.ei:] 

»aITE    (Lit.    1031    IJiJ.i.Klnl 
IFRINT-l 


199 


21C 


:zc 


C  SAUSS    JCAOAN    ELI^IHATIOH. 

10    6M10N.«PJ10N»N1 

CALL    IfL*     (M1»ITSR»ITR'1AA.S»SL3K> 
200  If    UTeR.LT.C)    S3    TJ    32 

mint  souiuaw  srsTSdi  if  ifrint  is  :et  to  i, 

IF    (IFRINT. NE.i)     U    TO    12 
209  HITS    ILO,    «C:T 

«RlTl    (LO.    *C3>    < J,J-l,*l.l) 


00     11     I»l»Nl,l 

11  MUTE    <L0»    »0<»l    IilllIiJIiJililIiUillI) 
ittlTE    (10.    903)    CJ.4-UHl.ll 

OEFINE    COEFFICIENTS    4*0    MIST    THEN,    ;F    :»»Ip»T    IS    ;»fc»fO    rn»™    3. 

12  CONTINUE 


00    13    Ll-l»Nt.l 

21S  13  cotrk!iM(i:i 

IF    (IF«INT.l.r.'))    iO    TO    l» 
•HITS    (L0»    606) 


00     15      Il'l.L.l 

vrite   <lo.  son   ii 

•  RITE     (lu.     60S)     (  13,  13-1,4.1) 
30    l*    12-1, n.l 
ICS'»T-<ll-l)»1N»t  12-lMN-l 
ICST0F-tll-l)»1N*l2«M 

;*  vrite  no,  604)  i2,(C3EF<t:..i:-icsTr-r, tcsTa?,!, 

229  19    .Rite    (LO,    SOS)    (!j.;3-l..<.ll 

C  STATISTICS.    IF    ISTAT    antHHG   ELSE   JUT   (I. 

16    IF    (ISTAT. EO.O)    S3    HI    3* 

210  IF    (fSTAT.LT.O)    ::j»t»-i 

IF    (JSTaT.GT.O)    iS»*0R-2 

tsTAf*iAisasTAn 
confute.  »esoiyTE  err3r. 

239  C      C  -(AFFROXIHAleO  VALUE   -   JIVE*  VALUE)  ] 

OO  17  I2«l»NFNt52,I 
00  17  I3-l»«FNTS3.l 
2*0  17  STAf(Il»{2,I3)--*(C1SF,<l(Il,  12,  I  3  ),X2  (  II,  12,  13  ) .  «3  til,  12  ,  t  3)  ,«J*0 

•il»r»0*5S2»NQ*0»l)-f(U»l2,l35 


C  CONFUTE    RELATIVE    ERR3R#    IF    ISTAT    44S    LSST     THAN    ZERO,    i.e.     iESSja-l. 

2*9  C  C    •(AFFROXIIIAtEQ    VALUE      -      SIVEI    VALUE  1 /t 5IVEN    VALUE)     ] 

If    1 1 ERROR. HE. 1)    30   TO   14 
30    IS    ll-L,r<f"mi,L 

290  '  88  S  I:tafii: 

IS  STAT  I II*  1 2#13) -STAT  til,  12,  1 1 1 /Tt  U,IZ,  13) 


299  C 


C  ERROR    STATISTICS!       CONFUTE    AVERAGE    ERROR    ANP    FINO    ISCATI3* 

C  OF    1AXINUN    cRRJJW 
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SUSROUTINE    *»T*  7A/171       3PW       »H0NP  FTN    <i.a»*'»4  c<./litl*.     i«.:,.<i  >A»; 

r 
19    ERRAVE'.O 

i«mi»x»*gs<ST*T«i*x#i)i 

INA*»1 

2M  <»«"•! 

30   20    IW.NPNTS1.1 

oa  2o  i2«i».*PNrs;,i 
so  20  {3«i.NP*)rs3.i 
if   nasi5FMm,;;.;)n.u.c>«Mii  si  ro  :o 

269  eRRNAX.STAKIl,I2,I3) 


2TC 


275 


HO 


320 


130 


1M1< 
KNA1>I3 

eRRAVE«EScRA»E/(NPNrSi«N»NrS2«NPNTS31 
If    »  IST»r..-.€.i  )    JO    fO    11 

PRINT    3R0INATES    0*    OAT*    POINTS    ANO    APPPOJUHATSO    /ALUES. 
IF    /  15  TAT/     IS    tOJAL    2. 


WRITE    (L3.    610) 

00   22    I3-1.NPNTS3.1 

write.   a;>  din   13 

wRlte    <l3,   ao3)   I  12.  I2«i»»phts2»II 

2(0  30    21  IW.NPNTS1.1 

21  WRITE  (to.    %♦!    [l.(r<I1.12.I3>»I2»:iMPtlTS2*ll 

22  WRITE  (LO.    403)    (  12 »  12*  i  .iPt r 52 , 1 1 
<RITE  (LO.    012) 

00    2*    I3«l»Ni>nrS3.  1 
2(9  WRITE    (LO.    oil)    13 

WRITS     <L3»    tCJ)     I  I2j 12-1.XP1T52.1) 
00    '3    Il'I.NPNTSl.r 

23  WRITE    (LO.     bG*l    Tl.lFMCISF.  11(11.  I2»  I  3  >»  *2(  11.  12.13).  UCl.12.  13) 
•.xgR0Rl,M0FOR2.N3JgR3).[2«l.N»NT3  2#l) 

290  2*    WRITE    (LO.    6CJ>    (  I  2. 1 1  ■  i ,  <•>  (  r  S2  .  i  J 

ISTAT-l 

S  PRINT    ERRORS.     IF    /ISTAT/    WAS    1   OR    2. 

299  29    IF    (ISTAT. NE.l)    33    TO    31 

IF    I  ItRHuR.cU.l)     :0    T]    29 


300  WRITfc    (LO.    613) 

"    J?    I3»l 


PRINT    A8SCLUTE    ERRORS. 

WRITE    (LO.    613) 

30   27    I3»1.*PNTS3.  1 


WRITE    (LO.    61*>    13 

»*IT£    (LO.    6C3)    I  I2.I2«I.WPNTS2»1) 

00    24    2l*l.NPNtst.r 
2ft    WRITeCLQ.     SO* i     II , ( STaT I II. 12. 13 ) . 1 
27    WRITE    (LO.    4031    (  I2»12«l.HPNTS2.l) 

SO  t5  31 


309  2ft  wRITelLO.  so*  i  II  .  (  STaT  ( II,  12.  13  )  .  I2«l.'t»NTS2»i) 


PRINT  RELATIVE  ERRORS. 


28  WRITE  (L3>  419) 

00  30  I3«1»hpnTS3»1 

WRITE  (LO.  9l*l  [3 

WRITE  110.  4t,3)  (  I2.I2»l»<PNTS2.1l 

119  00  29  Il-l#NP!iT3:,l 


29  WRITE  (LO.  60*1  1 1.  (  STAT1  U.  1  2.  13  )  »  I2«  1/')PNTS2*  1) 

30  WRITE  (LO.  903)  U2»i2»l»«PlTS2.1) 

PRINT    ERROR    STATISTICS. 


11  IF  (IERR0R.E3.2)  WRIT;  a:.  Ms)  :"»vfi!,">J,»'»',»'^",« 
IF  (  ItRROR.cJ.l)  wRITe  ILJ.  417)  :*»A/t..'»ii.ii>i«.i".c<<iii 
S3    TO    3* 

329  2  FATAL    ERROR   OETECTEO.      OUTPUT    ERROR    1ES:aGE    ANO    PRESET    COEFFICIENTS. 

C 

32    ITER  —  ITER 

[ERRCR—tERRQR 

3ll)«r<l/I.iI 


31    S?l|?..;S^"'1'1 

WRITE    1LI»    101)    ITER.EP! 


319  C 

C  RETURN. 

34    RETURN 
END 


i rite  aa*  6ia)  ;ter.cPsl3n 
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1  SUBROUTINE    NATV1     (sursi.NMrv.nJHT^.d.N.liij^i.i^M.n.'iJ^.'l 

t  preset  syan»r*ii  i.  aiss  sjsni  im  run  cystca  naTrii  »  fjk  a     1 

J  C  .  ao-poltnoiual   approxJiia  r[3N. 

c         l  turxoAt        -)4M|  zfSNti 

C  .  OATEi  :t:£n»t»    H.     198a 

10  '.  *A«l*aie         BlPLANATION  TTM 

:  ;«•!•••    suiii'SP'urt'iAt^rivsiMSfuit'fi    •«»*$;« 

.  np*ts2  nuiher  op  oata  joints  ;.r  constant  «2       iNreii* 

.  NPNTS3  NU.IdER    If    OATA    'OIlTi   A[    CUNSTiNf    (3         IHTjiSa  . 

19  t  .  n  tu/nca  qf  coefficients  ?3R  fng   poit-       itreii* 

.  NONIAL  IN  tl    DIRECTION, 

.  *  NUN8EA  OP  CJePrlCIENTS  rOR  THE  ?OLT-    IHTiieR      . 

-laiiit   i.i   u   omcriON.  . 

.  ipowri  ?otf£*  for  r^s  u  values  in  the  intej*» 

:o  6  .  3uni»rn«. 

.  SUIftllMUHNl    5U8N4NU    I.    CLASS    I*    T  l»    STJTSN  >EAk 

ItUK.  . 

.  AN  OtlENSIjNlNO    yA«[A«LE    fOP.    £UOB(.*H*N         tNTjac* 

.  SU8N2(N«N)     iu*^»r*Il    2.    CLiSS    IN    THE    ..JdHATPU  JeAL 

29  .  1.  :uii. 


JJ 


5J 


-■0 
30  *W«"-1 


«t»i.    Vj«nl<nn.iNi,  suinKn.ii 

traeil«c 


tf  ?UtWXSJ   T3    01 

T»    E06C    SECTION    SU8NATRICES    2.    CLASS. 


ISTART.l 

ISTOP-N 
jSTart.<<-i)»n»1 

jstop«ii«n 

*C  SO    TO   02 


RISMT   HANO   EOGS   SECTION    SJSNATRlCsS   Z.    CLASS. 
01    ISTART«CK-N>«N»1 


*9  ISTaP«UK»l)-N>»N  \ 

Jst*»r.((i-ii«N*i. 

•     JSTOR«N»N 

C  PRESET    cOSt    SECTION    SU8NATRICES    2.    CUSS    IN    SU8NATRH    1.    CLASS 

50  C  USIN6   IUTA2. 

02   C»LL    «*T*a    (NPNTSi»N»NTS2»NPNTS3.M»lra»«Rl,IPawR2»SU»N2) 

iaai  :iii;»»t,i;r]#,i 


ji»ji»i 


START, JSTOP.l 

,USNl(I»Jl»5UgN2C  Il»Jll 
60  J*    lP0*R2«ti>aMR2»l 


ccrt  oErtNco  susnatrices  z.  class  ousonallt. 


89  00    06    tl»l»«l»l 

00   06    JW.fl.l 

I5rof»«ii*N 

JSTART»<J1-1)»N»L 

TO  JST0P«J1»N 

og  os  i-ista«t,i5t:p,i 

09    Susnii  !-N,j-»|.ijWil  I.  Ji 
06   CONTINUE 
T9  *ETUAN 

END 
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T*/i7S     jti'i     »io*»  ft*  <,ti)«i  ;:i\.in*.    u« 

l  suaaouriNE  «4T<.2  iMP*rsi<*i»Nr52»H»'<T'?2,'i,if>aj*;»tPCH««<:.;j-jflZ) 

c  !    heset  su8«»r»r»  z.  cuss  sjjij  r-i  ne  3U9h»t«u  i.  cuss  1 

2  I     tuTMOti         i»ss   :?j,*e»  ! 

10  {  '•       VA«I»IIE         E»»lAa4TI0N  TY»t  I 

2  .     *p*rs2  su«9E*  a?  oir*  >:n>«rs  \r   :afisT*nr  »z       [NTe>£« 

C  •      NENTSI  HU13ES    3F    34T4    eaf.lfi    »T    CONSTANT    <3  !•»  Ti  vie »  . 

19  .      *  *J!M6»    IP    C3cP«  IC IE »* Tt    r"HP.    TuE    "Qtr-         [NTEui« 

c  .  -<Qiial  cm  u  3i»6cnon. 

C  .   tPOwdi      l>O.S«  ?-3«  THE  '1  VAlutS  !H  rHE  IMTcSeR 

c  .    ifOKgj  ton*  -~>i  rut  (2  v4lue:  :m  the  intese* 

20  c  .  jui.nrijN. 

2  .     subh2(n«ni   suaiirau  2..  :i*si   :n  r-ic  :usb*trix         »eai 

c  .  1.    CLASS* 

23  «E**L  "Iu6A5I»»Al 

[■au«3>a 


30 


JO   03    "lol.l 

IF    IK. 61. Ml    SO    T3    Ql 

TV    e0«*    SECTIT*    ELslEtTS 


35  .0    TO   02 

c         «isht  h*«o  eoss  seem*  clemnts. 
c 

01  I»««-*l»l 

C  nCSiT    SECTION    ELEMENTS    :■«    Sa8«A.T*I!l    2.    CLASS    JSINS    j*. 

02  ;u«B2(itj).s«ci»-<r;i.tF'«r?2.>««'«r5j.iPjH(«:/  i'a.m.i'i^i.ji 
*9                         33   ;»a««3«i»o»«i»i 

Ni.f»-1 

ao  :»   [-i,*.,! 
oo  o*  j«e»"». 

3*    SO«n2(l»l>.-i).SJS12t  I.  Jl 
50  •ETURN 

ENS 
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7*/i7s      0PT»2      pho.ip 


FTn    i.jhii 


SO/li/C*.     11..5. 


1 
) 

10 

15 

20 

15 

30 

35 

♦0 

«5 

50 

55 

M 


real  function  :*  {npntsi»npnts2.npnts3»ipimri, ipohr2.ipojr3.it> 


perforn  suhnations  ro  preset  me  cleaents  au.ji    in  r*: 

STSTEN  ,i»ri»  IX  »  ANO  THE  cLE'SNTS  »II1  !M  THO  JlirtT  ha»0 
SlOE  VECTOR  J  FOR  THE  LINEAR  E0UATI1M  irsr^l,  THAI  IUjI 
it    501*60   FOR    A    POLTNONIAL    O-APPP.flXIrtAT  ION. 


AU  THOR I 
OATEI 

variable 

npntsI** 
npnisj! 

NPNTS3 
IPQ»«1 

IP0MR2 

IP0WR3 
IT 


s* 

Xl(9.7»J> 

13(9.7, 51 

r«9. 7.  »i 


HANS     IE3NER 
OECErtSE*    31.     1910 

EXPLANATION 

!?uSSE***3?"5AfA*?aNTs"f"SoNsriAVii 

NUHiER    JF    OATA    POINTS    AT    CONSTANT    «2 

NU.IiffJ    tip    JaTa    POINTS   AT    CONSTANT    H 

PO*eJ    FOP.    flE    XI    VALUES    HI    ME 

SUNNATION. 

POWER    FOR    T-IE    (Z    VALUCS 

3U1NATI0N. 

"W»    FIR    T-IE    13    VALUES 

llflATIQN. 

INTRJL    PARAMETER 

T    VALUES    ARE    NOT    IMCLUOEO    IN 


& 


III    THt 

IN    THE 


•1 


summation. 

"    STSTEN    MATRIX    CLCfCNTSl 


THE 


.    MCLUOEO    IN    THE 
SUNNAflON. 

(->    RIS-iT    HANO    SIOC    VECTOR    EL6N.) 
RESULT    OF    SUNNATION. 


'-"Hlk*" 


IN 


IN 


*E 


Pe 

rlsg3 

TfciE* 
rile  J 
fEiE* 

TSGER 

Tti:» 

TESER 

rue* 


Xl-COORDINATSS  OP  THE  3A7*  POINTS 
rJ-COJROlNAr^S  OF  T^i  OATi  POINTS 
XO-COOROINAfES    OF    HE    DATA    POINTS 


•COJROINATcS    OF    THE    OATi   POINTS 
.-COORDINATES    OF    HE    DATA    POINTS 
r-CJOROIMATES    OF    THE    OATA    POINTS 


9EAL 
«feiL 
■•iAL 
<EAL 


IBPORrihri     THE    3ATA    POINTS    XI.    X2j     tj    vNO    -    HAVE    TO    Is    P»0- 
VI3E0    ThUQUCH    THE    IA4CIE0   CONNQN    SLOCK    3AI»*. 


COnnOM  /  oata%  /   u.x'.o.r 

REAL    tllft  7.S  I.X2C.  7.  J  1.4  31  3.7.  5»  .7(9.  7,5  » 
S«.S 

00    07   *l«l»NPNTSl»l 

00   07    <2«1»NPNTS2»I 
30   U7    «3»I»nPNTS3»j. 

if  (ipaMRi.c4.u»   sg  ra  oi 

Sl*XUrU.*2>M)"IP0vai 
50    TO    02 

a!  If  (Ipowr2.eo.oi  go  to  93 

;2««2Hl»fW,*ii«»lPO.R2 
SO    TO   0* 
03    S2«l. 

P    (IPOwR3.EO.01    SO    TO  05 


05 
34 


07 


f3«X}(<U.«2*K3l»«lP3WR3 


;T«r(Kl.<2.<3) 

IF    (IT.Ea.U    ST.l. 


}•$»(( 

RETURN 
ENO 


S1»S2)»S3»»ST 
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f»>Wi       1PT*t       MOW  FT*    <i.«*t«9 

1  ihlclei    FUHCTUH    IEL»*     (N> 


J 


10 


J4uss  jo*o»M  !Lini**TliN  <tr«t  *4T»ii  \  *no  vector  J.     me 

/»«[*iLc    16  L*    !j    SET    H    '..     Is     rut    lLI>1In*TIJn    1$    CUnFleTeO. 
Ifcl*    LESS     I  ham    3    [NOtCiftS    4*    ERROR. 

4UTH0RI  -UNS     J?3N-» 

OftTti  j*Ct«4e*    31.    1968 

ytsusiE      Em*N»rroN  ft»e 


J  .      N  'lUnSER    jF    63U4TI3NS  I'«'65s» 

c  .     m«  cjnt*3i  '4R4iete«»  retur'ico   "a  me  iireieR 

.;  j  .  :*i.l :■•«■■»  routine. 


20 


3*6- 


.»uiv,    ■<  ."ji  ine  • 

■I    54JS3    J3R041    ElINlHtriCM    *«S 

>E«FnR*S3    C3RReCTlT. 
<0    ERROR    ElCautTcaeo    }UR  I  IG    J4USS 

jo«qan  elimination. 

6*, 3*1  JTSTFN    «4TtI,X     t  «E4l 

6*1  Rl.nf    -UNO    S13E    /ECT3K   3  Re4t, 


{  .       IN»ORTANTi     «4f8IJ    4    »N3    VECTOR    9    -utfC    r0    56    •ROVIDEO    M*iju4-t 

IS  C  .  Th£    C496LE0    C0NN3N    9t,0Crt    ..rt*.       dd  H    4    4,13     >    rfll.1 

"LL    Jt0t»40    CNA.14ES    3U*:HG    EXECUTION,     [.e.     H; 

[HtllOD. 


3C 


*1U 

EU« 


C0H1QN    /    id*      /    4.4 
*£*l    4l»*»a*lt  Jl3*> 


is  :ei»*i 

m**«i 

EFS*l.00E-O3O 
30   Oft    l«l»«»l 
30   Ji    **l#»#i 
*0  If    h.c),  t|    JO    "1    02 

C0NST*-4C*.i)/AII,I) 
30   01    W»N»i 

4<X,J)-4K,J]«C:iST«4l!,JI 

IF    li.-9.li    4K..1..J 

♦9  oi  eanrfNu* 

02    CONTINUE 

const  ■4?i«ii 

00   03    .•l.»»l 

jo  03  4(i. v  i.4i  i.4i/ca*jr 

till  I  Ml. 

ill >••«;> /const 

30   3ft    J»l»N»l 

IE«ror*o 

J?  ('Jas'tJjlJu.sr.scsi  33  ra  o» 

l£RR0a»lSRS0»»i 
Oft    CONTINUE 

if   aisHijii.;i,i>',i  30  ra  33 
09  if  <IcJ»o«.ir.N»u  so  ro  37 

Oft   CONTINUE 
40    '3    US 
07    IEl»— I 
ft9  !a.«i 

OS    RETURN 
ENO 
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1 
1 

10 

13 
:o 
i) 

30 
M 
*0 
*3 

50 
J! 

»o 
bt 

TO 

ts 

10 

39 


>6/173      0PT«2      »MON» 

subroutine  ia»  in.iter.itrnax.epslom) 


FTN     -..!•■•.*'> 


9C/12/2?.    i*. >».->■* 


.  Gaussian  elimination  *tT*  linear  uuaTIOn  srsttN  *  it-ie  »ig« 

.  h4ho  sioe  vector  3  rs  I.«*rE5a»reo  tigctAep  *ith  t*e  srsr?N  m 

.  TRIX    4    IN    T  IE    irutX    VI.       ITE9     £J    THC    -IU19ER    OF    ITERATIONS    a« 

.  solution  so  Sany  nils,  jirlt  the  'cflAlf)Tr*<  hror  :s  ijThin 

.  THE    NAGNITUOfc     3F     SfSHJN.        SHO'JLO    f-lf     SUH8ER    UF     ITERATIONS    ei 

•  *S£8   il"n**t  r*6    subroutine    :e<ii«*r£:    r-.e   Calculation   ano 

.  sits  Iter  n  a  -iejaTive  <»lje. 


Tnoast 
te« 


KONRAO    VOGELER.    -iA-«S    :E3NCR 
DECeNSER    31.    1?IQ 


.      VARIA3LE         :«»UM«1IGH 


TYPe 


TER 


TRNAX 
PSION 


NUMBER    3F    ciu»ri3NS 

nunser  jf  iterations  n  bc  retusineo 
>o  iAussu*  elimination  Ftp.FORNto 

CORRECTLT. 

«o  no  convergence  of  rue  solution  fc 

THE    3ESUE0    »iS3luri5M    tPSLON. 

iteration  3io  ri*"ia*TE  after 

itrnax  iteration  sfc/». 
naxi1un  num3er  of  iterat'tm  steps 
oesireo  resolution  of  tn£  solution 


»<6«.6*> 


SYSTEM    1ATRU     A 

RIGHT  HANQ  SIOE  /ECTOR  3 


integer 

iNTeGfcR 


INTEGER 
REAL 


REAL 

REAL 


.       IMPORT AftT  I    MATRIX    A    itO    VECTOR    1   HAV  =    T3   3c    FROVIOEO    THROUGH 
•  Trig    LA3ELE3   COnnjn    HOCK    ABJ.      SOTh    A    a«j    i    .ill 

.  *ill  jn0er4j  changes  ourihg  execution,  i.e.  thc 

.  Slimution. 

• 

*25SA5N*/*A8i*/*i* **************************** ******** •••••• 

REAL  vta<..o5),i(i}).Al»«.4»(.aibJ) 

LOGICAL     IPRlNl 

DATA  [PRINT  /.FALSE./ 

OAlA  PI  /3.1<.U-)'.b3'./ 


FORMATS    IEI*    START 
601    FORMAT    <///}<a-i      iaTRU 


JNO    VEKTOR   3    VOR    ELININATI3N/) 

ERATIONSSCHRITT/ 


602    FORMAT  M *i»l3(aX. ;; I  I  I 

m  massif  i^ji^'n^sKiiciUs'jr-Ti^ffi.-nriwiT 
•• 

609  FORMAT  I///2X. 12, 2*H    «3EF= I2IENTEN    AUS    GAUSS/) 
606    FORMAT  (l2X»i2/lX»;PE4.2)) 

bat   FORMAT  I///39H      <LtMSTER    rcUeNVEKTnprCMNITTvrNilEL    •.E12.3) 

6d«    FURMaT  (27X.124      <E3I0IU^   •.Eli.J.iH   •UCH»I3»iin    I TERArTjNeN) 

40<»    FORMAT  1///20X.11-I      -i'll.'-*    <(»(2tlrtt;2/3H>     •tei2.j/l 

610  FORMAT  I///2JH      *TSTE*   MAC1*   E tl^lM ATlJft/ ) 
FORMATS    IEL*    STOP 

10*6 
N«M*1 
•      «»M-1 


01 


REARRANGE    VECTOR    3    IH    ARRAT    V*    IF    N    IS    LESS    THAN    66. 
IF    iN.EQ.6*)    GO   TO   02 

00    il     I«I>M>1 

VU.MMVI  I»«>l 

PRINT    STSTEM    MiTRIX    A    AN8   »;JHT    HANO    JIOE    VECTOR    3    3EF0RE     THE 
ELIMINATION,    IF   IArIHT    IS    Scf    .TRUE.    . 


92  \t  wnwi  rsh 


(L9>    601) 
.     (U»    )02)     !J.J«l»M»l) 

OJ    IF    flF*INt)"*RtTE     (L0»    603)     C»<  VI  t»  J)  #  J«l#«#l) 
IF    (IPRIND    vRITE    (L3>    602)    N»4«1*M«1> 

LINEAR    EQUATION  STSTErl      V(M.M  )»6(N  l»V(N»N»l) 

38  8:  SriiSli 

06    Ail,j)«v«  Iui 
00   03    t«l«N»l 
03    3(I)«.0 

CHICK    FOR    STAtUITY. 
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90 


SUMOUTInE    let*  /»/l7t       JFW      ,>11N»  fT«    *.a»t.a  3C/i2/c-<. 

FI»,0 

06    -llM.i 

30   94    [•1»N*1 
00  07    J»1#n.I 

00  ;o   :>i»«»i 


99 

00    10    J'JJ.N.i 

5x«.o 

00   09    I3«i»«.l 
i)9    S*-$**yit »  I3i«y(j»{  Ji 
100  e.s*/  (in  )»G(j  )  ) 

If     It8S  If  I  .iHUSIFIl  )    ei-f 
10    CONTINUE 
f l.»*S(f I  ) 

fi.acosifij 

109  f  l-lf  I'laO.J.'^l 

jt*«t  iteration  rs  fino  ».io  inmove  :jtuTinH. 

00    27    LX»l»999,l 
UO  ITS*-«.X 

00    13    J»1>»»1 

mvotisation  a.io  saussian  ST6>  foa  cntumt  j. 

113  |X»4 

e-vtJ.JI 
00   12    .•■>,", I 

If    (*»i(v<l»JI  >.ST.A9S(E>>    10  TO   11 
ia  ro  12 
no  n  1 1«  i 

tiXilril 

12  Continue 

IF    i  tl.fcJ.JI    SO   TO    l» 

oo  ii  ;>w»i 

129  iltl<«U<>I) 
y<i»»  i)«yij»n 

13  /Uili-Jili 
SHI«.0 

19   nnaj*i 

130  00  13    >nn»n,i 
»•«( i.ji/vi j.j  i 

OO    13     *»»J»« 

Vlt«ftftl>VtXfftKW(J*<M>*U 

19    CONTINUE 


139 


160 


If    I y(I#I).ia..4J     S3    TO    12 
SOLUTION* 


190  00    16    !-l,N.i 

16   ;;:)•. o 

(I»H»I».«I/»IN|1I 

00    19    [W,«.l 

I«N-U 

199  .     If     lV(I«l).E0..OI    iO    TO    12 

Sll)»VU»fll 

00   17    ij-l»«,l 

IF    (J.ES.I)    SO   n    19 
130  IT    ill  ].«.(  i)-SUi«y(  I.J) 

19    CONTINUE 
19  III  )  •-.(  II  /y(  i#  t) 

CnjUECT    iOLUTIOK    }T    AN3UNT    IF    =J»0«.    <H*--i    yAS     iNyOtVEO    W    TJJS 

199  c  m»»riQN  stef.    ?a<e  This  *ect:»  a*  JtW.t  vector  ?3r  fne  -mt 

C  ITERATION. 


C 


30   20    t-l.N.i 
20    ill  I- il Ii •Si!) 

COMPUTE    THE    »£N»ININ5    ERROR. 


00    22    I»1.N,1 

y<l«A>«.a 

169  30   21    J'l.N.l 

51    »Tj#a»«<MI*«»»*UI«»CX»4l 

11      KI..1IMI  !,NI-yC»N) 

»••  J 

30  :i   i*i>n,i 

170  23    4««*v(!.N>*y(t,f1| 


177 


susrjuuns  iet*  7»/i7s     ji»t^    ^ow  ftm  w. *♦*»«!  :c/i2/2- 

••SS«T(«( 

IF    (fUkT.EPSLONI    SO    TO    2S 

if  (Li.sr.n»n»x)   si  rn  3i 

so  2*   i»i»H,i 

179  2*    SII>««tl> 

print  systen  n»r*rx  a  *no  riont  hano  sidc  vector  3  »ft£»   rwe 
c  pivatiS»f ion. 

1*0  tf     (  I  P*  [NT  •     WRITE.     (L3»     40*1     IX 

IF    (IPRlNT)    .*!rj     llOi    902)     (J.-'l.H.i) 

oa  2S   :»i»n»i 

29    IF    (IFklNT)    WRITE     a).     503)     t  .  (  V  (  I  .  J  I  ,  J  •:»  1.  1 )  #0  (I  ) 
IF    (IminT)    ■«[[;    UJ»    4y2>    <J.J«ltN.l) 


1(9  00    26    l«i»».l 

03   26    J«l.N.[ 
26    <(I>J  l««l  I»Ji 

C 

190  C 


200 


209 


STO>    (.30*    TO   FINO    S31UTI0N. 

27  conrtNuf 

REARRA.iGS    SOUTH*    II    THE    mir    4, 


19!  28    CONTINUE 

oa  29  I'ui.i 


29  V(1<61I««(I) 
PRINT    SOLUTION.    IF    IPRINT    tS    SET    .TRUE.    . 

30  IF    ({MINT]    WRITE    U3»    403  1    i 

if  (ifrimt)  writc    iui  404i   (:»•;!  !)»:•:»«.:) 

IF    (IfRINT)    wRITe     ((.I*    407)    rl 
IF    (IPRINT)    .KITS     (LJ.     40a)    >.(.« 

ORUCKEK    9ER    "ATRH    *    OHO    )ES    VECTORS    1    HACH   EUH IN4T 131. 

IF    (IFRINT)    WRITE    (L3»    613) 


IF     (IFRINT)     WRITE      (LJ»     402)      (  J.J«l,.1.i) 

210  00    31    ['li*.l 

31    IF    (CFRIMT)    .R{T;     (L3.    40))     I t < v ( t . J ) . J •:/ i. 1) 

IF    (IFRINT)     wit  IT?     IlJ.    402)     (JU'li.tlii 
RETURN 

219  J  ERROR    ENCOUNTERED.       ST3». 


32 


220  RETURN 

C 

C 


ITER.-l 

WRjTt    (13.    4C9)    I.U'HiH 

S«I)«5U)/w(I»  I) 


13    ITER  —  ITER 
229  RETURN 


NO    CONVERGENCE    OF     THE     S01UTI3N. 

TER 

ETl 

NO 
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i 
* 


7*/17j       19T»l       ?nonp  FT*    S.9**ii  iU/ii/Z^.     if. 

X  DEAL     FUNCTION    F*     (  iZi  F.  «  I.  I  2  .  <  J.  M3R08  l«  fHKOP^,  itJ90«  3  I 

1      CONFUTE    FUNCTION    yALUe     r  •;  <  «1 ,  I  2,  »  3  I    CiF    4    fcO  — '  JL  YNJNJAC    'JSItG    I 

j  !  .     t-u  -hj»»t«   scheme  l.t  jttoEi  ra  inctca;;  :ntj  »«o  *c:j*act. 

C  •      AUTMO«l  HAHS    I£»NE«  * 

C  .       OiTci  JECENSEU     11/     i'80 

10  C  I     »»«Ullt       EXFtANATION  rw«  ! 

J         :    MUM      M£fAMUCaEMMK!l3K Hi*.*"  I 

C  .       »1  Xl-VALJE,    *Hc9£    T    *ILL    3C    C>1.1i»UTtO.         9E4L 

j  .       «2  «2-y4CJE.    K-teRE    r    <JLL    dl    C3"l»UT?d.         »e»l 

IS  C  .      li  «I-y»LUE»    «-t£»E    r    *iu    d?    C3rtPUTS0.         ietl, 

C  .     iq«o«i  jaoer    if   rug   •>jt»t3«UL  In  t  it  i«t-;cj 

C  .  U    3I»ECT'Qi. 

C  .  N0I082  3»0*«     IF    frt«    >3Lr>*0«[AL    :*l   THE  ixTiit* 

C  .  «Z    OIXcCTIJ*. 

20  :  .  SOUOO  ;«OcJ     JF     ME    ^OLTNOntAC    :rt    THE  INTeJE* 

c  .  x3  ji^Ecna-*. 

2  .       FA  1B0INATE    OF    THE    *5-»0crH01IAt.  *«AC 

£      • : 

IS  REAL    cSlFift*»f5(4*l»El64l 

NCQfcFj.«»0f>0i»i»i 

NCQSF3-t»O»0«3»i 

10  CONFUTE    COEFFICIENTS    9    JASEQ    On    COEFFICIENTS    COEF    4N0    13 

30  02    ll«l#NCO§Fl#l, 

lC»Ul-tMNCOeF2«tCJSF3*I2«N:3EF3 

T«CQEF(iC) 

00   i»     rJl»l,MJ*0».J.l 
I3-NCOEF3-131 

IC«(Jl-ll*'"CaEFJ«HC3EF3»«I2-l>«MCC)EFi*:3 
40  91     fCOeFttd  i«ij»T 

02    DUOI-T 
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